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Figuur voorpagina: Schematische voorstelling van de gebeurtenissen tijdens een 
zonnevlam. se | 


De eerste aflevering van het NTvN-A in 1980, dat thans voor u ligt, begint 
met een korte bijdrage van het bestuur van de NNV, waarin een recente re- 
dactie-wijziging wordt gemeld. Daarin worden ook de verdiensten van de 
aftredende redactieleden voor dit tijdschrift gememoreerd. Zodoende kan 
ik mij hier beperken tot de opmerking dat m.i. de oude redactie er volledig 
in is geslaagd om de enkele jaren geleden nieuw geformuleerde doelstelling 
van het NTvN-A — het vele interessante speurwerk dat in Nederland of 
elders gebeurt op natuurkundig of verwant gebied in korte overzichtsar- 
tikelen onder de aandacht van de leden van de NNV te brengen — te rea- 
liseren. Het was verheugend dat zoveel onderzoekers bereid zijn geweest 
voor het NTvN-A over hun werk te schrijven. leder die wel eens voor de 
taak heeft gestaan om de essentie van een verricht onderzoek duidelijk 
weer te geven voor een lezerskring zoals die van dit tijdschrift met een 
wijd uiteenlopende fysische belangstelling, weet dat dit geen eenvoudige 
opgave is. Het vereist misschien wel meer inspanning en inzicht dan het 
schrijven van een artikel voor een vaktijdschrift. De redactie hoopt ook in 
de toekomst op dezelfde enthousiaste medewerking te kunnen rekenen. 
Dit nummer bevat een zestal artikelen op het gebied van de halfgeleiderfy- 
sica en -technologie, een in eerste instantie doelgericht onderzoeksgebied 
met vele interessante fundamentele en praktische problemen, zoals u bij 
lezing van de artikelen duidelijk zal worden. Deze artikelen zijn tot stand 
gekomen met de medewerking van de sectie Halfgeleiderfysica en -elek- 
tronica van de NNV. 

De overige artikelen zijn gewijd aan een viertal onderwerpen uit de astro- 
fysica en geven een duidelijke illustratie van hoe boeiend dit gebied van 
onderzoek is. 


A. van der Woude 








De redactie van het NTvN-A heeft enkele veranderingen ondergaan. Dr. J. 
Polman heeft afscheid genomen als redactievoorzitter. Hij wordt als zo- 
danig opgevolgd door dr. A. van der Woude, KVI, Rijksuniversiteit Gro- 
ningen. Ook dr. P.W.M. Glaudemans heeft zijn taak als lid van de redactie 
neergelegd. Als nieuwe leden treden toe dr. A. Broese van Groenou, Philips 
Natuurkundig Laboratorium, Eindhoven en dr. H.M. Rietveld, ECN, Pet- 
ten, waardoor de redactie dus met één lid is uitgebreid. Bij deze gelegenheid 
past een woord van oprechte dank van het NNV-bestuur aan Glaudemans 
en Polman. 

Piet Glaudemnmans, 48, hoogleraar in de kernfysica aan de Rijksuniversiteit 
te Utrecht, trad in oktober 1970 toe tot de toenmalige redactieraad. Door 
problemen van verschillende aard werd het in het midden van de zeventiger 
jaren steeds moeilijker bij de toenmalige opzet van het tijdschrift de gewen- 
ste kwaliteit te handhaven. Tegelijkertijd begon de groei van het ledental 
van de vereniging te stagneren, ja zelfs trad er een lichte daling op (zie 
NNV-jaarboek 1978/79, p. 8). Een samenhang met de moeilijkheden van 
het tijdschrift werd vermoed. Glaudemans, in 1974 secretaris van de NNV 
geworden, heeft zich toen met hart en ziel ingespannen voor een verbete- 
ring in de situatie, hetgeen in januari 1977 heeft geleid tot de ons nu ver- 
trouwde splitsing van het tijdschrift in de delen A en B. Hij werd lid van 
de redactie van het NTvN-A. Het goede onthaal van dit blad bij de leden 
en de snelle stijging van het ledental in de jaren 1977 en 1978 hebben hem 
terecht veel voldoening gegeven. 

Jan Polman, 44, verbonden aan het Philips Natuurkundig Laboratorium te 
Eindhoven, werd in januari 1972 lid van de redactieraad en twee jaar later 
algemeen redacteur. Ook hij was nauw betrokken bij het overleg over de 
reorganisatie. Zijn brede belangstelling voor alle aspecten van de natuur- 
kunde maakte hem bij uitstek geschikt voor het redactievoorzitterschap. 
In deze functie, die hij in januari 1977 aanvaardde, is hij erin geslaagd een 
mooi evenwicht te vinden tussen wetenschappelijke strengheid en popula- 
risering, tussen aandacht voor fundamenteel onderzoek en voor technische | 
toepassingen, terwijl daarnaast ook het onderwijs en de problemen van na- 
tuurkunde en samenleving niet werden vergeten. 

Glaudemans en Polman hebben in de afgelopen jaren in belangrijke mate 
het gezicht van het NTvN bepaald. Hun werk wordt nu overgenomen 
door anderen. Het bestuur prijst zich zeer gelukkig met de nieuwe beman- 
ning en wenst haar veel succes. Het doet tevens een beroep op de leden de 
taak van de redactie te verlichten door een enthousiaste en spontane me- 
dewerking. 


O.H.F. Vrehen 
voorzitter redactie-adviesraad 


Lithografie ten behoeve van 
Very Large Scale Integration 


R.P. Kramer 


Het toenemen van de complexiteit van geintegreerde schakelingen gaat enerzijds ge- 
paard met het groter worden van het chip-oppervlak, anderzijds met het kleiner wor- 
den van de minimale afmetingen op de chip. De eisen die worden gesteld aan de pa- 
troondefinitie worden daardoor strenger. In dit artikel wordt een overzicht gegeven 
van verschillende lithografische methoden die ons in deze ter beschikking staan; spe- 
ciale aandacht zal worden besteed aan de licht-optische repeteerprojector en de elek- 


tronenoptische beeldprojector. 


De vervaardiging van een transistor of ge- 
integreerde schakeling bestaat uit de op- 
eenvolging van een aantal lithografische 
stappen, afgewisseld door een combinatie 
van etsen, oxyderen, implanteren, diffun- 
deren en deponeren van diëlektrische of 
geleidende laagjes. Onder lithografie ver- 
staan we de verzameling processtappen 
die ertoe leidt dat een gewenst patroon 
op het oppervlak van een substraat wordt 
overgebracht. In de IC(Integrated Cir- 
cuits)-technologie gebeurde dat tot voor 
kort vrijwel uitsluitend door een masker, 
met daarop het gewenste patroon gerepe- 
teerd aanwezig, in nauw contact te bren- 
gen met een plak silicium waarop een 
laagje fotogevoelige lak was aangebracht 
en deze plak-maskercombinatie met ultra- 
violet licht te belichten (het zogenaamde 
contact-afdrukken). Na ontwikkeling van 
de laklaag is het maskerpatroon op de 
plak overgebracht. Aangezien voor de ver- 
vaardiging van een transistor of IC meer- 
dere patronen na elkaar moeten worden 
aangebracht, moeten de verschillende 
maskers ten opzichte van elkaar worden 
uitgericht. 


Aan dit contact-afdrukken kleeft een aan- 
tal bezwaren zoals maskerbeschadiging 
(ten gevolge van het op elkaar drukken 
van masker en plak kunnen maskers in de 
produktie slechts een keer of vijf worden 
gebruikt) en buigingsverschijnselen (ten 
gevolge van het niet overal in perfect con- 
tact met elkaar zijn van masker en plak). 
Bovendien vereist de vervaardiging van 
het masker vele processtappen. 


Hoewel er geen principieel verschil be- 
staat tussen de vervaardiging van (bijvoor- 
beeld) een enkele transistor en (bijvoor- 
beeld) een microprocessor, worden bo- 
vengenoemde bezwaren groter als we naar 
meer complexe circuits gaan. Deze com- 
plexiteit komt namelijk, voor wat de 
technologie betreft, vooral in de lithogra- 
fie tot uitdrukking, en wel door het ge- 
lijktijdig optreden van alsmaar kleiner 
wordende minimale afmetingen (2 à 3 
um), de daarbij behorende grotere uit- 
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_ richtnauwkeurigheidseisen en het grotere 


chip-oppervlak (> 5 x 5 mm? ). De uit- 
richtnauwkeurigheid moet groter worden 
omdat het nuttig effect van kleinere af°- 
metingen goeddeels verloren gaat als deze 
niet gepaard gaan met een nauwkeurige 
passing van de verschillende lagen. Als 
vuistregel geldt dat de uitrichtnauwkeu- 
righeid beter dan een kwart van de klein- 
ste afmeting moet zijn. 


De categorie circuits met bovengenoemde 
complexiteit en met meer dan 10° tran- 
sistoren wordt Very Large Scale Integra- 


Fig. 1. Schematische voorstelling van 
verschillende lithografische methoden. De 
magneetband boven in de figuur bevat in 
digitale vorm de geometrische gegevens 
van het ontwerp. Aangegeven is op welke 
manieren vanaf de magneetband het ont- 
werp op de Si-plak wordt afgebeeld. 








Drs. R.P. Kramer (36) studeerde experi- 
mentele natuurkunde aan de Gemeente 
Universiteit van Amsterdam. Hij is sinds 
1970 werkzaam als wetenschappelijk me- 
dewerker van het Philips Natuurkundig 
Laboratorium te Eindhoven met als vak- 
gebied halfgeleiderfysica en -technologie. 
De laatste jaren houdt hij zich vooral be- 
zig met problemen die samenhangen met 
het verkleinen van MOS-schakelingen, in 
het bijzonder met de lithografie hiervan. 


tion (VLSI) genoemd. De combinatie van 
bovengenoemde bezwaren maakt het 
steeds moeilijker om met contactafdruk- 
ken te voldoen aan de eisen die aan de 
vervaardiging van VLSI-circuits worden 
gesteld, en om die reden wordt er druk 
gewerkt aan een aantal alternatieve me- 
thoden. Een (onvolledig) overzicht van de 
verschillende mogelijkheden treft U aan 
in figuur 1. In de loop van dit verhaal zul- 
len verschillende begrippen uit deze fi- 
guur nader worden toegelicht. De aard 
van de oplossingen waaraan wordt ge- 
werkt is soms praktisch (geen contact, 
minder vervaardigingsstappen), soms fun- 
damenteel (kleinere golflengten) en soms 
beide. 
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Het is in dit bestek niet doenlijk om al 
deze verschillende methoden op hun me- 
rites te beoordelen; wel is het nuttig om 
een aantal criteria te noemen die aan deze 
beoordeling ten grondslag liggen: 1) 
kleinst af te beelden afmeting; 2) uitricht- 
nauwkeurigheid; 3) repeteermogelijkheid 
(in verband met plakvervorming); 4) 
automatiseringsaspecten; 5) golf- en 
deeltjesverschijnselen (staande golven, in- 
terferentie, verstrooiing); 6) focusdiepte ; 
7) uniformiteit van belichting; 8) moge- 
lijkheid om voor plakvervorming te corri- 
geren;9) invloed van plakgrootte; 10) 
onderlinge uitwisselbaarheid van de ma- 
chines; 1 1) maskervervaardigingsproble- 
men en 12) aantal te belichten plakken 
per uur ofwel: kosten! De meeste van de- 
ze punten zullen in de loop van dit ver- 
haal nog ter sprake komen. 


Als voorbeeld van een methode die dui- 
delijke vooruitgang heeft gebracht, zou- 
den we licht-optische 1: 1-beeldprojectie 
(lens of spiegel) willen noemen. Deze 
heeft als voornaamste voordeel dat mas- 
ker en plak niet meer met elkaar in con- 
tact komen, waardoor geen maskerbe- 
schadiging meer optreedt. Met deze me- 
thode kunnen 4-inch-plakken worden be- 
licht tot in het gebied van 2 à 3 um mini- 
male afmetingen. Mede vanwege de voor- 
treffelijke werking van dit type apparaat 
is deze methode binnen enkele jaren de 
industrie-standaard voor de produktie van 
geavanceerde circuits geworden. 


In plaats van met optische technieken 
wordt tegenwoordig meer en meer het 
masker vervaardigd door met een elektro- 
nenbundel het gewenste patroon in elek- 
trongevoelige lak te schrijven. Deze bun- 
del kan een zeer kleine diameter (onge- 
veer 0,1 um) hebben en patronen kunnen 
zeer nauwkeurig worden aangebracht. De 
combinatie van elektronenstraalmaskers 
en optische 1:1-beeldprojectie is dan ook 
een zeer aantrekkelijke methode. Toch 
heeft ook deze methode zijn beperkingen 
(minimale afmetingen, focusdiepte) en we 
willen daarom twee andere methoden na- 
der toelichten, waarvan de eerste zeker, 
en de andere misschien, een belangrijk 
deel van de lithografie ten behoeve van 
VLSI voor zijn rekening zal nemen. 


De Repeteerprojector 

Thans wordt het gerepeteerde masker, in- 
dien door optische technieken verkregen, 
vervaardigd door een 5 of 10 maal ver- 
groot circuitpatroon (het reticle) te ver- 
kleinen en te repeteren op het uiteinde- 
lijke (moeder-)masker. Het ligt nogal voor 
de hand om dit uiteindelijke masker over 
te slaan en het reticle direct op de Si-plak 
te projecteren en te repeteren. Dit nu is 
het principe van de repeteerprojector (zie 
figuur 1). Dat deze methode pas sedert 
korte tijd wordt toegepast heeft te maken 
met de verschillen die optreden wanneer 
we in plaats van een masker een Si-plak 
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belichten. Het voornaamste verschil be- 
treft de snelheid, waarmee een gehele 
plak moet kunnen worden belicht om de 
methode ook economisch te maken. Ook 
het feit dat de plak een minder ideaal 
beeldvlak vormt dan het masker zorgt 
voor complicaties. Bovendien kan de 
vorm van een Si-plak in de loop van het 
vervaardigingsproces van een geiïntegreer- 
de schakeling veranderen; hierop komen 
we later nog terug. 


De volgende ontwikkelingen hebben het 
gebruik van repeteerprojectoren in massa- 
produktie binnen bereik gebracht. 


1) Lenzen met hoog scheidend vermo- 
gen. 

De aan de lenzen te stellen eisen zijn zeer 
zwaar voor wat betreft de combinatie van 
scheidend vermogen, veldgrootte en 
beeldvervorming. Een lens van topkwali- 
teit van dit moment heeft een Numerieke 
Apertuur van 0,3 (minimum dimensie 1 à 
1,5 um) bij een beeldveld van 10 x 10 
mm? en een beeldvervorming van onge- 
veer 0,1 um. 


2) Sterke en homogene lichtbronnen. 
Er zijn nu lichtbronnen die in repeteer- 
projectoren een belichtingstijd van % tot 
\ sec geven bij een homogeniteit van + 
2%. 


3) Snelle, met behulp van laser-interfe- 
rometersystemen gecontroleerde, x-y-po- 
sitionering. 

De precisie van deze positionering ligt in 
de orde van £ 0,1 um en de verplaats- 
ingstijden liggen in de orde van % tot % 
seconde. 


4) Autofocus. 

Omdat het oppervlak van een Si-plak over 
vele microns in de hoogte kan variëren is 
het noodzakelijk dat de projectiekolom 
voor iedere belichting wordt gefocus- 
seerd. Er zijn hiervoor verschillende con- 
cepten ontwikkeld. 


5) Uitrichtsystemen. 

Het is in het licht van bovengenoemde 
minimale afmetingen en eerdere opmer- 
kingen over vereiste uitrichtnauwkeurig- 
heid duidelijk dat een adequaat uitricht- 
systeem van het hoogste belang is en dat 
een automatisch uitrichtsysteem welhaast 
noodzakelijk is. Het enige automatische 
uitrichtsysteem dat tot nu toe met succes 
is toegepast is dat van de repeteerprojec- 
tor die in het Philips Natuurkundig Labo- 
ratorium is ontwikkeld. 


Samenvattend zouden wij de volgende 
gunstige eigenschappen van de repeteer- 
projector willen benadrukken: 1. Als 
masker wordt een reticle gebruikt. Voor 
de vervaardiging hiervan zijn slechts één 
of twee stappen nodig en dit zal zeker op- 
brengstverhogend werken; bovendien is 
dit reticle eenvoudig te inspecteren en te 


repareren. 2. Het scheidend vermogen ligt 
in de orde van 1 tot 1,5 um. 3. Het sys- 
teem is in staat zich aan de vorm van de 
plak aan te passen. Dit wil zeggen dat de 
grootte van de plak in eerste benadering 
niet relevant is en dat voor lokale hoogte- 
variaties kan worden gecorrigeerd. Boven- 
dien kan bij toepassing van een adequaat 
uitrichtsysteem voor laterale plakver- 
vorming (waarover later) worden gecorri- 
geerd, hetzij door de grootte van de repe- 
teerstap te variëren, hetzij door op meer- 
dere plaatsen uit te richten. 


Het is vanwege de combinatie van deze 
voordelen dat wij ervan overtuigd zijn dat 
er voor de repeteerprojector een zeer 
ruime plaats zal zijn ingeruimd in het 
VLSI-tijdperk. 


De Silicon-Repeater 

Als voorbeeld van een repeteerprojector, 
zouden wij enige aandacht willen beste- 
den aan de Philips repeteerprojector, de 
Silicon-Repeater (Si-Re). _ 

Figuur 2 toont een schematische voorstel- 
ling van de Si-Re, dit apparaat werkt met 
een verkleiningsfactor van 5. 


Er is uitgegaan van twee hydrostatisch ge- 
lagerde sleden die door hydraulische li- 
neaire motoren langs onderling lood- 
rechte banen bewogen kunnen worden. 
De ene slede draagt, op een speciale in- 
stelbare houder, de met fotolak bedekte 
Si-plak, de andere draagt de projectieko- 
lom. Twee laser-interferometersystemen 
maken het mogelijk om de positie van 
beide sleden met een nauwkeurigheid be- 
ter dan £ 0,1 um te meten. De projectie- 
kolom bevat als lichtbron een superhoge- 
drukkwiklamp. Een elliptische spiegel 
achter de lamp en een speciaal condensor- 
systeem zorgen voor een voldoend sterke 
en homogene verlichting van het te pro- 
jecteren masker. Het lenzensysteem is te- 
lecentrisch, dat wil zeggen dat bij projec- 
teren op een oneffen beeldveld de vergro- 
ting niet varieert. Een bijzonder sterk 
punt van deze Si-Re is het automatische 
uitrichtsysteem, dat inmiddels is beschre- 
ven! . Ten behoeve van dit uitrichtsys- 
teem bevatten zowel masker als plak twee 
uitrichtkenmerken bestaande uit rasters 
in de x- en y-richting. De rasters op de 
plak worden door een laserbundel (links 
in de figuur) belicht en door de projectie- 
lens afgebeeld op de overeenkomstige ras- 
ters op het masker. Voor verdere details 
zij naar genoemde publikatie verwezen. 
Met behulp van de informatie die tijdens 
de automatische uitrichtcyclus wordt ver- 
kregen is het in principe mogelijk om de 
door het laser-interferometer-systeem ge- 
controleerde repeteerstap-grootte auto- 
matisch aan te passen en zo voor even- 
tuele plakvervorming te corrigeren. Als 
bijvoorbeeld de plak groter is geworden 
gedurende het IC-vervaardigingsproces, 
dan moet een grotere slede-verplaatsing 
worden uitgevoerd. Het is ook mogelijk 
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Fig. 2. Schema van de Silicon-Repeater. 





om meerdere uitrichtkenmerken aan te 
brengen, en dus navenant vaker uit te 
richten. De uitrichtnauwkeurigheid van 
dit systeem is in de orde van £ 0,1 um. 
In de praktijk vinden we een 30-waarde 
van £ 0,25 um. We hebben nu enige jaren 
ervaring opgedaan met dit apparaat en op 
grond van deze ervaring wordt nu gewerkt 
aan een tweede generatie, met als voor- 
naamste verbetering een hogere produk- 
tiesnelheid. Nog enkele getallen, die een 
indruk geven van de prestaties van de hui- 
dige Si-Re: maximaal maskerbeeld: 





Fig. 3. Foto, waarop de automatische 
uitrichtnauwkeurigheid van de Silicon-Re- 
peater te zien valt. Het verschil in grootte 
tussen twee op elkaar passende vierkant- 
jes, die in twee verschillende belichtingen 
zijn aangebracht, is steeds 1 à 2 um. Het 
betreft hier een 5-masker-proces. 
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| 50 x 50 mm? ; beeldveld bij maximale 


scherpte: 7 x 7 mm? ; kleinste detail in 
schakeling: 1,5 tot 2 um; belichtingstijd 
per beeldveld 0,5 s; verplaatsingstijd: 0,5 
seconde. 


Figuur 3 laat een voorbeeld zien van de 
uitrichtnauwkeurigheid die wordt behaald 
bij de vervaardiging van een IC. De paars- 
gewijze in elkaar passende vierkantjes zijn 
hulpfiguren die met de verschillende mas- 
kerlagen op elkaar worden gelegd en 
bieden de mogelijkheid om de uitricht- 
nauwkeurigheid visueel te controleren; 
het verschil in grootte tussen twee op el- 
kaar passende vierkantjes is l à 2 um. 


De elektronenoptische beeldprojector® 


Principe 

De methode berust erop dat het patroon 
van een in vacuüm geplaatst masker door 
middel van elektronen via l:1-projectie 
op (bijvoorbeeld) een siliciumplak wordt 
overgebracht. Hiertoe wordt een chroom- 
masker, in essentie identiek aan een 
masker zoals dat voor conventionele fo- 
to-lithografie wordt gebruikt, bedekt met 
een fotokäathode (ongeveer 200 Á cesium- 
jodide). Dit masker wordt van achteren 
beschenen met ultraviolet licht, zodat 
elektronen worden geëmitteerd in al die 
gebieden die niet door Cr zijn gemas- 
keerd. Deze elektronen worden door een 
elektrisch veld versneld naar een met een 
elektrongevoelige laklaag (PMMA, poly- 
methylmetacrylaat) bedekte Si-plak en 
gefocusseerd door een magnetisch veld. 
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Fig. 4. Schema van de elektronenop- 
tische beeldprojector. 


Omdat we met UV licht van 184,9 nm 
werken, moet het maskersubstraatmateri- 
aal kwarts zijn. Het uitrichten van de op- 
eenvolgende lagen gebeurt automatisch 
met een uitrichtnauwkeurigheid van + 0,1 
um. Translatiefouten worden gecorrigeerd 
door het aanpassen van twee, als afbuig- 
spoelen fungerende, hulpmagneetvelden, 
die loodrecht op elkaar en op het focus- 
seerveld staan; rotatiefouten worden 
mechanisch geëlimineerd door het masker 
te roteren. In figuur 4 is een schematische 
doorsnede van de projector te zien. 


Hierna wordt op enkele deelaspecten van 
dit systeem nader ingegaan. 


Fotokathode 

Csl is in een elektronenoptische beeld- 
projector een geschikte fotokathode. De- 
ze stof voldoet aan alle hoofdeisen voor 
gebruik in een machine voor massapro- 
duktie: eenvoudige bereiding en verwijde- 
ring, stabiel gedrag in droge lucht, be- 
trouwbare werking in een betrekkelijk 
slecht vacuüm (10 mPa) en een stabiele 
fotostroom met een stroomdichtheid van 
tenminste 1 UA/cm*. 


De combinatie van deze gunstige eigen- 
schappen zorgt ervoor dat de laag Csl in 
de machine zelf aangebracht kan worden, 
ongeveer 50 belichtingen kan verrichten 
en bij gebruik van PMMA (gevoeligheid 
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ongeveer 100 uC/cm?) een belichtingstijd 
in de orde van 10 seconden oplevert. Het 
masker zelf lijdt niet onder het gebruik en 
heeft, in principe, het eeuwige leven. 


De gebruikte combinatie van UV fotonen 
(lagedrukkwiklamp met een golflengte 
van 184,9 nm, foto-energie = 6,7 eV) en 
Csl als fotokathode (uittree-arbeid = 6,5 
eV) levert een belangrijk voordeel op. Dit 
voordeel is dat de energie van de geëmit- 
teerde elektronen maximaal ongeveer 0,5 
eV bedraagt. Hierdoor is er slechts een ge- 
ringe spreiding in deze beginenergie en 
hebben we in het beeld geen last van het 
elektronenoptisch equivalent van chro-_ 
matische aberratie en dus van verlies aan 
scherpte en scherpte-diepte. 


Beeldvorming, scheidend vermogen en 
beeldvervorming 

De elektronen, die door de fotokathode 
worden geëmitteerd, worden naar de Si- 
plak (de anode!) toe versneld door een 
elektrische spanning van ongeveer 20 kV 
en gefocusseerd door een magnetisch veld 
van ongeveer 80 kA/m. De plaats van het 
focus hangt af van de sterkte van het 
magnetisch veld; teneinde eventuele ver- 
tekeningen in het beeld kleiner dan 0,1 
um te houden, mogen afwijkingen in de 
homogeniteit van het veld ten hoogste 
van de orde van grootte van 10” zijn, 
over het gehele oppervlak van de plak. 
Met behulp van een samengestelde sole- 
noïde is een magnetisch veld verkregen, 
waarvan de homogeniteit over een gebied 
van 50 mm diameter beter dan 3 op 10° 
is. 

De elektronen worden niet alleen door 
het magnetisch veld gefocusseerd, maar 
in feite door de combinatie van het mag- 
netisch veld en het elektrisch veld. Het 
focusseren van het systeem geschiedt in 
de praktijk dan ook door het veranderen 
van de versnelspanning. De combinatie 
van elektrisch en magnetisch veld moet 
ervoor zorgen dat de elektronen, als ze de 


afstand masker-plak (in dit geval 16 mm) 


hebben afgelegd, precies een hele omwen- 
teling rond de magnetische veldlijnen heb- 
ben gemaakt. 


Als de machine in focus is, is aan de 
volgende formule voldaan: 


cos « VVoleV —1 +(m/Bd) 
* V2mV/e =0 


waarin: « _ = uittree-hoek van de ge- 
emitteerde elektronen 
Vo _= uittree-energie (onge- 


veer 0, 2 eV) 

= elementaire lading 

= versnelspanning 

= magnetische inductie 

= afstand tussen anode en 
kathode 

= massa van het elektron. 


[all 


= 
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Het scheidend vermogen van dit systeem 
wordt niet zozeer bepaald door elektro- 
nenoptische factoren (homogeniteit mag- 
neetveld, uittree-energie van de elektro- 
nen — als we ons zouden beperken tot de- 
ze zuiver elektronenoptische overweging- 
en, dan zou dit leiden tot een scheidend 
vermogen van 0,01 um!), als wel door 
factoren als lakdikte, onvlakheid van de 
plak, elektronenverstrooiing in lak en sub- 
straat, dikte van het Csl, randscherpte van 
het Cr-masker enzovoort. Hoewel het niet 
mogelijk is om een exact getal te geven 
voor het gezamenlijke effect van deze fac- 
toren, is het uit ervaring duidelijk gewor- 
den dat het scheidend vermogen van de 
‘opvallende’ elektronen in de orde van 
enkele tienden van microns ligt. Eén van 
de gunstige aspecten van dit systeem is de 
grote scherptediepte: als we een randon- 
scherpte toelaten van 0,1 um dan is die 
scherptediepte 100 um en dit betekent 
onder andere dat de onvlakheid van de 
plak ons in eerste benadering niet erg 
hindert. 


De term ‘opvallende’ elektronen is hier- 
boven met opzet gekozen om aan te ge- 
ven dat het scheidend vermogen van dit 
systeem (evenals dit het geval is met alle 
elektronenoptische systemen) veel meer 
wordt bepaald door wat er met de elek- 
tronen gebeurt als ze de Si-plak met een 
eindige laklaagdikte hebben bereikt. In 

de lak, tegen en in het substraat worden 
de elektronen namelijk verstrooid en 
gedeeltelijk teruggekaatst. En dit kan, af- 
hankelijk van lakdikte en substraatmate- 
riaal, aanleiding geven tot een significante 
lijnverbreding. Bovendien, anders dan in 
elektronenoptische patroongeneratoren, 
worden elektronen die na terugkaatsing 
uit de laklaag te voorschijn komen, door 
het daar heersende elektromagnetische 
veld weer naar de plak teruggebogen en ge- 
ven daar aanleiding tot (ongewenste) be- 
lichting. Dit proces kan zich over verschil- 
lende mm’s uitstrekken en veroorzaakt 
verdere randvertekening. Hier dient nog 
verder onderzoek te worden verricht. 
Niettemin lijkt het, mede gezien onze er- 
varing, verantwoord om te stellen dat dit 
systeem in staat is om in het submicron- 
gebied te werken. 

Hierboven is al aangegeven dat elektronen- 
verstrooiing in lak en substraat aanleiding 
kan geven tot beeldvervorming. Er zijn 
ook andere faktoren die beeldvervorming 
veroorzaken. Beeldvervorming is een ge- 
compliceerd onderwerp en moet worden 
onderverdeeld in reproduceerbare en niet- 
reproduceerbare beeldvervorming. Het zal 
duidelijk zijn dat de niet-reproduceerbare 
beeldvervorming, dat wil zeggen die 
beeldvervorming die van belichting tot be- 
lichting ofwel van laag tot laag verschil- 
lend is, zeer ernstige consequenties kan 
hebben bij de vervaardiging van geïnte- 
greerde schakelingen, waarbij 6 tot 10 
lagen met een grote nauwkeurigheid op 
elkaar moeten passen. 


Zoals we hierboven al gezien hebben, kan 
het magnetisch veld voldoende uniform 
worden gemaakt om de vervorming bene- 
den de 0,1 um te houden, maar met het 
elektrische veld ligt dat heel anders (de 
plak is de anode!). Het elektrische veld 
wordt bepaald door masker, plak en de 
houders hiervan: er zijn onvermijdelijke 
randeffecten en discontinuiteiten die het 
veld verstoren. In de praktijk vinden we 
meerdere microns vertekening aan de 
rand van de plak, maar deze vertekening 
is reproduceerbaar en gaat niet ten koste 
van de geëiste passingsnauwkeurigheid. 
Een ernstige oorzaak van beeldvervorming 


is verandering van plakkromming. Gedu- 


rende het vervaardigen van geintegreerde 
schakelingen verandert de onvlakheid van 
de Si-plak ten gevolge van oxidaties, dif- 
fusies en gedeponeerde lagen. Uit bereke- 
ningen en experimenten is gebleken dat 
op Z2-inch-plakken ongeveer een twintigste 
van de verandering van de plakkromming 
als een dwarsverplaatsing van het beeld 
zichtbaar wordt. Aangezien deze verande- 
ring enkele tientallen microns kan bedra- 
gen is het duidelijk dat dit tot ontoelaat- 
bare passingsfouten leidt. De voor de 
hand liggende oplossing is om de plak te 
dwingen vlak te liggen. Dit kan niet ge- 
beuren, zoals bij optische lithografie, 
door een vacuümaanzuigplaat, maar aan- 
getoond is dat we gebruik kunnen maken 
van de hoogspanning die al in het systeem 
aanwezig is om de plak vlak te houden 
door elektro-statische aantrekkingskracht. 


Door de plak vlak te houden vermijden 
we wel deze vorm van beeldvervorming, 
maar heffen we nog niet de spanning in 
de plak op. Oxidaties, gedeponeerde lagen 
en diffusies veroorzaken namelijk tevens 
spanningen die, in eerste benadering, line- 
aire uitzettingen of inkrimpingen van de 
plak veroorzaken die met het krapper 
worden van de uitrichttoleranties te groot 
zijn om te verwaarlozen. Ook tempera- 
tuurverschillen tussen de opeenvolgende 
belichtings- en uitrichtstappen veroorza- 
ken dit effect, dat overigens alle 1:1-be- 
lichtingssystemen hindert. De elektronen- 
optische beeldprojector kan dit effect 
echter corrigeren. Er is aangetoond dat 
een kleine aanpassing van het magneet- 
veld een geringe vergroting of verkleining 
van de 1:1-maatstaf kan bewerkstelligen. 
Om dit automatisch te regelen, gebruiken 
we de informatie die reeds beschikbaar is 
in de röntgensignalen die vanaf de uit- 
richtkenmerken komen. Hiermee zijn we 
op het aspect van het automatisch uit- 
richten gekomen. 


Het automatisch uitrichten 

Het werkgebied van de elektronenop- 
tische beeldprojector is de produktie van 
patronen met minimale afmetingen in de 
orde van 1 um. Dit betekent dat de ge- 
wenste uitrichtnauwkeurigheid in de orde 
van 0,25 um ligt. In een machine die is 


bestemd voor massaproduktie moet men 
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een uitrichtsysteem hebben dat behalve 
nauwkeurig, zowel snel als automatisch 
werkt. Dientengevolge is er een methode 
ontwikkeld, die gebruikmaakt van de 
remstraling die ontstaat wanneer elek- 
tronen een zwaar element treffen. In deze 
methode wordt een raster van marke- 
ringen (van tantaaloxyde) op de Si-plak 
aangebracht en een identiek patroon in het 
masker. Wanneer het beeld van het raster 
op het masker precies samenvalt met het 
raster op de plak, vallen al de elektronen 
van het masker-raster op het tantaaloxy- 
de; wanneer er echter een fout in het uit- 
richten zit, zullen sommige elektronen op 
het silicium terechtkomen, dat een klei- 
ner rendement heeft voor de opwekking 
van röntgenstralen. Zo zal dus bij perfect 
uitrichten de opbrengst aan röntgenstra- 
len maximaal zijn. Achter de plak gemon- 
teerde detectors zetten de röntgenstraling 
die door de plak heengekomen is om in 
elektrische signalen. De instelling op maxi- 
male intensiteit — d í. het uitrichten — 
geschiedt automatisch. De door de detec- 
tors geproduceerde signalen worden fase- 
gevoelig gedetecteerd, wat tot eenvoudige 
elektronische schakelingen leidt. Er wordt 
zowel aan de linker- als aan de rechterzij- 
de van de plak/het masker een raster van 
markeringen gebruikt, hierdoor kunnen 
niet alleen translatie-, maar ook rotatie- 
fouten worden gedetecteerd en gecorri- 
geerd. De translatiecorrecties worden uit- 
gevoerd door middel van de afbuigspoe- 
len, maar de rotatiecorrecties kunnen niet 
op dezelfde wijze worden uitgevoerd; 

dat zou een onaanvaardbare vertekening 
van het beeld opleveren. In plaats daarvan 
wordt de maskerhouder gedraaid. Het uit- 
richten kan, door gebruik van snelle half- 
geleiderdetectors, binnen de eerste 20% 
van de totale bestralingstijd plaatsvinden. 
De rest van de plak wordt gedurende het 
uitrichten eveneens bestraald, maar dit 
komt, omdat de gebruikte elektronenlak 
PMMA een vrij hoog contrast heeft, niet 
tot uitdrukking in het ontwikkelde beeld. 
Daarmee hebben we het grote voordeel 
dat een aparte belichting voor de uitricht- 
kenmerken niet nodig is. 


Het doel van deze uitrichtmethode is het 
behalen van een uitrichtnauwkeurigheid 
van + 0,1 um ter plekke van de marke- 
ringen. Deze nauwkeurigheid wordt be- 
reikt als we uitgaan van ongeproceste 
plakken en de temperatuur goed in de 
hand hebben. In de praktijk echter, tij- 
dens de vervaardiging van geïntegreerde 
schakelingen, is deze nauwkeurigheid, 
mede ten gevolge van de in de vorige pa- 
ragraaf vermelde plakvervorming, wel iets 
minder, maar is toch vrijwel altijd beter 
dan + 0,25 um. Dit betreft uiteraard de 
uitrichtnauwkeurigheid ter plekke van de 
markeringen: elders op de plak kan de 
passingsnauwkeurigheid door een aantal 
eerder vermelde oorzaken (verandering 
van plakkromming, plakvervorming) min- 


der zijn. 
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Resultaten en conclusies 

Er is nu enige jaren ervaring met een aan- 
tal van deze elektronenoptische beeldpro- 
jectoren, er zijn zowel bipolaire als 
MOST-circuits mee gemaakt en samenvat- 
tend kan worden gesteld dat de machines 
aan hun doelstelling, het automatisch be- 
lichten en uitrichten van structuren tot in 
het submicrongebied, beantwoorden. Als 
sterke punten gelden dus het hoge schei- 
dend vermogen en de hoge automatische 
uitrichtnauwkeurigheid. Bovendien is het 
systeem relatief goedkoop. Als nadelen 
moeten worden genoemd de gevoeligheid 
voor plakvervorming, de effecten ten ge- 
volge van elektronenverstrooiing en de 
noodzaak van het invoeren van tantaaluit- 
richtkenmerken. 


De huidige machines hebben een cyclus- 
tijd van ongeveer 3 minuten, maar het 
moet niet moeilijk zijn om dit in een pro- 
duktiemachine aanmerkelijk te verkorten, 
aangezien van deze 3 minuten slechts on- 
geveer 40 seconden wordt belicht. Het is 
nog te vroeg om uitspraken over proces- 
opbrengst te doen. 


Wel is duidelijk dat de spectaculaire voor- 
uitgang in de laatste jaren van optische 
belichtings- en uitrichtapparatuur de be- 
hoefte aan deze elektronenoptische beeld- 
projector naar de toekomst verschoven 
heeft. Er lijken echter geen fundamentele 
redenen te zijn die dit systeem beletten 
om in het submicrongebied een levensvat- 
bare en economisch haalbare methode 
voor massaproduktie te zijn. 


Enige slotopmerkingen 

Niet alle methoden die in het schema van 
figuur 1 voorkomen zijn even uitgebreid 
besproken. 

Er is bijvoorbeeld nog niets gezegd over 
röntgenstraallithografie, een methode die 
voor het submicrongebied interessant kan 
worden omdat er noch buiging, noch ver- 


strooiing optreedt. 

Daarnaast zijn er nog methoden in ont- 
wikkeling, die op het schema niet voor- 
komen, zoals de elektronenoptische repe- 
teerprojector en, meer recent, ionen- en 
laserstraallithografie. 


Er is in dit artikel niet gepoogd om de 
verschillende lithografische methoden te 
behandelen met een uitgebreidheid die 
evenredig is aan hun belang, voor zover 
dat althans te overzien zou zijn. We heb- 
ben een overzicht gegeven van methoden 
die ons in het VLSI-tijdperk ter beschik- 
king staan. Daarbij hebben we speciale 
aandacht besteed aan licht-optische repe- 
teerprojectie, als voorbeeld van een me- 
thode die zeer zeker een belangrijk deel 
van de lithografie voor zijn rekening zal 
nemen en aan elektronenoptische beeld- 
projectie, als voorbeeld van een methode 
die, hoewel zeer interessant, veel specula- 
tiever is voor wat betreft zijn feitelijke 
toepassing. 


Anderzijds is, gezien tegen het belang van 
de techniek, te weinig aandacht besteed 
aan de methode die elektronenstraalmas- 
kers en licht-optische 1: 1-beeldprojectie 
combineert. Met deze methode zal in de 
komende jaren de grote meerderheid van 
de geavanceerde circuits worden vervaar- 
digd. Dit te meer omdat er ontwikkeling- 
en gaande zijn om deze projectie in het 
verre ultraviolet te laten plaatsvinden, 
waardoor het mogelijk zal zijn om nog 
kleinere afmetingen af te beelden. 





Fig. 5. SEM-opname van een bipolaire 
hoogfrequent transistor. De verschillende 
(4) lagen van deze transistor zijn gedefi- 
nieerd door met behulp van een elek- 
tronenoptische patroongenerator direct 
op de plak te schrijven. Elk van de witte 
streepjes stelt 10 micron voor. 
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De algemeen aanvaarde visie op het VLSI- 
tijdperk voor wat betreft de lithografie is 
dan ook dat tot in het 1-microngebied 
licht-optische 1:1-beeldprojectie en licht- 
optische repeteerprojectie de feitelijk toe- 
gepaste methoden zullen zijn. 


Wat er daarna, in het submicrongebied, 
zal gebeuren is veel minder duidelijk. Zo- 
wel röntgenstraallithografie als elektro- 
nenoptische beeldprojectie maken daar 
een kans, het vereist echter een apart ar- 
tikel om de verschillende fysische proble- 
men, moeilijkheden en mogelijkheden van 
deze (en andere!) methoden evenwichtig 
tegen elkaar af te wegen. Er is voor dit 


submicrongebied nog minstens één andere 
kandidaat en dat is de methode die op het 
schema van figuur 1 als ‘direct op de plak’ 
staat aangegeven. 


We willen dit artikel besluiten met een re- 
sultaat van deze methode. In dit geval 
wordt het patroon, dat in digitale vorm 
op de magneetband staat, direct door de 
elektronenoptische patroongenerator op 
de plak geschreven. Het is duidelijk dat 
deze methode vele voordelen biedt; het 
grootste nadeel is dat bij de huidige 
schrijftijden deze manier van patroonde- 
finitie veel te duur is om voor gebruik in 
massaproduktie in aanmerking te komen. 


Figuur 5 toont een SEM-opname van een 
bipolaire hoogfrequent-transistor met 
submicronbasis en emittercontacten; 
naast de scherpte van het patroon valt 
vooral de perfecte (automatische) uit- 
richtnauwkeurigheid op. Ook voor deze 
methode zal er in de toekomst zeker een 
plaats zijn. 
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Schottky-diode versus pn-diode 


M. Kleefstra 


De stroom-spanningskarakteristieken van Schottky-dioden en die van pn-dioden in 
matig en homogeen gedoteerde halfgeleiders worden gegeven door vrijwel dezelfde 
formules. Bij de gebruikelijke pn-dioden wordt de effectieve snelheid van de ladings- 
dragers begrensd door diffusie van ladingsdragers. In Schottky-dioden wordt de snel- 
heid van de meerderheidsladingsdragers bepaald door de thermionische emissie van die 
ladingsdragers en moet de verlaging van de potentiaalberg door de beeldkracht worden 
meegenomen. Of het ene type diode te verkiezen is boven het andere type hangt af van 


de toepassing en van de materiaalkeuze. 


In geïntegreerde schakelingen komen 
alom metaal-halfgeleider- en p-n-overgan- 
gen voor. Bij matig gedoteerde halfgelei- 
ders werken deze overgangen als gelijk- 
richter en spreekt men van respectievelijk 
Schottky- en pn-diode. Beide typen 
dioden worden vaak apart en onvolledig 
behandeld!,2,3, zodat de overeen- 
komsten en verschillen moeilijk te zien 
zijn. De Schottky-diode alleen wordt dui- 
delijk door Rhoderick* behandeld, de 
pn-diode door Grove”. 

Hier zal worden aangetoond dat voor 
dioden met matige dotering, van 10° — 
102* m”° in silicium, de klassieke for- 
mules voor het /-V-verband dusdanig kun- 
nen worden geschreven dat overeen- 
komsten en verschillen zo te zien zijn. 
Voor dioden met hogere doteringen zijn 
de formules zeer verschillend en deze 
zullen hier buiten beschouwing blijven. 
Dioden met matige dotering worden 
gebruikt in gelijkrichters, zonnecellen, 
basis-collectorovergangen in transistoren 
en in CCD's. 

De beschouwingen zullen worden beperkt 
tot het gedrag rond kamertemperatuur, 
de dotering van de halfgeleiders wordt 
homogeen genomen, de rand van de ont- 
ruimingslaag wordt oneindig scherp geno- 
men en de weerstand van contacten op de 
niet-diode-kant van de halfgeleider wordt 
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steeds verwaarloosd. De gebruikte symbo- 
len worden in tabel I toegelicht en voor 
de halfgeleiders silicium en gallium- 
arsenide worden waarden voor de para- 
meters aangegeven. 

Het /-V-verband zal met behulp van drie 
mechanismen worden berekend: 

— de recombinatie en generatie van 
ladingsdragers in de ontruimingslaag, 

— de thermionische emissie van ladings- 
dragers over een potentiaalberg, 


Fig. 1. Energieschema van de Schottky- 
diode. 
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Prof. dr. M. Kleefstra (49), sinds 1 ja- 
nuari 1980 hoogleraar, studeerde schei- 
kunde te Groningen en promoveerde 
aldaar in 1961 op een proefschrift over 
kleurcentra. Van 1960 tot 1969 werkte 
hij aan elektronische componenten bij 
Thomson-CSF in Frankrijk. Sindsdien is 
hij verbonden aan de Afdeling der Elek- 
trotechniek van de Technische Hoge- 
school te Delft, vanaf 1977 als lector in 
de elektrotechnische materialen. Hij 
houdt zich daar bezig met Schottky- 
dioden en met CCD's die door pn-dioden 
worden aangestuurd. 





— de diffusie van ladingsdragers als min- 
derheden in een stuk halfgeleider nadat 
de ladingsdragers de potentiaalberg heb- 
ben ‘genomen’ dankzij de thermionische 
emissie. Met behulp van de Fermi-Dirac- 
statistiek en de diffusievergelijking J, = 
qD dn/dx zijn de formules voor het /-V- 
verband af te leiden. Pas na numeriek 
uitrekenen kan worden beslist welk 
mechanisme overheerst bij de gekozen 
parameters. 


De I-V-karakteristiek van de Schottky- 
diode 

De toegelaten energieniveaus van een 
Schottky-diode op een n-halfgeleider zijn 
in figuur 1 aangegeven. De bandbuiging 
tengevolge van de spiegelbeeldkracht is 
voor E‚ en E, aangegeven. In figuur Ì 
wordt het potentiaalverschil Vr — V 
opgebouwd door de positieve ruimte- 
lading in de ontruimingslaag met dikte 
W4. Door de werking van de spiegelbeeld- 
kracht wordt de potentiaalberg verlaagd 
met 








symbool definitie eenheid waarde(n) 

GaAs 
D, diffusiecoëfficiënt van elektronen m?°/s 3 x 10 2x 107? 
Dy diffusiecoëfficiënt van gaten " EE KE ie 8x 10% 
Ee onderkant van de geleidingsband eV 
Er Fermi-niveau „ 
E‚ bovenkant van de valentieband 7 
E‚ breedte van de verboden band " ne 1,43 
J stroomdichtheid Am”? 
bin difussielengte van de elektronen m Ln =VPnr 
Lp difussielengte van de gaten Z Lo=VDpT 
jà absolute temperatuur K 300 300 
Na concentratie van acceptoratomen m” 10 10 (a) LOE — 102 
Np concentratie van donoratomen oe Z Z 
Vv aangelegde voorwaartse spanning V 
Vp diffusiespanning n° 0-1 n 0-1 
Wa dikte van de depletielaag m 0O-10f 0O—10® 
W. dikte van de geleidende n-laag " 0O-10® 0 —10 
Wo dikte van de geleidende plaag z 0-10 O—10” 
k constante van Boltzmann Joule/K 1,38 x 1072 1,38 x 1072 
n concentratie van elektronen üm 010% O0 
ni = n Xp bij evenwicht m$ 2 x 10% 1 x 10% 
p concentratie van gaten m$ 0-10 (a) 0 — 10% 
Vin thermische snelheid ladingsdragers mes” 1 x 105 1 x 105 
e permittiviteit F-m 1,0 x 1071 1,0 x 107 
7 levensduur van ladingsdragers s 1010 _ 10 z10® (b) 


(a) 





TABEL I. Gebruikte symbolen. 


(a) min of meer instelbaar. 


(b) zie referentie 6. 





Ap= == (oe) 7 [Wp VND] *. 


ni/2 \2e 


Wordt de stroom in deze diode bepaald 
door generatie of recombinatie van 
ladingsdragers in de ontruimingslaag, dan 
is een benaderde formule voor de stroom- 
dichtheid | 


HM: 
Jg =4: Wa En (EIVI2KT_ 1). (1) 
T 


Hierbij wordt de dikte van de ontrui- 
mingslaag gegeven door 


Wa = (Ze d, Is | 


qNp 
Kan de generatie-recombinatie in de ont- 
ruimingslaag verwaarloosd worden, dan 
wordt de stroom van elektronen bepaald 
door hun thermionische emissie over de 
potentiaalberg en wordt verkregen 


Vith 


‚Np: e-WVp-ApP)/KT 
V 6 





Jn=4: 


(UVIT _ 1). (2) 


De gatenstroom wordt bepaald door dif- 
_fusie van gaten in het n-gebied en hierin 
treedt de grootheid W‚ op indien Lp >> 
W‚, maar Lp treedt op wanneer Lp $<< 
W. Voor de diode met eindige lengte, 
dat wil zeggen Lp >> W‚, is de gaten- 
stroomdichtheid 

Dp: ni 

(EU /T_Ĳ. 


3a 
W-Dr (3a) 


J,=q" 
10 


Voor het geval Lo << W‚ geldt: 
Do: ni 


: SUI 1) 
Ken ( ) 


Jp = q (3b) 
De totale stroomdichtheid is òf J‚g òf Jo 
+J,. Van die twee moet de grootste 
waarde bij een bepaalde spanning worden 
genomen. Alleen bij het punt waar J‚g = 
J„ +Jp zijn deze formules niet te gebrui- 
ken. 


De I-V-karakteristiek van de pn-diode 
Het bandenschema van de pn-diode is in 
figuur 2 aangegeven. Omdat hier geen 
overgang perfecte geleider-ontruimde 
halfgeleider bestaat, treedt er geen beeld- 
krachtverlaging op. De stroomdichtheid 
te danken aan generatie-recombinatie in 
de ontruimingslaag is hier weer 


(4) 


n; 
Jard We (eIVLKT _ 1). 
27 


Hierbij is ook weer 
2e( Vp haag V) 
De Ee 
qNp 


| 1/2 


Bij verwaarlozing van de stroom volgens 
(4) wordt de elektronenstroom bepaald 
door de diffusie van elektronen in het 
p-gebied volgens de formules: 





D 
Jg Np eIVDIRT (UVIT 1) 
Wo 
voor Ln >> W‚ (5a) 


n 


‚Np: e IVDIKT (EUVT _j) 
En 





J„=4: 


voor Ln << Wp. (5b) 


Het is gebruikelijk deze formules te schrij- 
ven als bijvoorbeeld 
Dn ’ ni 


Ter dn 
p “A 


Jn=4" 


maar hierbij is vergelijken met formule 
(2) moeilijker. De gatenstroomdichtheid 
is weer voor de diode met eindige lengte 


2 
Dy: ni 


e4V/KT | 
We No ( ) 


Jp=q" (6) 


Wordt een zeer grote stroom door een 


Fig. 2. Energieschema van een pn-diode. 
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pn-diode geleid waarin voornamelijk de 
gaten de stroom verzorgen dan moet, 
vanwege de elektroneutraliteit op het 
punt P de elektronenconcentratie toene- 
men en zal aldaar gelden p =n =n; 
e4V/2KT De gatenstroomdichtheid wordt 
dan 
2 
5 Hi 
W‚ . Hi; 





Je =q. ga VIa 


en verandert dus minder snel met de aan- 
gelegde spanning dan in formule (6). De 
totale stroom door de diode is weer J‚, of 
Jn +Jp waarbij de grootste waarde moet 
worden genomen. 


Overeenkomsten en verschillen 

Bij pn-dioden is de diffusiespanning een- 
voudig te berekenen uit de doteringen in 
de p- en de n-gebieden. Omdat er slechts 
één fase aanwezig is geldt: diffusiespan- 
ning = verschil in uittree-arbeiden = con- 
tactpotentiaal. Bij een metaal-halfgelei- 
derovergang is er sprake van een contact 
tussen twee fasen. De grootte van de dif- 
fusiespanning hangt zowel af van het ver- 
schil in uittree-arbeiden als van het aantal 
oppervlaktetoestanden (Tamm-toestan- 
den) op de halfgeleider en de ligging van 
de oppervlaktetoestanden ten opzichte 
van de verboden band. De variatie die in 
de diffusiespanning mogelijk is, is bij 
beide typen dioden beperkt en globaal 
kan gesteld worden dat een Schottky- 
diode wat dit betreft nogal overeenstemt 
met een p*n-diode waarbij door variëren 
van hetzij het metaal, hetzij de dotering 
in het p'-gebied de diffusiespanning kan 
worden veranderd van ongeveer 0,4 tot 
ongeveer 1,0 volt. 

In de /V-karakteristieken zijn de meeste 
stroomcomponenten voor beide typen 
dioden dezelfde, nl. de stroom Jg en de 
gatenstroom J. Bij de elektronenstroom 
geeft de Schottky-diode een hogere 
stroomdichtheid om twee redenen: 

— de ingebouwde potentiaalsprong Vp 
wordt door de beeldkracht met Aw ver- 
laagd wat vooral bij sperren invloed heeft, 
— bij de Schottky-diode hebben de elek- 
tronen in figuur Ì naar links de hoogst 
mogelijke snelheid, nl. vn /N/6r terwijl in 
het p-gebied van de pn-diode de elektro- 
nen door diffusie hun gemiddelde effec- 
tieve snelheid tot D/W gereduceerd 
zien. 

Bij gelijke Vj) en doteringen heeft de 
Schottky-diode een veel hogere verzadi- 
gingsstroom dan de gebruikelijke pn- 
diode maar een pn-diode met zeer kleine 
W‚, zal een bijna evengrote verzadigings- 
stroom bezitten. Vanwege de hogere ver- 
zadigingsstroom kan bij Schottky-diodes 
de stroom J‚, vaak worden verwaarloosd 
en bovendien treedt ook geen term exp 
(qV/2KT) in plaats van exp (q V/KT) op 
bij hogere stroomdichtheden. 

De opslag van minderheidsladingsdragers 
in pn-dioden verhindert het zeer snelle 
schakelen van deze dioden. Om snelle 
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schakeltijden in pn-dioden te bereiken 
wordt de levensduur 7 verkort door 
inbouwen van recombinatie-generatie- 
centra. 

Voor toepassingen in gelijkrichters is de 
Schottky-diode te verkiezen boven de 
gewone pn-diode, bij dezelfde stroom kan 
met een kleiner oppervlak worden vol- 
staan, is de schakeltijd korter en geldt 
over een groot gebied dat 


J=J, (eFV/T _1). 


Bij gebruik in zonnecellen is het gewenst 
dat de diodestroom volgens de hier ge- 
geven formules zo klein mogelijk is. 
Omdat de zonnecellen veelal zodanig wor- 
den gebruikt dat de diode tamelijk ver 
voorwaarts staat is J,, te verwaarlozen. 
Het trage diffusiemechanisme is bij uit- 
stek geschikt om de diodestroom laag te 
houden. In zonnecellen zijn dan ook pn- 
dioden te verkiezen boven Schottky- 
dioden. In CCD’s wordt gebruikgemaakt 
van gesperde dioden waarvan de lek- 
stroom zo laag mogelijk moet zijn om de 
CCD's bij lage frequenties te kunnen 
bedrijven. In het algemeen zijn hier dus 
pn-diodes te verkiezen boven Schottky- 
diodes. 


Praktische dioden 

Bij het maken van dioden treden de vol- 
gende problemen op: 

— het doteren van de halfgeleider door 
donor- of acceptoratomen op rooster- 
plaatsen in te bouwen, 

— het maken van een zuiver metaal-half- 
geleidercontact, | 
— het lateraal begrenzen van de diode, 
— het aanbrengen van een goed gelei- 
dend contact op de niet-diodekant van de 
halfgeleider. 

Het doteren van halfgeleiders en in het 
bijzonder van silicium behoort tot de 
relatief gemakkelijkst oplosbare proble- 
men. 

Voor het maken van een zuiver metaal- 
halfgeleidercontact zijn de mogelijkheden 
beperkt. Bij het opbrengen van een 
metaal door sputteren of opdampen 
onder vacuüm is het moeilijk een schoon 
en oxydevrij halfgeleideroppervlak te 
verkrijgen. Schoonsputteren van het 
oppervlak leidt vaak tot beschadiging van 
het kristalrooster. Het neerslaan van 
metalen uit vluchtige verbindingen voor- 
afgegaan door etsen in de gasfase lijkt de 
meest aantrekkelijke methode maar 
wordt weinig gebruikt. Meestal worden 
Schottky-dioden gemaakt zodanig dat het 
originele onzuivere halfgeleideroppervlak 
binnen het metaal komt te liggen. Twee 
soorten metalen komen hierbij in aanmer- 
king. Metalen die SiO, kunnen reduceren 
zoals Al en Ti kunnen worden gebruikt 
om de SiO, -laag op Si te doen verdwij- 
nen. Metalen die silicides vormen zoals Pt 
kunnen worden gebruikt om een gedeelte 
van het silicium in silicide om te zetten. 
Bij beide groepen metalen zijn temperatu- 
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Fig. 3. Verschillende dioden met 
kromtestraaleffect. a. pn-diode met r > 0. 
b. Schottky-diode met r = 0. c. Schottky- 
diode omringd met pn-diode. d. 
Schottky-diode omringd met MOS-capaci- 
teit. | 





ren tussen 500 en 1000 °C nodig zodat 
diffusie van het metaal in het silicium, 
afscheidingen bij rekristallisatie en derge- 
lijke niet altijd te voorkomen zijn. Ook 
omdat het grensvlak hier niet een ideaal 
vlak wordt is het met deze methoden 
moeilijk perfecte Schottky-dioden te 
maken. | 
Het lateraal begrenzen van de dioden is 
nodig vanwege het kromtestraaleffect wat 
maakt dat het elektrische veld aan de 
rand van de diode hoger is dan midden _ 
onder de diode. Pn-dioden kunnen ge- 
makkelijk met een grote kromtestraal 
worden gemaakt en dotering en geometrie 
van losse dioden kunnen zo worden geko- 
zen dat het veld aan de rand voldoende 
klein blijft om hier geen doorslag te ver- 
oorzaken. Bij Schottky-dioden kan de 
laterale begrenzing worden verkregen 
door omringing met een MOS-capaciteit 
of met een pn-diode. Bij zeer kleine 
Schottky-dioden, zoals in geïntegreerde 
schakelingen, is omringing met een MOS- 
capaciteit het gemakkelijkst. Voorbeelden 
van dioden zijn in figuur 3 schematisch 
aangegeven. 


Slotopmerkingen 

Om met succes dioden te kunnen gebrui- 
ken is kennis van de eenvoudige theorie 
aanbevelingswaardig. Naast theoretische 
overwegingen moet ook rekening worden 
gehouden met het gemak van fabricage. 
Dit laatste heeft gemaakt dat er momen- 
teel meer pn-dioden worden vervaardigd 


dan Schottky-dioden, hoewel de laatste 
| 11. 


honderd jaar ouder zijn en in theorie 
soms zijn te verkiezen, vooral als het om 
gelijkrichters gaat. 
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SiO-lagen op Si 


M.W. Hillen en J.F. Verwey 


De kwaliteit van geïntegreerde schakelingen wordt voor een belangrijk deel bepaald 
door de eigenschappen van de SiO,-laag en van het Si/SiO-grensvlak. Daarom vormen 
SiO,-lagen op Si een belangrijk onderzoekgebied van de fysica van halfgeleiderelemen- 


ten. 


Het belang van SiO, -lagen en van de kwaliteit van deze lagen in moderne geïntegreerde 
schakelingen wordt aangegeven. De MOS (metaal-oxyde-silicium)-condensator blijkt bij 
uitstek een meetinstrument voor de eigenschappen van de SiO,-laag. Beweeglijke ionen 
veroorzaken tal van ongewenste effecten in SiO, lagen. Op dit gebied is interessant 
onderzoek te doen aan de wijze waarop de ionen zich in de SiO, -laag (al of niet in een 
chlooromgeving gegroeid) voortbewegen en hoe zij worden ingevangen aan het 


Si/SiO, -grensvlak. 


Het belang van SiO-lagen in 
geïntegreerde schakelingen 

De grote vlucht die de elektronica heeft 
genomen, is te danken aan de ontwikke- 
ling van de monolitische geïntegreerde 
schakeling. Hierbij zijn de elementen van 
een circuit in één substraat aangebracht. 
In de meeste gevallen is dit substraat een 
schijf monokristallijn silicium en vindt de 
fabricage plaats met de planaire techno- 
logie, waarin siliciumdioxyde (SiO, ) een 
belangrijke rol speelt. 

In de planaire technologie wordt de sili- 
ciumschijf (de ‘plak’) voorzien van een 
laag SiO, (dikte ca. 100 nm), meestal 
door oxydatie bij hoge temperatuur 


(1000-1200 °C) in een zuurstofomgeving. 


In de oxydelaag worden vervolgens met 
een foto-lithografische techniek vensters 
geëtst, waardoor doopstoffen in het Si- 
kristal kunnen worden gebracht. Op deze 
wijze ontstaan dan n- of p-type gedo- 
teerde gebieden. Dit proces wordt even- 
tueel herhaald en tenslotte wordt een 
metaalpatroon over het oppervlak aange- 
bracht dat de verbindingen vormt tussen 
de verschillende circuitelementen. 


Voor de hier beschreven maskerende 
werking van het oxyde worden niet zo 
hoge eisen aan de kwaliteit van het Si0 
gesteld. Om andere redenen gebeurt dat 
wel. 


Allereerst 1s het voor een goede werking 
van transistors van belang dat de ruimte- 
ladingsgebieden van de pn-juncties niet 
12 


aan het reactieve siliciumoppervlak 
komen, maar door een laag SiO, worden 
beschermd (stabiliserende werking). In de 
tweede plaats is het belangrijk dat er 
weinig oppervlaktetoestanden zijn (passi- 
verende werking). Oppervlaktetoestanden 
zijn plaatsen aan het Si/SiO, -grensvlak, 
waar vrije ladingsdragers uit het Si inge- 
vangen kunnen worden. Dat wil zeggen, 
zij introduceren energieniveaus in de ver- 
boden band van het silicium en worden 
daarom gemeten als functie van de ener- 
gie (cm”?eV”*). Hun aantal is zeer afhan- 
kelijk van de wijze van oxyderen en van 
de nabehandeling. 


Vooral de MOSFET (Metal-Oxide-Semi- 





Fig. 1. Schematische doorsnede van een 
p-kanaal-MOSFET. Enkele typische afme- 
tingen zijn aangegeven. Let wel, de figuur 
is niet op schaal getekend. De metaallaag 
is zwart weergegeven. 
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conductor Field-Effect Transistor) is hier 
zeer gevoelig voor. In figuur 1 is een 
p-kanaal-MOSFET getekend; het kanaal 
ontstaat doordat een voldoend grote 
negatieve spanning op de veldelektrode 
(‘gate’) gaten naar het oppervlak trekt. Er 
is nu een geleidende verbinding tussen het 
source- en het draingebied. Oppervlakte- 
toestanden kunnen het stroomtransport 
in het kanaal van een MOSFET ernstig 
belemmeren. Dit probleem was in het 
verleden zo groot, dat, hoewel in principe 
de werking van de veldeffect-transistor in 
1930 bekend was, het toch tot de jaren 
zestig heeft geduurd voordat deze kon 
worden gefabriceerd. 


Hoe het aantal oppervlaktetoestanden in 
de loop der jaren drastisch is verminderd, 
is te zien in figuur 2. Deze afname is voor- 
namelijk bereikt door schonere fabricage- 
methoden en betere materialen. Voor een 
MOSFET in een modern geïntegreerd 
circuit is het aantal bewegende ladings- 
dragers in het kanaal ongeveer 10! per 
cm?. Het aantal oppervlaktetoestanden 
mag maar een fractie hiervan zijn om een 
ongestoord transport mogelijk te maken. 
Het minimaal te bereiken niveau is op dit 
moment ongeveer 108 cm? Dit betekent 
dat slechts 1 op de 10” oppervlakte-ato- 
men werkzaam is als oppervlaktetoestand. 
Voor de MOS-veldeffecttransistor is het 
bovendien van belang dat het oxyde 
onder de veldelektrode weinig defecten 
(minder dan 10 per cm?) en weinig lading 
bevat. Defecten (zwakke plekken in het 





Fig. 2. Afname van de oppervlaktetoe- 
standsdichtheid Ns, in het midden van de 
verboden band als functie van het jaartal. 
De gegevens zijn afkomstig van Goetzber- 
ger, e.a.* Jen uit recente publicaties. 
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Fr. M.W. Hillen (31) studeerde technische 
natuurkunde aan de Rijksuniversiteit te 
Groningen. In 1976 behaalde hij het inge- 
nieursdiploma na een onderzoek naar 
recombinatiestromen in epitaxiale diodes. 
Bij de onderzoekgroep Fysica van Half- 
geleiderelementen verricht hij sindsdien 
onderzoek naar de invloed van beweeg- 
lijke lading in SiO, op de basisstroom van 
bipolaire transistoren, in het kader van 
een promotieopdracht onder leiding van 
prof. dr. J.F. Verwey. 


oxyde) veroorzaken doorslag bij kleine 
elektrodespanningen en aangezien één 
defecte transistor in principe een gehele 
geïntegreerde schakeling (met honderden 
elementen) onbruikbaar maakt, is het 
onderzoek naar defectvrije oxydes in 
industriële sfeer van groot belang voor de 
opbrengst. 


Lading in het oxyde kan beweeglijk zijn 
en vast. Vaste oxydelading bestaat uit 
positieve ladingen (bijvoorbeeld extra Si), 
die zich in de eerste 20 nm van het oxyde 
(gerekend vanaf het Si/SiO,-grensvlak) 
bevinden. Zij kunnen geen lading met het 
silicium uitwisselen, zijn onbeweeglijk in 
een sterk elektrisch veld en onvermijdelijk 
aanwezig. Door een geschikte nabehande- 
ling is hun aantal gering genoeg te maken. 
Tengevolge van deze ladingen is bijvoor- 
beeld in een n-kanaal-MOSFET het kanaal 
al aanwezig zonder dat een elektrodespan- 


Fig. 3. Doorsnede van een p-kanaal- 
MOS-condensator. De elektrodespanning 
is altijd gedefinieerd ten opzichte van het 
substraat. Een negatieve elektrodespan- 
ning veroorzaakt een positieve ruimte- 
lading, omdat de vrije elektronen van het 
grensvlak worden weggeduwd en er dan 
geïoniseerde vaste donoratomen in het 
kristal achterblijven. 
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ning is aangelegd. Beweeglijke lading in 
het oxyde blijkt ook positief te zijn en 
bestaat uit ionen van verschillende veront- 
reinigingen die tijdens het technologisch 
proces worden ingebouwd. Deze ionen 
bewegen onder invloed van een elektrisch 
veld en veroorzaken zo instabiliteiten in 
de werking van een transistor. Zij zijn 
slechts te vermijden door de oxydatie 
onder zeer schone omstandigheden uit te 
voeren. Wij komen op de beweeglijke 
lading nog uitvoerig terug. 


Onderzoek met behulp van de MOS- 
condensator 

Uit het vorige is duidelijk geworden dat 
de MOS-veldeffecttransistor zeer gevoelig 
is voor de kwaliteit van het SiO, en van 
het Si/SiO,-grensvlak. Een eenvoudiger 
structuur waarmee de eigenschappen van 
de SiO, -laag en het Si/SiO -grensvlak 
kunnen worden bepaald, is de MOS-con- 
densator (figuur 3). De spanning (VG ) die 
op de veldelektrode van de MOS-conden- 
sator wordt aangelegd verdeelt zich over 
het oxyde (Vox) en de ruimteladingslaag 
in het silicium onder de elektrode (de 
oppervlaktepotentiaal W‚): Vo = Vox + 
Ws. (De spanningsval over de rest van het 
silicium is verwaarloosbaar.) 


Er is één elektrodespanning waarbij Ws 

= 0 (de energiebanden in het silicium bui- 
gen dan niet) en dat is de vlakkebanden- 
spanning Verg. Voor een ideale MOS-con- 
densator is Vr B = 0, maar tengevolge van 
enerzijds het verschil in uittreepotentiaal 
tussen metaal en silicium en anderzijds 
oppervlaktetoestanden, vaste en beweeg- 
lijke lading in het oxyde, zal Vrg in de 
praktijk ongelijk aan nul zijn. 


De invloed van een lading Q ter plekke x 
in het oxyde (figuur 4) op de vlakkeban- 
denspanning is: Vr =-xQ/dox Co, 
waarin C‚ de capaciteit van de oxyde- 
laag is en dx de dikte. Voor een bewe- 
gende lading Q geldt dan dat deze aan het 
Si/SiO, -grensvlak (x =d x) een maximale 
invloed uitoefent: Vr =— Q/C, en aan 


Fig. 4. Een lading Q ter plekke x 
(oSx<d,) in het oxyde induceert beeld- 
ladingen in de metaallaag (Qp) en in het 
Si (Qs;). Uit Om kan Verg worden bere- 
kend. 
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het metaal/SiO, -grensvlak totaal geen 
invloed uitoefent (Vr g = 0). Duidelijk is 
dan ook dat vaste lading aan het grensvlak 
en lading in oppervlaktetoestanden een. 
bijdrage aan Vrg leveren. 


Het is van belang de hoeveelheid lading te 
kennen. De hoeveelheid beweeglijke 
lading is op de volgende manier te meten. 
De MOS-condensator wordt opgewarmd 
tot een temperatuur 7 (bijvoorbeeld 
200°C). Dan wordt er een positieve span- 
ning op de veldelektrode aangesloten 
Va = V, > 0. De beweeglijke ionen zullen 
nu onder invloed van het elektrisch veld in 
het oxyde naar het Si/SiO, -grensvlak 
gaan. Het blijkt experimenteel dat de 
ionen in het oxyde blijven; zowel het 
Si/SiO,- als het SiO; /metaal-grensvlak 
laten geen ionen door. 
Na een tijd t, wordt de plak weer afge- 
koeld, waarbij de spanning V, gehand- 
haafd blijft om de ionen op hun plaats te 
houden. Bij kamertemperatuur zijn de 
ionen minder beweeglijk en hier bepalen 
we de vlakkebandenspanning Verg. Dan 
wordt dit experiment herhaald met ge- 
bruikmaking van een negatieve spanning 
Va = -V,, waardoor de ionen naar het 
metaal/SiO, -grensvlak worden getrokken. 
Na afkoeling meten we een vlakkebanden- 
spanning Vrg en uit de verschuiving 
AVrgB=|VrB — Vgl volgt het aantal 
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beweeglijke ionen: 


N= — Co: AVrp 


Al —= 


Dit BTS (‘Bias-Temperature Stress’)-ex- 
periment is een standaardmeting in vele 
laboratoria. 


Ongewenste effecten veroorzaakt door 
beweeglijke ionen in SiO, 


Instabiliteit van de drempelspanning in 
een MOSFET 


De drempelspanning is die spanning die 
op de veldelektrode van een MOSFET 
moet worden aangelegd, om juist een 
kanaal te laten ontstaan. Aangezien de 
drempelspanning afhankelijk is van de 
vlakkebandenspanning, kan deze variëren 
als er beweeglijke ionen in het oxyde aan- 
wezig zijn. De eis is nu voor een moderne 
MOSFET: een spanningsverschuiving ten 
gevolge van instabiliteiten van maximaal 
0,05 V. Dit komt neer op minder dan 
2x10!° beweeglijke ionen per cm? 
(minder dan 1 op de 100.000 atomen in 
een monolaag). 


In de technologie zijn verschillende 
methoden in gebruik om tot zo’n lage 
concentratie te komen. We noemen hier 
slechts een zeer bekende en misschien wel 
de meest effectieve: schoon werken. Eén 
aspect hiervan is bijvoorbeeld het schoon- 
maken van de ovenbuis vóór de oxydatie. 
Dit schoonmaken gebeurt door gedurende 


Fig. 5. De invloed van de tijd, dat de 
oxydatie-oven wordt schoongespoeld op 
de concentratie van beweeglijke ionen in 
het oxyde (N;) verkregen door oxydatie 
na het schoonspoelen. De chloorverbin- 
ding, waarmee wordt schoongespoeld is 
C‚H3Cl3, dat bij hoge temperatuur equi- 
valent is met het veel agressievere HCI. 
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reiniging : 1200°C 
in 2% C‚H,Cl,/0, 
oxydatie : 1200°C 
in O,, 30 min. 
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enkele uren HCI door de ovenbuis te lei- 
den. Hierdoor wordt de ovenwand uitge- 
loogd en gedurende enige tijd daarna ver- 
wachten we geen verontreinigingen uit de 
ovenwand. Dat dit effect inderdaad 
optreedt blijkt uit figuur 5, waar de con- 
centratie van beweeglijke ionen in het 
oxyde is weergegeven als functie van de 
tijd dat de oven werd schoongemaakt 
vóór de oxydatie. 


Creatie van extra oppervlaktetoestanden 
door beweeglijke ionen 


Al in 1968 veronderstelde Goetzberger' 


dat een positief ion aan het Si/SiO» -grens- 


vlak een elektron in het Si kan aantrek- 
ken. Dit elektron kan dan niet meer vrij 
door het Si-kristal bewegen en er is een 
oppervlaktetoestand ontstaan. Voor de 
beweeglijke ionen betekent dat, dat er 
meer oppervlaktetoestanden zouden moe- 
ten zijn als de ionen zich aan het Si/SiO,- 
grensvlak bevinden dan als de ionen aan 
het metaal/SiO, -grensvlak zijn. 


Dat een dergelijk effect inderdaad 


optreedt is pas zeer onlangs aangetoond. 


Het probleem hierbij is dat het aantal 
oppervlaktetoestanden moet worden 
bepaald bij een temperatuur, waarbij de 
ionen niet meer beweeglijk zijn (7 < 
—30 °C). Snel? heeft een lineair verband 
gevonden tussen het aantal beweeglijke 
ionen aan het Si/SiO, -grensvlak en het 
aantal oppervlaktetoestanden in de mid- 
delste 0,8 eV van de verboden band van 
Si (met een methode die werkt bij 7 = 
—196 °C). Het is echter niet zo dat elk 


ion juist één oppervlaktetoestand veroor- 
zaakt (figuur 6). 


Veroudering 


Geen enkel halfgeleiderelement heeft het 
eeuwige leven, ook al wordt het geheel 
volgens de voorschriften in een schakeling 


Fig. 6. Verband tussen de oppervlakte- 
toestandsdichtheid en het aantal beweeg- 
lijke ionen N,. N; is bepaald met de TSIC- 
methode (zie verderop ). 
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gebruikt. De levensduur van een MOSFET 
wordt vaak bepaald door de tijd die het 
duurt voor een gate-oxyde elektrisch 
doorslaat. Aangezien de levensduur jaren 
kan en moet bedragen, worden levens- 
duurtests vaak onder verzwaarde omstan- 
digheden (hoge temperatuur en hoge elek- 
trische spanning) uitgevoerd. Het BTS- 
experiment is hier zeer geschikt voor. 
Veroudering is het verschijnsel dat de 
kwaliteit van een SiO, laag afneemt bij 
langdurig normaal gebruik. Het uit zich in 
een fatale elektrische doorslag van het 
oxyde. Veroudering kan worden onder- 
zocht met BTS-experimenten, die worden 
uitgevoerd bij een spanning, die lager is 
dan de doorslagspanning. De doorslag- 
spanning is de spanning waarbij het oxyde 
direct doorslaat, dus vóórdat het BTS- 
experiment wordt toegepast. Figuur 7 
toont experimentele resultaten, waaruit 
blijkt dat als de concentratie van beweeg- 
lijke ionen in het oxyde hoger is, de ver- 
oudering sneller optreedt. 


Een mogelijke verklaring is dat bij posi- 
tieve elektrodespanning het gecombi- 
neerde effect van positieve elektrode- 
lading en positieve ionlading aan het 
Si/SiO»-grensvlak elektronen in het 
oxyde trekt. Deze elektronen, versneld 
door het veld, veroorzaken stoot-ioni- 
satie, waardoor weer extra elektronen 
ontstaan, enzovoort. Dit lawine-effect 
veroorzaakt dan doorslag van de SiOs - 
laag. 


Of dit inderdaad het verouderingsproces 
is, is echter nog niet duidelijk. Dat er 


Fig. 7. Maximale levensduur tax van 
een MOS-condensator tijdens een verou- 
deringsexperiment als functie van de aan- 
gelegde spanning V; tijdens het BTS-ex- 
periment. Parameter is de concentratie 
van beweeglijke ionen N; ae 
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ionen bij betrokken zijn lijkt waarschijn- 
lijk uit de experimenteel gevonden afhan- 
kelijkheid van de veroudering van de toe- 
gepaste elektrodespanning en de ge- 
bruikte temperatuur. Toenemende kans 
op veroudering blijkt op te treden bij het 
toepassen van steeds dunnere oxydenen j 
bij het gebruik van elektronen- en laser- 
‚bundels in de technologie. Dit maakt ver-: 
| der onderzoek van de verouderingspro- 
cessen noodzakelijk. 


_De eigenschappen van beweeglijke ionen 
in SiO,-lagen 

Uiteraard kunnen de ongewenste effecten 
van beweeglijke ionen worden onder- 
drukt, door te zorgen dat zij niet in de 
oxydelaag worden ingebouwd. Er is dan 
ook veel onderzoek gedaan naar fabricage 
van halfgeleiderelementen in een zeer 
zuivere omgeving, maar enige verontreini- 
ging is vaak niet te vermijden, gelet op het 
gewenste lage niveau van 1 beweeglijk ion 
op de 100.000 oppervlakte-atomen. 
Onderzoek naar de aard en het gedrag van 
de ionen kan ons in dit verband helpen 
om de ongewenste effecten te onder- 
drukken. 

Uit de hierboven beschreven BTS-experi- 
menten bleek niet alleen dat ionen onder 
invloed van een elektrisch veld door het 
oxyde bewegen maar ook dat een 
bepaalde hoeveelheid energie moet wor- 
den toegevoerd (elektrisch of thermisch) 
om ze van de grensvlakken vrij te krijgen. 
Kennelijk zijn er aan de grensvlakken 
plaatsen waar de ionen een minimale 
energie bezitten (invangcentra). 

Deze beide eigenschappen kunnen met 
afzonderlijke experimenten nader worden 
onderzocht. Overigens merken we op dat 
het experimenteel blijkt dat de ionen in 
het oxyde blijven; zowel het Si/SiOs- als 
het SiO, /metaal-grensvlak laten geen 
jonen door. Omdat we in de eerste plaats 
geinteresseerd zijn in de eigenschappen 
van het Si/SiO,-grensvlak, beschouwen 
we in het volgende de beweging van ionen 
vanaf het Si/SiO,-grensvlak naar het 

SiO /metaal-grensvlak toe (dus onder 


Fig. &. Het TSIC-experiment. De stroom 
ten gevolge van de beweging van de ionen 
I; als functie van de temperatuur’. 
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invloed van een negatieve elektrodespan- 
ning). 


De invangcentra aan het 
SifSiO,-grensvlak en de aard van de ionen 
Een methode om de energie te bepalen 
waarmee de ionen aan een grensvlak inge- 
vangen zitten is de thermisch gestimu- 
leerde ionenstroom (TSIC)-methode. 
Hierbij wordt na een BTS-experiment 
onder positieve elektrodespanning (alle 
ionen bevinden zich hierna aan het 
Si/SiO,-grensvlak) de temperatuur van de 
MOS-condensator langzaam verhoogd, 
van 0 tot 350 °C, met een constante 
opwarmsnelheid (O,1-1 °C/s). Er wordt 
een zo grote negatieve spanning op de 
elektrode gezet, dat elk ion dat uit het 
invangcentrum weet te ontsnappen direct 
naar de metaalelektrode wordt getrokken. 
De kans om te ontsnappen is sterk tem- 
peratuurafhankelijk en neemt toe met 
toenemende temperatuur. Aangezien voor 
elk ion dat het oxyde oversteekt, zich één 
elektron in het uitwendige circuit ver- 
plaatst, meten we de laadstroom van de 
condensator. Een TSIC-meting is getoond 
in figuur 8. In het bewuste temperatuur- 
gebied zien we twee stroompieken. Dat 
de stroom weer nul wordt betekent dat 
alle ionen zijn overgestoken. 


Wat de aard van de mobiele ionen betreft: 
algemeen wordt aangenomen dat 
Na'-ionen de voornaamste oorzaak van 
instabiliteiten zijn. Na'-ionen kunnen 
vrijkomen tijdens de fabricage (bijvoor- 
beeld uit de wand van de oxydatieoven of 
uit het schuitje vanwaaruit het metaal 
wordt opgedampt). Het blijkt dat de piek 
bij lage temperatuur in de TSIC-meting 
door Na* wordt veroorzaakt. Het ligt voor 
de hand voor de tweede piek aan een 
groter ion te denken (heeft meer ther- 
mische energie nodig om uit het invang- 
centrum te worden geëmitteerd). Zeer 
waarschijnlijk hebben we hier te maken 
met K*-ionen, zoals proeven met opzette- 
lijk ingebracht K* (uit een K* bevattende 
oplossing of door K*-implantatie) aan- 
toonden? ). Deze mening wordt echter 


Fig. 9. Schematisch model van invang- 
centra voor ionen aan het Si/SiO-grens- 
vlak. Als een ion is ontsnapt, wordt het 
direct door het elektrisch veld afgevoerd. 
(a) centra met één enkele energie E‚, (b) 
centra met spreiding in energie 

Eo -AESESE +AE. 












ionen energie 
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Si0, afstand SiO, afstand 
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niet door iedereen gedeeld. Het probleem 
is dat bij het opzettelijk aanbrengen van 
de ene verontreiniging, de andere vaak 
moeilijk is te vermijden. 


We kunnen ons de ionen in de invang- 
centra als volgt (figuur 9a) voorstellen: n, 
ionen zitten in de potentiaalkuil aan het 
grensvlak gevangen. Er iseen bepaalde _ 
kans dat zij een energie E, opnemen, 
waarmee ze uit de potentiaalkuil kunnen 
ontsnappen, waarna ze door het elek- 
trisch veld direct worden afgevoerd. Deze 
kans P(E) wordt beschreven met een 
Boltzmannfactor. Het aantal ontsnapte 
ionen per tijdeenheid hangt af van deze 
kans en van het aantal aanwezige ionen: 


ke P(E Zen 
= Ny. = Ny. S exPp(— 
df t ) t p( KT 


s heeft hier de dimensie van frequentie en 
wordt geassocieerd met de atomaire tril- 
lingsfrequentie (orde 10!2 Hz). Dit 
model blijkt de TSIC-curves heel aardig te 
verklaren. Een verfijning van het model 
bereiken we door aan te nemen dat niet 
alle ionen precies dezelfde ontsnappings- 
energie E‚, nodig hebben, maar dat daar 
een zekere spreiding AF in mogelijk is 
(figuur 9b). Het is mogelijk om uit de 
TSIC-curve de energieverdeling van de 
ionen in de invangcentra te bepalen. 


De beweeglijkheid van de Na*- en 
K'-ionen 

Bekeken we met de TSIC-methode het 
emissiegedrag, het is ook mogelijk om het 
transport van de ionen van de ene kant 
van het oxyde naar de andere kant te 
bekijken. Brengen we bijvoorbeeld na een 
BTS-experiment de MOS-condensator op 
een zo hoge temperatuur dat een groot 
aantal ionen voldoende energie heeft om 
uit de invangcentra te ontsnappen en 
variëren we langzaam de elektrodespan- 
ning, dan zullen de ionen het oxyde over- 
steken zodra het elektrisch veld van teken 
omkeert. Dit gebeurt in de 7VS(‘Trian- 
gular Voltage Sweep’)-methode. Voor de 
elektrodespanning wordt hier een lang- 
zaam verlopende driehoekspanning (hel- 
ling O,1-1 V/s) gebruikt. 


Bepalen we nu met de TSIC-meting die 
temperatuur, waarbij de Na'-ionen wel en 
de K*-ionen nog niet genoeg energie heb- 
ben bereikt om te ontsnappen, dan kun- 
nen we het transport van alleen de 
Na'-ionen onderzoeken met een TVS- 
meting op deze temperatuur (225°C, zie 
figuur 8). 


Om het transport van de ionen onder 
invloed van een elektrisch veld te kunnen 
beschrijven, moeten we allereerst hun 
beweeglijkheid u kennen. Beschouwen we 
het oxyde als een aaneenschakeling van 
potentiaalputjes met gemiddelde diepte W 
(figuur 10), dan kan de beweeglijkheid 
worden geschreven als 
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Fig. 10. Schematische voorstelling van 
de ionenbeweging in SiO,, onder invloed 
van een elektrisch veld. 


H=(HolT) exp [—W/KT]. 

Zo’n verband wordt experimenteel inder- 
daad gevonden. Het ion springt als het 
ware van het ene putje in het volgende. W 
is de activeringsenergie van de drift (bewe- 
ging onder invloed van een elektrisch veld). 
De beweeglijkheid van Na'- en K*-ionen 
kan op twee manieren worden gemeten. 
In de eerste plaats door bij een bepaalde 
temperatuur de elektrodespanning plotse- 
ling om te keren en de oversteektijd te 
meten van die ionen die genoeg energie 
hadden om niet ingevangen te zijn. Het 
nadeel van deze methode is dat bij lage 
temperatuur dit aantal ‘vrije’ ionen erg 
klein kan zijn. In de tweede plaats met 
behulp van de TVS-methode, waarbij de 
elektrode-spanning langzaam wordt 
gevarieerd. Zo’n TVS-meting is in 

figuur 11 weergegeven. Door de hoge 
temperatuur hebben zowel de Na'- als de 
K*-ionen voldoende energie om uit de 
invangcentra te ontsnappen. De beweeg- 
lijkheid van Na'-ionen is zo groot dat zij 
het oxyde direct oversteken als het elek- 
trisch veld in het oxyde van teken 
omkeert (bij VG = 0,1 V). In de TVS- 
meting veroorzaken zij een stroompiek bij 
die spanning. De K*-ionen beginnen bij 
deze spanning ook aan hun overtocht. 
Hun beweeglijkheid is lager en het duurt 
enige tijd voor zij zijn overgestoken. Hun 
snelheid wordt wel steeds groter omdat 
de elektrodespanning en dus het elek- 
trisch veld in het oxyde steeds toeneemt. 
Als de eerste K*-ionen de overkant berei- 
ken is de stroom maximaal. Uit de plaats 


van dit maximum als functie van tempera- 


tuur en helling van de langzaam verlo- 
pende driehoekspanning kan de beweeg- 
lijkheid worden bepaald* . Een overzicht 
van recente bepalingen van de beweeglijk- 
heid van Na'- en K*-ionen volgens beide 
methoden is weergegeven in figuur 12. De 
resultaten sluiten goed bij elkaar aan. 


Beweeglijke ionen in chlooroxyden 

Een nieuwe ontwikkeling in de fabricage 
van MOSFET’s is de toepassing van 
chloorverbindingen (HCI of 

1.1.1-C‚ H3Cl3 ), niet alleen om de oven- 
buis te reinigen, maar ook als toevoeging 
tijdens de oxydatie. Mits dit niet in te 
grote hoeveelheden gebeurt, treden heil- 
zame effecten op: het aantal defecten 
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waaruit de beweeglijkheid van K*-ionen op Na’. Het blijkt dat het chloor in het 
kan worden berekend“. De laadstroom 
ten gevolge van de capaciteit van het 


oxyde is onderdrukt. Fig. 13. Defectdichtheid in het oxyde 


als functie van de HCI-toevoeging tijdens 


„6 
wordt minder (figuur 13) en de veroude- de oxydatie”. 


ring wordt minder. Of het aantal opper- 50 
vlaktetoestanden groter of kleiner wordt 
hangt sterk af van de oxydatie-omstandig- 
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oxyde wordt ingebouwd dicht bij het 
Si/SiO, -grensvlak. Worden nu de Na'- 
ionen in een TVS-experiment bij 225 °C 
naar het Si/SiO, -grensvlak gedreven (met 
een positieve elektrodespanning) dan 
komen de ionen als de spanning van teken 
wordt omgedraaid niet meer terug! 
Bovendien blijkt uit een vlakkebanden- 
spannings-meting dat het effect dat 
Na'-ionen aan het Si/SiO, -grensvlak op de 
vlakkebandenspanning hebben geheel 
verdwenen is: de Na'-ionen hebben ken- 
nelijk hun lading verloren, ze zijn geneu- 
traliseerd. Dat dit effect niet definitief is, 
zien we aan TSIC-metingen. Deze metin- 
gen worden uitgevoerd na een zodanig 
BTS-experiment dat de Na'-ionen naar 
het Si/SiO, -grensvlak zijn gebracht. In 
figuur 14 worden TSIC-metingen getoond 
aan MOS-condensators met en zonder 
chloor in het oxyde. Beide oxydes zijn 
opzettelijk met Na'-ionen verontreinigd. 


In de TSIC-curve van de MOS-conden- 
sator zonder chloor in het oxyde, herken- 
nen we de stroompiek ten gevolge van 
‘vrije’ Na'-ionen, die ook te zien was in 
figuur 8. Dat de plaats, vorm en grootte 
van deze Na*-piek in beide figuren niet 
gelijk is, komt door verschillen in opwar- 
ming, veld in het oxyde en aantal Na*- 
ionen. Van het stroompiekje bij 250 °C 
vermoeden we dat het K*-ionen zijn, inge- 
bracht tijdens de metallisatie. Hun aantal 
is echter veel kleiner dan het aantal Na*- 
ionen. 


In de TSIC-curve van de MOS-conden- 
sator met chlooroxyde ontbreekt de bij- 
drage van de ‘vrije’ Na'-ionen. Pas bij ruim 
300 °C komen de geneutraliseerde Na'- 
ionen weer in ionvorm vrij en steken het 
oxyde over. Deze hoge temperatuur bete- 
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kent dat de invangcentra diep zijn, dat wil | Fig. 14. TSIC-experiment aan MOS-con- 


zeggen: er is een grote energie E‚ nodig 
om de ionen uit de centra te laten ont- 
snappen. 


De conclusie is, dat bij toevoeging van 
chloor aan oxyden de invangcentra voor 
Na'-ionen aan het Si/SiO»-grensvlak zo 
worden gewijzigd, dat het Na*-ion dieper 
wordt ingevangen en bovendien zijn 
lading verliest. Het mechanisme van deze 
vorm van invangen is nog een onbegrepen 
zaak. Omdat het chloor aan het grensvlak 
wordt ingebouwd, is de beweeglijkheid 
van Na'-ionen in het oxyde onafhankelijk. 
van de toevoeging van chloor. | 
Dat er voor K*-ionen in chlooroxyden ook! 
merkwaardige effecten optreden, blijkt 
uit TVS-metingen. Wat er aan de hand is 
zal nog nader moeten worden onder- 
zocht. > 


densators met oxyden gegroeid in Oy 
(curve a) en in 3% C‚H3 Cl3/0Oz (curve b). 


|-Na de oxydatie zijn de oxyden met Na*- 


ionen uit een Na' bevattende oplossing 
verontreinigd. De ionen bewegen vanaf 
het Si/SiO-grensvlak naar het metaal/ 
SiO,-grensvlak toe. 
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Allereerst wil ik trachten de twee grond- 
gedachten van de moderne differentiaal- 
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meetkunde te isoleren. Daarmee heb ik 
dan tevens het gemeenschappelijk funda- 
ment van de thans te bespreken werken 
aangegeven. 


Enerzijds wenst men het begrip inhoud 
(volumen) zodanig te formuleren dat het 
geldig is voor willekeurige n-dimensionale 
ruimten. Het is bekend dat in IR” de de- 
terminant functie det(v;,..., Vn) op het 
teken na het volumen geeft van het paral- 
lelepipidum dat opgespannen wordt door 
de vectoren vi ,..., vn. De determinant- 
functie is lineair in elk van zijn argumen- 
ten, en bovendien is ze anti-symmetrisch. 
Men zegt dat de determinant een alter- 
nerende n-vorm is. Indien men nu iets 
analoog wil doen voor k vectoren in een 
n-dimensionale ruimte (k <n) dient men 
het begrip alternerende k-vorm te intro- 
duceren. 


Maar behalve de notie van volumen wil 
men òok het begrip oppervlak generali- 


gelijkt op wat een fysicus zich traditio- 


seren. Intuitief is een cirkel een één- 
dimensionaal oppervlak in IR? en een bol 
een twee-dimensionaal oppervlak in IR*. 
De vraag is natuurlijk: hoe parametriseren 
wij deze oppervlakken? Dit gebeurt mid- 
dels het begrip deelvariëteit. De cirkel is 
een één-dimensionale deelvariëteit van 
R?,en in het algemeen spreekt men van 
een m-dimensionale deelvariëteit in IR” 
(m <n). Een deelvariëteit beschrijft men 
met behulp van ‘kaarten’, die dan in een. 
‘atlas’ worden verenigd. Een nette deel- 
variëteit zal men tevens willen voorzien 
van een rand en een oriëntatie. Indien 
zulks mogelijk is verkrijgt men een fraaie 
synthese van de begrippen alternerende 
k-vorm en (k + 1)-dimensionale deelva- 
riëteit met oriëntatie en rand via de stel- 
ling van Stokes. 


De lezer dient zich evenwel te realiseren 
dat de moderne differentiaalmeetkunde 
op het eerste gezicht nog maar in weinig 
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neel had voorgesteld, en ze bovendien een 
zwaar en zeer formeel (zo niet abstract) 
begripsapparaat benodigt. Kortom, men 
dient zich te realiseren dat dit vak niet 
eenvoudig is, doch tevens dat het in zijn 
huidige vorm stormachtige ontwikkelin- 
gen doormaakt. Wil men het vak voor fy- 
sici tot leven laten komen dan is een goe- 
de, uitgebreide, en zorgvuldig motiveren- 
de inleiding van het allergrootste belang. 
Voor mathematici is zulk een inleiding 
gegeven door Spivak [1, 2], maar tot op 
heden ontbrak het aan literatuur [3] die 
duidelijk gericht was op fysische proble- 
men. In deze leemte hebben nu Thirring 
en Arnold willen voorzien. Bij een ge- 
slaagde opzet zouden we hun daarvoor 
zeer erkentelijk moeten zijn. Maar we die- 
nen wel te bedenken dat hun boeken be- 
doeld zijn als leerboeken, niet als weten- 
schappelijke monografieën. 


De klassieke mechanica met haar ‘neven- 
voorwaarden’ en corresponderende gege- 
neraliseerde coördinaten geeft op natuur- 
lijke wijze aanleiding tot de introductie 
van deelvariëteiten, raakruimten, en meer 
van dat schoons dat dient tot opluistering 
van de bewegingsvergelijkingen. Ook het 
begrip inhoud komt op natuurlijke wijze 


naar voren via het theorema van Liouville. 


Een passend kader ter introductie van de 
nieuwe begrippen is hier dus duidelijk 
aanwezig. 


Thirring’s boek ‘Klassische dynamische 
Systeme’ bestaat uit drie gedeelten: een 
derde is gewijd aan het mathematische 
formalisme, een derde aan Hamilton sys- 
temen (met een bewijs van het KAM- 
theorema!) en niet-relativistische meer- 
deeltjes problemen, en het laatste deel 
aan de relativiteitstheorie. De eruditie van 
de schrijver is imponerend, de veelheid 
van problemen is fascinerend, maar de 
presentatie lijkt mij volstrekt ongeschikt 
om dit boek als leerboek te laten functio- 
neren. We zullen ter staving van deze be- 
werking eens het eerste deel (het mathe- 
matische formalisme) bezien dat absoluut 
essentieel is voor al wat volgt. Niet alleen 
is de stijl een verschrikking (praktisch 
iedere zin is genummerd) maar ook ont- 
breekt elke motivatie. Het geheel is een 
dorre opeenvolging van definities, lem- 
ma’s, stellingen en opmerkingen, en er 
wordt een algemeenheid en een abstractie 
nagestreefd die, me dunkt, zelfs Bourbaki 
tot schuchterheid en bescheidenheid zou- 
den hebben gestemd. Ook de formulering 
is veel te kort en te strak gehouden (eg. 
het lemma van Poincaré is in een vraag- 
stuk ondergebracht). De stijl en opzet van 
het ‘samenvattende’ collegedictaat waar- 
bij de collegegever een essentiële rol 
speelt, is overal nog speurbaar; te licht- 
vaardig is deze stijl in het boek gehand- 
haafd. 


Ik heb mij afgevraagd hoe dit boek in de 
nederlandse context zou kunnen functio- 
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neren. Thirring zelf geeft een suggestie: 
hij werkt ‘auf einführende Vorlesungen 
aufbauend’. Goed, ik heb eens een inlei- 
dend (prekandidaten) collegedictaat [4] 
ter hand genomen en nog eens te meer 
ontdekt wat aan Thirring’s boek ont- 
breekt: leesbaarheid. Hoe dan ook, zelfs 
met een inleidend college à la [4] is, vrees 
ik, de nederlandse student nog onvol- 
doende toegerust om het mathematisch 
formalisme snel en doeltreffend te kun- 
nen verwerken. Ik kan derhalve ‘Klassi- 
sche dynamische Systeme’ alleen aan die 
docenten (en studenten) in overweging 
geven die èn zeer gecharmeerd zijn van 
een formele, ‘mathematisch lückenlose’ 
opbouw, èn over een zee van tijd beschik- 
ken, èn de steun hebben van een vakgroep 
Wiskunde die het vele, wiskundig beno- 
digde grondwerk willen verrichten. 


Thirring’s tweede recente boek, ‘Klassi- 
sche Feldtheorie’, is in bijna dezelfde stijl 
geschreven als het eerste deel. Na een in- 
leidend hoofdstuk komen aan de orde het 
elektromagnetische veld van gegeven la- 
dingsverdelingen, het veld in de aanwezig- 
heid van een geleider, en tenslotte gravita- 
tie. Beide delen bevatten vele opgaven 
met de bijbehorende oplossingen. 


Alle vier delen van de door Thirring ge- 
plande reeks over Mathematische Fysica 
zullen in het engels worden vertaald; het 
eerste deel is reeds verschenen, eveneens 
bij Springer (nu Berlin, etc.). Voor slechts 
weinige guldens meer krijgt men dan een 
gebonden boek, een betere typografie en 
een iets betere formulering. Ik betreur 
echter dat hierdoor deze boeken worden 
gefixeerd in een vorm die nog te zeer 
tijdsgebonden en daardoor nog te weinig 
definitief is. Maar wellicht is dat een ken- 
merk van baanbrekend werk. 


VI. Arnold’s boek, Mathematical 
methods of classical mechanics’, ademt 
een geheel andere sfeer. Ook hier moder- 
ne differentiaalmeetkunde, maar wat een 
verschil tussen de mathematicus Arnold 
en de fysicus Thirring! Arnold heeft zijn 
boek zeer zorgvuldig geschreven [5,7]. 
Hij illustreert zijn definities en bewerin- 
gen uitgebreid aan de hand van voorbeel- 
den en figuren en tracht een en ander zo- 
veel mogelijk te visualiseren. Het boek is 
met grote zorg opgebouwd (het heeft een 
‘gelaagde’ structuur) zó dat de lezer 
steeds goed gemotiveerd blijft. En een 
stelling wordt vaak niet in haar abstract- 
ste en algemeenste formulering bewezen 
maar juist zo (bijv. door een eenvoudig 
geval te bekijken) dat de bewering begrij- 
pelijk en doorzichtig wordt. Toegegeven, 
ook hier wordt af en toe enige zelf werk- 
zaamheid verlangd, maar in de praktijk is 
reeds gebleken dat er met dit boek uitste- 
kend te werken valt. Het laat zich ook 
voortreffelijk als aanvulling op Goldstein 
[6] gebruiken: begrippen als Legendre- 
transformatie, canonieke transformatie, 


integraalinvariant van Poincaré-Cartan, of 
‘the principle of least action’ vindt men 
bij Arnold op heldere en, zo men wil, 
alternatieve wijze behandeld. Dit boek be- 
staat eveneens uit drie gedeelten: de me- 
chanica volgens Newton, volgens Lagrange 
(introductie van deelvariëteiten), en vol- 
gens Hamilton (introductie van alterne- 
rende vormen, etc). In een uitgebreide 
reeks appendices worden tenslotte geavan- 
ceerde onderwerpen behandeld. Ik zou 
(theoretisch) fysici dit boek ten zeerste 
willen aanbevelen. 
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Beeldopnemer met 


weerstandselektroden 


M.G. Collet 


Het vastleggen van bewegende beelden blijft veel mensen boeien. Het ligt voor de hand 
dat de elektronische industrie zou proberen met elektronische middelen optische 
beelden weer te geven. Dit heeft geresulteerd in televisie en videorecording. De volgen- 
de stap is de ontwikkeling van een handzame goedkope TV-camera. Een essentieel on- 
derdeel hiervoor is een vastestofbeeldopnemer. 

Sinds de opkomst van de geïntegreerde circuits is gezocht naar een circuit waarmee een 
beeld omgezet kan worden in een elektronisch signaal. Dat circuit zou veel goedkoper 
moeten zijn dan de opneembuizen gebruikt in studiocamera’s. Het zou bovendien veel 
kleiner en betrouwbaarder kunnen zijn. Nieuwe principes (emmertjesgeheugen' ), 
ladingsgekoppelde inrichting* en de versnelde technologische ontwikkeling (‘very 
large scale integration’) hebben de vastestofbeeldopnemer en daarmee een goedkope 


TV-camera bijna mogelijk gemaakt. 


In deze bijdrage worden kort enkele ontwerpen van beeldopnemers geschetst. Daarna 
volgt een iets uitgebreider beschrijving van de beeldopnemer met weerstandselektroden 
waaraan op ons laboratorium gewerkt wordt. 


Integratie van fotostroom, lekstroom 
Door het aanleggen van een positieve 
spanning aan het oppervlak van een p-ge- 
leidend siliciumplaatje kunnen de gaten 
van het oppervlak worden weggeduwd. 
Zo ontstaat een laagje waaruit alle gaten 
verwijderd zijn. Deze depletielaag heeft 
een dikte van enkele microns (10° m). 
Hierin heerst een veld dat elektronen de 
laag intrekt en gaten wegduwt. Wanneer 
door invallend licht gat-elektronparen ge- 
genereerd worden zal het veld de gaten en 
elektronen van elkaar scheiden. Zo ont- 
staat een fotostroom. 

In een vastestofbeeldopnemer met de re- 
solutie van een TV-beeld moet de foto- 
stroom op ongeveer 250.000 plaatsen ge- 


Fig. Ll. Matrix van metaal-oxyde-silici- 
um-condensatoren en het potentiaalpro- 
fiel dat bij het siliciumoppervlak wordt 
geïnduceerd. a: metaalvlekje; b: silici- 
umdioxyde;c: p-type silicium. 
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meten worden. Door de fijne onderverde- 
ling wordt de fotostroom per beeldpunt 
zeer klein. Daarom wordt de fotostroom 
in ieder beeldpunt enige tijd geintegreerd. 
Zo ontstaan ladingspakketjes die door het 
silicium of via geleidersporen op het silici- 
um naar een uitgangsversterker worden 
getransporteerd. Voor het integreren van 
de fotostroom zijn — zeer kleine — con- 
densatoren nodig. Een hiervoor veel ge- 
bruikte structuur (figuur 1) bestaat uit 
metaalvlekjes op een isolerende silicium- 
dioxyde huid op het p-type kristal. Een 
positieve spanning op de vlekjes zorgt 
voor lokale depletiegebieden. Hierin wor- 
den de gegenereerde elektron-gatparen 
van elkaar gescheiden. De foto-elektronen 
worden verzameld op de vlakke-plaatcon- 
densator gevormd door de metaal-oxyde- 
silicium (MOS)-structuur. Het transport 
van de ladingspakketjes naar de uitgang 
vormt een hoofdstuk apart. Gat-elektron- 
paren kunnen behalve door licht ook 
door thermische excitatie worden gegene- 
reerd. Deze thermische generatie wordt 
door bepaalde verontreinigingen en kris- 
talfouten enorm versneld. Om te voor- 
komen dat deze lekstroom het beeld be- 
derft moeten zeer hoge eisen aan de sili- 
ciumtechnologie worden gesteld. Op 

102 siliciumatomen mag er slechts één 
ongewenst verontreinigingsatoom aanwe- 
zig zijn. De gemiddelde concentratie kris- 
talfouten mag niet veel groter zijn dan 
één per vierkante centimeter. 


Een uiterst summiere schets van drie door 
veel laboratoria intensief onderzochte 
beeldopnemers 

De X-Y opnemer is het oudste type op- 
neemcircuit (figuur 2). Het maakt gebruik 
van een schakelmatrix. Elk beeldelement 
wordt op zijn beurt galvanisch verbonden 





Ir. M.G. Collet (38) is natuurkundig inge- 
nieur (Delft, 1966). Sinds zijn afstuderen 
werkt hij op het natuurkundig laborato- 
rium van Philips. Tot 1972 bestudeerde 
hij het siliciumoxyde-siliciumgrensvlak en 
de eigenschappen van kristalfouten in si- 
licium.. 

Daarna hield hij zich bezig met de ‘device’ 
fysica van ladingsoverdracht-inrichtingen 
(CTD'’s). Vanaf 1977 maakt hij deel uit 
van een werkgroep die het realiseren van 
een vastestofbeeldopnemer tot taak heeft. 





met een verticale Y-lijn en een horizonta- 
le X-lijn en via deze weg met een uitgangs- 
versterker. Voor de schakelaars worden 
MOS-transistoren gebruikt. Het grote 
voordeel van deze X-Y-opnemer is dat hij 
in een al op grote schaal toegepaste tech- 
nologie kan worden gemaakt. Met de 
X-Y-opnemer worden op het moment dan 
ook de beste beelden gemaakt”. Het gro- 
te nadeel is de capaciteit van de verticale 
en horizontale geleidersporen. Deze is 
veel groter dan de capaciteit van het 
beeldelement. Het kleine pakket foto- 
elektronen geeft maar een heel kleine po- 
tentiaalverandering op de geleiders. Ten- 
gevolge van deze verzwakking is dit soort 
opnemer erg ongevoelig. | 


Dit bezwaar kan worden ondervangen 
door gebruik te maken van ‘charge | 
coupled devices’ (CCD). Een CCD (figuur 
3a) kan bestaan uit een n-type laag op een 
p-type kristal. De n-type laag is bedekt 
met een dunne oxydelaag. Hierop zijn 
elektroden aangebracht. De n-laag is ZO 
dun en licht gedoteerd dat alle elektronen 


Fig. 2. Basisschema van X-Y-beeldop- 
nemer. 
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X b. 
Fig. 3a. Dwarsdoorsnede door CCD. a: 
elektrode van polykristallijn silicium; b: 
siliciumdioxyde; c: n-type silicium; d: 
| p-type silicium. 
b. Potentiaalverloop en ladingsover- 
dracht in de n-laag van figuur 3a. 


door het aanleggen van een positieve 
spanning aan deze ndaag kunnen worden 
verwijderd. Daarna worden de elektroden 
op verschillende spanningen gebracht. 
Hierdoor ontstaan in de n-laag potentiaal- 
maxima, potentiaalputten voor elektro- 
nen. In de potentiaalputten kunnen pak- 
ketjes elektronen, geïsoleerd van elkaar, 
worden opgeslagen. Als elke CCD-cel drie 
elektroden heeft kan één elektrode (op 
negatieve potentiaal) de pakketjes ge- 
scheiden houden. Een positieve elektrode 
zorgt voor een ontvangende potentiaal- 
put, een van positief naar negatief gescha- 
kelde stuwt het pakket voort (figuur 3b). 
Door dit proces van elektrode tot elektro- 
de te herhalen kunnen pakketjes elek- 
tronen door een silictumlaag worden ge- 
transporteerd. De verliezen zijn laag, 
minder dan een duizendste procent per 
stap, het transport is bovendien snel. 

In de zogenaamde “interline transfer” 
beeldopnemers (figuur 4) bevindt zich 


Fig. 4. ‘Interline transfer’ beeldopnemer. 
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naast elke kolom beeldelementen een 
CCD-transportregister. Gedurende de in- 
tegratietijd wordt in elk beeldelement een 
pakket foto-elektronen gevormd. Aan het 
eind van de integratietijd worden alle la- 
dingspakketten in de verticale CCD-regis- 
ters geplaatst. Deze sluiten aan op een 
gemeenschappelijk horizontaal CCD. Dit 
transporteert telkens één lijn informatie 
naar de uitgangsversterker en ontvangt 
dan weer een nieuwe lijn uit de verticale 
registers. In de tijd nodig voor het uitle- 
zen van de registers worden in de beeld- 
punten weer nieuwe ladingspakketten ge- 
vormd. Tijdens het transport mogen de 
ladingspakketten niet meer door opval- 
lend licht worden vergroot. Daarom wor- 
den de verticale registers bedekt met on- 
doorzichtig materiaal. Dit betekent dat 
een flink deel van het opvallende licht 
niet wordt gebruikt. Bovendien wordt de 
horizontale afstand tussen de beeldele- 
menten vergroot door de ruimte vereist 
voor de verticale registers. Dit beperkt 
het scheidend vermogen dat met een cir- 
cuit van gegeven afmetingen bereikbaar 
is. 


In de ‘frame transfer’ opnemer® wordt 

de fotostroom geintegreerd in de poten- 
tiaalputten van een CCD-register (figuur 
5). Het lichtgevoelige deel van dit circuit 
bestaat uit parallelle CCD-registers. Gedu- 
rende de integratietijd worden de elek- 
trodepotentialen constant gehouden. Ver- 
volgens wordt het hele ladingspatroon 
snel naar een buffergeheugen getranspor- 
teerd. Ook dit buffergeheugen bestaat uit 
parallelle CCD-registers maar afgeschermd 
tegen licht. Het aantal cellen per register 
moet hier even groot zijn als in het licht- 
gevoelige deel. Maar hét oppervlak per cel 
kan kleiner zijn, omdat de opslag van la- 


Fig. 5. ‘Frame transfer’ beeldopnemers. 


1: lichtgevoelig gedeelte; 
2: buffergeheugen. 
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Fig. 6a. Dwarsdoorsnede door transport- 
kanaal met weerstandselektrode. a: hoog- 
ohmige elektrode; b: siliciumdioxyde; c: 
n-type silicium; d: p-type silicium. 

b. Potentiaalverloop en ladingsover- 
dracht in de n-laag van figuur Óa. 


ding minder oppervlak vraagt dan de ge- 
neratie. De uitlezing van het buffergeheu- 
gen is analoog aan de uitlezing van de 
‘interline transfer’ beeldopnemers. 

Voor het lichtgevoelige deel van de frame 
transfer-schakeling zijn transparante elek- 
troden nodig. Polykristallijn silicium vol- 
doet heel goed als elektrodemateriaal en 
laat groen en rood licht door. Een kleu- 
rencamera moet ook gevoelig zijn voor 
blauw licht. Indium- en tinoxyde lijken 
hiervoor bruikbaar maar kunnen nog niet 
in de siliciumtechnologie worden inge- 
past. 


Beeldopnemer met weerstandselektroden 
We willen een gevoelige beeldopnemer, 
dus geen X-Y-uitlezing. We willen ook 
blauw licht detecteren, dus niet het ge- 
heel met elektroden bedekte oppervlak 
van frame transfer. Resteert de interline 
transfer-organisatie, maar liefst met smal- 
le banen voor het verticale transport. 

Een nieuw transportmechanisme lijkt 
hiervoor mogelijkheden te geven (figuur 
6): een CCD-structuur met slechts één 
elektrode. De elektrode is gemaakt van 
hoogohmig polykristallijn silicium, zodat 
een flink spanningsverschil tussen de uit- 
einden weinig dissipatie geeft. Net als in 
het CCD is de n-laag en een deel van het 
onderliggende p-materiaal gedepleerd. Het 
spanningsverschil langs de elektrode indu- 


ceert een gradiënt in de n-laag (figuur 6b). 


Rekening houdend met de vertragende ef- 
fecten van diffusie en coulombafstoting 
kan men berekenen dat een elektronen- 
pakket in 25us langs een 0,4 cm lang ka- 
naal kan worden getransporteerd met een 
potentiaalverschil van 10 V. Het tijdsche- 
ma gebruikt bij televisie- en video- 
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Fig. 7a. Dwarsdoorsnede door beeldele- 
ment en transportkanaal van een beeld- 
opnemer met weerstandselektroden. 

a: integratie-elektrode; b: weerstands- 
elektrode; c: n-type silicium; d: zwaar ge- 
doteerd p-type silicium; e: p-type silici- 
um. 

b. Potentiaalverdeling tijdens integratie. 
c. Potentiaalverdeling bij ladingsover- 
dracht. 


recording staat 64 us toe. Dit mechanisme 
is dus voldoende snel voor het verticale 
transport van de ladingspakketten, mits 
we voldoende beeldlijnen op 0,4 cm kun- 
nen plaatsen. 

De koppeling van de beeldelementen met 
het transportkanaal is geschetst in figuur 
7. In één n-gebied worden beeldelement 
en transportkanaal gemaakt. Ze zijn van 
elkaar gescheiden door een ondiep zwaar 
gedoteerd p-gebied. De potentiaal in het 
gedepleerde n-gebied onder de p-baan 
wordt door de vaste potentiaal van de p- 


Fig. &. Beeldopnemer met weerstands- 
elektroden. 1: digitaal register; 2: verticale 
transportbaan; 3: integratie-elektrode; 

4: bufferelektrode; 5: CCD-register; 6: 
uitgangsversterker. 
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baan op een kleine positieve waarde ge- 
houden. De potentiaal van het integratie- 
gebied wordt bepaald door de integratie- 
elektrode. Tijdens de integratie is deze 
flink positief zodat een elektronenpakket 
kan worden verzameld. Daarna wordt de 
elektrode naar nul geschakeld. De poten- 
tiaalput wordt ondieper, de elektronen 
worden geïnjecteerd in het transportka- 
naal (figuur 7c). De andere begrenzing 
van het integratiegebied wordt gevormd 
door een zwaar gedoteerd p-gebied in di- 
rect contact met het p-kristal. De poten- 
tiaalbarrière aan die kant is zo hoog dat 
de elektronen alleen in het ‘eigen’ trans- 
portkanaal geïnjecteerd worden. 

Met deze elementen hebben we de beeld- 
opnemer met weerstandselektroden ont- 
worpen (figuur 8)°. De transportkanalen 
lopen verticaal, de integratie-elektroden 
horizontaal. Met een digitaal schuifre- 
gister kunnen de integratie-elektroden op 
eenvoudige wijze sequentieel naar nul volt 
worden geschakeld. De geïnjecteerde la- 
dingspakketten worden na transport on- 
der een bufferelektrode verzameld. Tij- 
dens het transport wordt een voorgaande 








Fig. 9. Microscoopopname van een ge- 
deelte van de beeldopnemer met weer- 


standselektroden. 


lijn uitgelezen door het horizontale CCD- 
register. Als dit register leeg is geeft de 
bufferelektrode een nieuwe lijn over aan 
het CCD. Zo kan de hele lijntijd worden 
gebruikt voor het verticale transport. 

De elektroden laten een flink deel (25%) 
van het oppervlak onbedekt. Dit type 
sensor heeft dus een behoorlijke gevoelig- 
heid voor blauw licht. 


Bereikte resultaten 

Om onze ideeën te testen hebben we een 
opnemer ontworpen met 300 lijnen en 
200 kolommen. Het beeldformaat is 0,42 
cm bij 0,56 cm. Het totale circuit meet 
0,55 em bij 0,75 cm. 

Een microscoopopname (figuur 9) toont 
rechts de verticale transportkanalen met 





Fig. 10. Beeld opgenomen met het cir- 
cuit van figuur 9. 





loodrecht daarop de integratie-elektro- 
den, links daarvan het digitale schuifregis- 
ter voor de lijn-aansturing. Onder de: 
transportkanalen het CCD-register, links 
aansluitend op de uitgangsversterker. De 
bereikte beeldkwaliteit is behoorlijk (fi- 
guur 10). Met 50 000 elektronen per sig- 
naalpakket wordt een signaal-ruisverhou- 
ding van 40 dB bereikt. Om elektronen- 
pakketjes van deze grootte te genereren 
is, bij gebruik van een F2-lens, een ver- 
lichtingssterkte van 120 lux voldoende. 
Dit komt overeen met het lichtniveau in 
een huiskamer. 

Voor voldoende horizontale resolutie is . 
een groter aantal beeldpunten per lijn no- 
dig. Wanneer we bovendien kleur- en hel- 
derheidsinformatie met één opnemer wil- 
len verkrijgen, dan moeten we naar 600 
lijnen van elk 400 kolommen. Dat vraagt 
grotere circuitafmetingen en dus moeten 
de verticale transportkanalen met hun 
weerstandselektrode ook langer worden. 
Om na te gaan hoever we kunnen gaan 
hebben we de transportsnelheid van de 
huidige kanalen gemeten. Daartoe injec- 
teren we de bovenste lijn en meten we de 
stroom aan het eind van de 0,42 cm lange 
kanalen als functie van de tijd (figuur 11). 
De getrokken lijn geeft het theoretisch 
verwachte resultaat voor een aangelegd 
potentiaalverschil van 10 V en een elek- 
tronenbeweeglijkheid van 750 cm°-V”! 
“st, De onderbroken lijnen geven de 
meetresultaten voor toenemende grootte 
van de ladingspakketten. Voor grote pak- 
ketten komt ongeveer de helft van de la- 
ding met de verwachte tijdsverdeling naar 
buiten. De rest van de lading is verdeeld 
over een relatief lange staart. De staart 
verzadigt met toenemende pakketgrootte. 
Het lijkt erop dat een deel van de elek- 
tronen wordt ingevangen en met varië- 
rende tijdconstante weer wordt vrijgelaten. 
Het aantal invangplaatsen is beperkt, ge- 


L-42mm 
is _AV=10V 
H=750 cm2/ vs 





t lus) 
Fig. 11. De uitgangsstroom (in arbitraire 
eenheden) van de verticale transportkana- 
len als functie van de tijd na injectie van 
de bovenste rij ladingspakketten. 


zien de verzadiging van de staart met toe- 
nemende pakketgrootte. 

Het is bekend dat in silicium verontreini- 
gingen aanleiding kunnen geven tot gelo- 
kaliseerde energieniveaus tussen gelei- 
dings- en valentieband. Dergelijke niveaus 
kunnen elektronen invangen en — door 
thermische excitatie — weer vrijgeven. 
Om de waargenomen effecten te verkla- 
ren zijn ongeveer 5 x 10 '? van deze cen- 


tra per kubieke centimeter silicium nodig. 


Het is ook mogelijk dat variaties in do- 
teringsconcentratie of afmetingen van het 
kanaal kuilen veroorzaken in het poten- 
tiaalverloop langs het kanaal. Als daar- 
door gelokaliseerde potentiaalputjes ont- 
staan, dan kunnen daarin elektronen 
worden ingevangen en door thermische 
excitatie weer worden vrijgemaakt. Nader 


onderzoek zal ons moeten leren welk van 


deze effecten ons parten speelt en welke 
tegenmaatregelen kunnen helpen. Pas 
daarna zullen we de beeldopnemer signifi- 
cant kunnen vergroten. 


Slotopmerkingen 

Dit verhaal laat iets zien van de mogelijk- 
heden voor het ontwerp van vastestof- 
beeldopnemers. Daarbij ontbrak de ruim- 
te om op deelproblemen in te gaan: hoe 
verdelen de foto-elektronen zich over de 
verschillende beeldelementen, hoe voor- 
kom je het overlopen van overbelichte 
elementen, welke kleurfilterontwerpen 
zijn optimaal, welke bewerkingen moet je 
nog op het uitgangssignaal toepassen, hoe 
kun je de optimale doteringsconcentraties 
en laagdiktes berekenen. De opgegeven 
litteratuurplaatsen geven daar gedetail- 
leerde informatie over. 

We moeten hier volstaan met de conclusie 
dat in verschillende laboratoria is aange- 
toond dat met vastestofbeeldopnemers 
goede kleurencamera’s gemaakt kunnen 
worden. Na verdere perfectionering van 
de silicitumtechnologie zal het ook moge- 
lijk worden de opnemers foutvrij en goed- 
koop te maken. 
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De ontwikkeling van ionengevoelige 
veldeffecttransistoren (ISFET) op basis 
van de MOS-IC-technologie 


P. Bergveld en N.F. de Rooij 


De ISFET (ion-sensitive field effect transistor) is een ionengevoelige sensor, ontwik- 
keld door integratie van de ionengevoelige glaselektrode en de MOSFET (metal-oxide- 
silicon field effect transistor). Werking, vervaardiging en toepassing van de ISFET wor- 


den beschreven. 


Historie 

Reeds uit de jaren twintig stamt het idee 
de elektrische geleiding van een plaatje 
_halfgeleidermateriaal te veranderen met 
een elektrisch veld, loodrecht op het 
plaatje. Het zou daardoor mogelijk zijn de 
22 


elektrische stroom energieloos te sturen 
in een kleine robuuste structuur, zonder 
gebruikmaking van een gloeidraad 


Uitvoerige metingen van Bardeen en Brat- 
tain van Bell Telephone Labs. in de 


periode 1947-50 aan germanium-éénkris- 
tallen leidden niet tot het genoemde veld- 
effect. De oorzaak lag in het feit, dat de 
geïnduceerde ladingsdragers immobiel 
werden door vangst in vrije valentiebin- 
dingen: de oppervlaktetoestanden. Het 
bewijs, dat de oppervlaktetoestanden 
kunnen worden geneutraliseerd, leverden 
Garrett en Brattain bij metingen aan een 
germanium/elektrolytcel. Ook bleek dat 
onder kathodische voorspanning de gelei- 
ding van het germanium aanzienlijk kon 


worden beïnvloed. Praktische betekenis 


hadden deze experimenten niet. Wel 
brachten deze onderzoeken Shockley tot 
de uitvinding van de pn-overgang (1949) 
en de junctietransistor (1950). 


In de jaren die volgden onderging de junc- 
tietransistor een grote ontwikkeling. De 
oppervlaktetoestanden bleven echter van 
storende invloed op de eigenschappen van 
de transistor. Veel werk werd gestoken in 
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Dr. ir. P. Bergveld (39) studeerde elektro- 
techniek aan de TH-Eindhoven. In 1965 
trad hij in dienst bij de Technische Hoge- 
school Twente als medewerker van de 
vakgroep Bio-informatica. Hij heeft zich 
daar gespecialiseerd in het ontwikkelen 
van meet- en stimulatiemethoden ten 
behoeve van de biomedische research. 
Zijn speciale aandacht heeft de koppeling 
tussen de “ionenwereld’ en de ‘elektronen- 
wereld’ onder andere via de in dit artikel 
beschreven ionengevoelige veldeffect-tran- 
sistor. Op dit onderwerp is hij in 1973 
gepromoveerd bij prof. ir. M.P. Breedveld 
en prof. ir. O.W. Memelink. Hij geeft col- 
lege in de biomedische instrumentatie en 
is projectleider van de werkgroep Biosen- 
soren. 


de beheersing van deze toestanden. In- 
middels begon silicium germanium te 
verdringen als uitgangsmateriaal voor de 
transistorfabricage. Een beslissende stap 
werd gedaan door Atalla (1958), die ont- 
dekte dat een thermisch gegroeide SiO, - 
laag in staat is het grootste deel van de 
oppervlaktetoestanden te neutraliseren. 
Snel werden toen ook de basisprincipes 
van de planaire siliciumtechniek vastge- 
legd, gevolgd door de uitvinding van de 
MOSFET door Kahng en Atalla (1960). 


Werking van MOSFET en CHEMFET 

De MOSFET heeft de volgende opbouw 
(figuur 1). Een p-type siliciumplaatje 
(chip) bevat twee n-gebiedjes (source en 
drain). Tussen source en drain bevindt 
zich een MOS-capaciteit met een dikte 
van het diëlektricum (SiO, ) van 

ca. 0,1 um. Een positieve spanning op het 
metaal (gate) ten opzichte van het sili- 
cium induceert inversie van het geleidings- 
type van het silicium waardoor een elek- 
tronenkanaal ontstaat dat een geleidende 
verbinding tussen source en drain vormt. 
De kanaalgeleiding wordt gegeven door: 


Fig. 1. Doorsnede van de MOSFET. 


source gate 
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waarin u de beweeglijkheid der elektro- 
nen is, Cox de MOS-capaciteit per opper- 
vlakte-eenheid, L de lengte van het kanaal 
tussen source en drain, W de breedte van 
het kanaal, V, de gate-potentiaal en V, de 
drempelspanning is. Uitdrukking (1) geldt 
alleen wanneer V, de waarde V, over- 
schreden heeft en inversie van het gelei- 
dingstype optreedt. 


ge HCx (W/L) (Vs me V‚) 


De drempelspanning V, blijkt chemisch 
beïnvloedbaar te zijn. V‚ kan worden 
geschreven als: 


Vi = Dm -®s — (Oss FO ord Cox E 


en Op/Cox (2) 
Pen ®, zijn resp. de uittreepotentialen 
van gate-metaal en silicium. (Oss + 
Oox)lCox is de bijdrage tot V, veroor- 
zaakt door lading Qs in oppervlaktetoe- 
standen en lading Q,x in het oxyde, die 
eveneens aan het grensvlak Si/SiO, 
gedacht is. De eerste drie termen van (2) 
geven de vlakkebandenspanning Ve: 


Ve= Dn — Ds — (Oss + Oox) Cox (3) 


Ve is de gatespanning waarbij de elektro- 
nenenergiebanden in de halfgeleider tot 
aan het oppervlak vlak zijn.” De laatste 
twee termen van uitdrukking (2) stellen 
de extra gate-spanning voor die nodig is 
om het halfgeleideroppervlak te inverte- 
ren. ®s is de zg. fermipotentiaal bepaald 
door de kristaldotering en Qg is de lading 
per oppervlakte-eenheid in de verarmings- 
laag onder het oppervlak, die eveneens 
van de dotering afhankelijk is. 

Chemische sensoren (CHEMFET's) berus- 
tend op het MOS-systeem, opereren met 
het oog op het te gebruiken meetprincipe 
met een zodanige lading Q,‚ dat een in- 
versielading —Q‚x wordt geïnduceerd, zo- 
dat voor V, = 0 reeds kanaalgeleiding 
aanwezig is. Door chemische beinvloeding 
wordt, onder de conditie V, =0, 
gevarieerd, wat zich uit in een verandering 
van de kanaalgeleiding. Voor een goed 
overzicht van soorten CHEMFET’s wordt 
verwezen naar Croset.” In het volgende 
zullen wij ons in het bijzonder bezighou- 
den met de ISFET (ion-sensitive field 
effect transistor). Voor een overzicht van 
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Dr. N.F. de Rooij (28) studeerde schei- 
kunde aan de RU-Utrecht. Na zijn afstu- 
deren trad hij in ZWO-dienst om aan de 
TH-Twente een promotieonderzoek te 
verrichten naar de fysische en chemische 
verschijnselen, die ten grondslag liggen 
aan de werking van de ionengevoelige 
veldeffecttransistor (ISFET). In 1978 
werd dit onderzoek afgesloten met een 
promotie. De promotoren waren prof. dr. 
J. Th.G. Overbeek en prof. dr. MJ. Spar- 
naay. 





recent ISFET-werk wordt verwezen naar 
Janata en Huber” 


Theorie van de ISFET 

Om tot een ionengevoelige structuur te 
komen wordt het gate-metaal van de 
MOSFET van figuur 1 vervangen door een 
elektrolytoplossing nadat source- en 
draincontacten genoegzaam zijn afgedekt. 
Via een referentie-elektrode kunnen wij 
de elektrolytoplossing een bepaalde span- 
ning geven ten opzichte van het silicium 
(figuur 2). 

In de ESO (elektrolyt-oxyde-silicium)- 
structuur wordt de vlakkebandenspanning 
gegeven door 


KB ABBO teln 
(4) 


waarin E, de referentie-elektrodepoten- 
tiaal voorstelt en °lA°*d het elektrisch 
potentiaalverschil tussen het oxyde- 
oppervlak en de elektrolytoplossing. Ten- 
gevolge van de interactie met water ont- 


Fig. 2. ISFET met elektrolyt en referen- 
tie-elektrode. 
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Fig. 3. Gehydrateerd oxyde-oppervlak 
met silanolgroepen. 





staan silanolgroepen aan het oppervlak 
die gedeeltelijk gedissocieerd zijn volgens 
het evenwicht 


= SiOH + H,O 2 = SiO0 "+ H30* 
(oplossing) (5) 


Schematisch is dit evenwicht weergegeven 
in figuur 3. Dit betekent dat het SiO,- 
oppervlak elektrisch geladen zal zijn, re- 
sulterend in een potentiaalverschil °lA9*® 
ten opzichte van de oplossing. Daarom zal 
de H3 O*-concentratie in de oplossing, of 
m.a.w. de zuurgraad (pH), van invloed 
zijn op de oppervlaktelading en dus ook 
op °lA°*®. In de praktijk blijkt dat 
elAoxP van meer ionsoorten dan alleen 
H*-ionen afhankelijk is, zoals bijvoorbeeld 
van Na'- en K' ionen. 


Vervaardiging van de ISFET 

De voor de vervaardiging toegepaste 
planair silicium techniek heeft zijn voor- 
delen reeds bewezen bij de realisatie van 
geïntegreerde schakelingen (IC's). Type- 
rend voor deze techniek is de massafabri- 
cage van zeer kleine gelijkvormige elektro- 
nische componenten of systemen. Dit 
betekent dat de ISFET eventueel met de 
elektronische signaalverwerking op één 
chip kan worden geïntegreerd. 

De op de TH-Twente vervaardigde 
ISFET’s zien eruit als gegeven in figuur 4. 
Tijdens het thermisch groeien van het 
gateoxyde ontstaat door positieve elek- 
tronenvacatures in het overgangsrooster 
Si/SiO,/SiO, een positieve oxydelading 
Oox + Qss. Met het oog op de reeds 
genoemde eis, dat inversie aanwezig moet 
zijn voor Vg =O is gekozen voor een 
n-kanaal -ISFET. 


De vervaardiging van de ISFET geschiedt 
als volgt. Op een siliciumplak van 3 cm 
doorsnee waarin vele gelijkvormige 
ISFET-structuren zullen ontstaan wordt 
door middel van chemical vapour deposi- 
tion (CVD) een met borium gedoteerde 
SiO, -laag (SiO, (B)) van 0,5 um dikte 
gedeponeerd. Fotolithografisch wordt dit 
oxyde verwijderd daar waar later source 
en drain van de ISFET komen. Figuur 4a 
geeft een bovenaanzicht waarin de 
SiO,(B)-laag is gearceerd terwijl figuur 4b 
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Fig. 4. a: bovenaanzicht ISFET-chip. b, 
c, d: doorsneden over AB. e, f: doorsne- 
den over CD. 


een doorsnede door de chip laat zien 
langs de lijn A-B. 

Hierna wordt fosforgedoteerd oxyde door 
middel van CVD over de hele plak gede- 
poneerd. Een tweede fotolithografie stap 
opent het oxyde daar waar het ionen- 
gevoelige oxyde komt (figuur 4c). Dit 
geschiedt door verhitten van de plak in 
een zuurstofatmosfeer en daarna stikstof- 
bij een temperatuur van 1200 °C. Simul- 
taan treden de volgende effecten op 
(figuur 4d): 

— groei van 0,05 um gateoxyde; 

— vorming van source en drain met lange 
n-type geleidende banen (eigen weerstand 
ca. 100 Q) door fosfordiffusie vanuit het 
oxyde; 

— verhoging van acceptordotering buiten 
de eigenlijke ISFET-structuur door 
boriumdiffusie vanuit het oxyde waar- 
door parasitaire inversielagen worden 
onderdrukt. 


Vervolgens worden de contactgaten in het 
fosforgedoteerde oxyde geëtst (figuur 
4e). Opdampen van aluminium en foto- 
lithografisch etsen verschaft de contact- 
vlakjes voor source, drain en bulk. De 
plak met de complete ISFET’s wordt dan 
gekrast en gebroken in individuele chips. 
ledere chip wordt op een stripvormige 


0,520 um 
©) 


doorsnede over A-B 


0,300um 


doorsnede over C-D 


Al, O3-drager bevestigd. Bovendien wor- 
den de drie contactvlakjes met dunne 
gouddraadjes met de geleidingssporen op 
het Al, 03 -substraat verbonden. 

Tot besluit wordt de ISFET geheel met 
siliconenrubber afgedekt, behalve het 
gategebied (figuur 5). 


Meetresultaten 

In principe berust de ionenconcentratie- 
meting met de ISFET op het meten van 
de kanaalgeleiding bij een gegeven refe- 
rentie-elektrodepotentiaal ten opzichte 
van de source. De gemeten ionenconcen- 
traties kunnen worden weergegeven als 
een equivalente gatespanning Ve eq. Een 
voorbeeld geeft figuur 6 waaruit blijkt, 
dat een thermisch oxyde zonder verdere 
bedekking gevoelig is voor H'-, Nd- en 

K* ionen. 

Om de selectiviteit van een ISFET te 
bevorderen zijn experimenten gedaan met 
een laagje parafine, gemengd met koolstof 
en chloorgevoelig Ruzicka-poeder, aange- 
bracht op het gateoxyde. Inderdaad bleek 
de verkregen sensor selectief te zijn voor 
halogenide-ionen met een perfecte 
Nernst-helling van 59 mV per decade voor 
Cl” ionen (figuur 7). 

Veel onderzoek is momenteel gaande naar 
laagjes of membranen met bijzondere 
selectieve eigenschappen. Behalve dat op 
de beschreven wijze bestendige en kleine 
ionengevoelige sensoren zijn te vervaardi- 
gen moet nog worden vermeld dat ook 
veel aandacht is besteed aan elektronische 
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meetmethoden waarbij storende invloe- 
den als temperatuurgevoeligheid worden 
gecompenseerd 

Dergelijke ISFET-systemen zijn toegepast 
voor snelle zuur-basetitraties. Het blijkt 
dat bij titratiesnelheden van 5 ml/minuut, 
het ISFET-meetsysteem afwijkingen klei- 
ner dan 1% oplevert, terwijl de normale 
glaselektrode door zijn traagheid slechts 
bij titratiesnelheden kleiner dan 0,9 ml/ 
minuut binnen de 1% fout blijft. Bij het 
continu meten van min of meer constante 
pH-waarden vertoont de ISFET enige 
meetdrift. Een drift van 1 mV/uur (> 0,2 
pH-eenheid/uur) is geconstateerd, hoewel 
in de literatuur ISFET’s worden vermeld 
met een drift van slechts 2 mV per 

4 dagen. Diverse technologische proces- 
variaties worden thans bekeken op hun 
merites om de drift te verkleinen. 


Nabeschouwing 

Met de in dit artikel beschreven ontwik- 
kelingen op het gebied van ISFET'’s is 

een nieuwe generatie sensoren geïntrodu- 
ceerd, die evenwel niet alleen staat. Ook 
ten aanzien van het meten van niet-chemi- 
sche grootheden (druk, flow, snelheid, 
versnelling, enzovoort) is een ontwikke- 
ling gaande van microminiaturisatie der 
transducers, gebruikmakende van IC-tech- 
nologie. 

De ontwikkeling van ISFET’s is aanvanke- 
lijk gestart ten behoeve van biomedische 
toepassingen. Er bestaat een grote 
behoefte aan intracorporale detectie van 
chemische grootheden, omdat het meta- 
bole evenwicht van de patiënt van belang 
is. Het zal duidelijk zijn dat hiertoe multi- 
sensoren dienen te worden ontwikkeld, 
waarbij de IC-technologie, mogelijk ge- 
combineerd met dunne- en dikke-film- 


Fig. 5. Van rechts naar links: een con- 
ventionele glaselektrodetip, een micro- 
glaselektrode, en een ISFET-elektrode. 
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Fig. 6. H'-K*-en Na'-gevoeligheid van 
ISFET met onbedekt oxyde. 





technologie, de aangewezen weg tot reali- 
satie lijkt. De beschreven ISFET is daar 
nog slechts een eenvoudig voorbeeld van, 
waarbij echter wel blijkt dat naast de 
voordelen van de microminiaturisatie ook 
geheel nieuwe meetconcepten mogelijk 
worden. De ontwikkeling zal dan ook 
verder gaan in de richting van transducer- 
systemen in plaats van transducers sec, 
met automatische parameterbewaking 

en calibratie.” Het ligt voor de hand 

dat in de toekomst delen van de meet- en 
proceselektronica ook bij de transducer- 
chip worden ingebouwd. Het lijkt niet 
uitgesloten dat mettertijd gekozen zal 
worden voor ‘wegwerp’ sensoren, een 
tendens die in de medische wereld hoe 
langer hoe meer gebruikelijk wordt. Voor 
een eenmalig gebruik van de sensoren is 


ISFET4O Vas =100 mV 





Fig. 7. Cl -gevoeligheid van halogenide- 
gevoelige ISFET. 


uiteraard massaproduktie noodzakelijk, 
hetgeen juist door toepassing van de IC- 
technologie mogelijk wordt, tegen relatief 
lage kostprijs. De laatste tijd valt ook de 
belangstelling van de auto-industrie op 
ten aanzien van het ontwikkelen van 
vastestofsensoren die bij de carburateur 
kunnen worden ingebouwd ten behoeve 
van de ontwikkeling van zuiniger en 
milieuvriendelijker verbrandingsmotoren. 
Er zal echter nog een groot aantal proble- 
men dienen te worden opgelost voordat 
praktisch bruikbare transducersystemen 
gemeengoed zijn. Afhankelijk van de 
gebruiksomstandigheden zijn dit proble- 
men op het gebied van de inkapseling der 
sensoren, de in situ calibratie en een inge- 
bouwd systeem voor een vast referentie- 
punt dat voor elke meting noodzakelijk 
is. 


De ISFET zoals beschreven in dit artikel 
is ontwikkeld aan de TH-Twente als 
vrucht van samenwerking tussen de vak- 
groepen Bio-informatica, Vastestofelek- 
tronica, Vastestoffysica en Chemische _ 
Analyse. Behalve de auteurs droegen aan 
dit artikel bij J. Holleman wat betreft de 
ISFET-technologie en O.W. Memelink 
wat betreft de historie van MOSFET -com- 
ponenten en de eindredactie. Een woord 
van waardering is op zijn plaats aan 

A.J. Verloop, die de afmontage van 
ISFET-chips verzorgt. 
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De actieve radiozon 


A.D. Fokker en J.M.E. Kuijpers. 


De corona van de zon is een magnetoplasma waarin (bundels van) snelle elektronen in- 
stabiele situaties doen ontstaan die leiden tot de opwekking van plasmagolven. Deze 
geven, volgens verschillende mechanismen, aanleiding tot de opwekking van radiostra- 
ling die de corona verlaat. De diverse soorten radio-uitbarstingen of ‘“transients’ zijn ge- 
voelige indicatoren voor de processen die zich in de corona afspelen. Uit kenmerken 
van die radioverschijnselen kan plausibel worden gemaakt wat voor processen dat zijn. 


De titel suggereert al dadelijk dat er ver- 
schillende ‘zonnen’ bestaan. Kenden wij 
aanvankelijk alleen de optische zon, met 
zijn scherpe, exact cirkelvormige rand die 
overeenkomt met het fotosfeer-niveau, in 
de jaren -40 leerde men de radiozon ken- 
nen: veel uitgebreider dan de optische 
zon en met meer vage en niet cirkelsym- 
metrische contouren; vervolgens in de ja- 
ren -60 de sterk inhomogene röntgen- 
(figuur 1) en ver-ultraviolet zon. 


De radiostraling van de zon is overwegend 
afkomstig uit de corona, dat is het ijle 
maar hete en sterk geïoniseerde plasma 
rondom de zonneschijf waarvan het be- 
staan vanouds uit totale zonsvetduisterin- 
gen bekend was. 


Opmerkelijk is de abrupte overgang, over 
een interval van slechts enkele duizenden 
kilometers, van de relatief koele (7 = 
5800 K) en dichte (MN = 10° m*) foto- 
sfeer naar de hete (1 à 2 x 10° K) en ijle 
(NS 10! m°) corona. De hoge tempe- 
ratuur van de corona en de daartoe beno- 
digde verhitting wordt in verband gebracht 
met de aanwezigheid van een circa 200 000 
km dikke zone onder de fotosfeer waarin 
het energietransport hoofdzakelijk door 
convectie plaatsvindt. De zogenaamde 
granulatie, waarmee het zonoppervlak in 
wit licht overdekt is, is van die convectie 
een manifestatie en men moet aannemen 
dat mechanische energie in de vorm van 
(akoestische en magnetische) golven en 
schokken aan de corona wordt toege- 
voerd. Wellicht zijn in het verhittingsme- 
chanisme ook magnetische velden in het 
spel zoals die veelvuldig door het zonop- 
pervlak uitbreken (waar dat het geval is 
ontstaat een zonnevlek of een kleine zgn. 
‘pore ”). De mechanische en de eventuele 
magnetische energieflux vertegenwoordigt 
maar een kleine fractie (10°) van de stra- 
lingsflux van de zon waarin wij ons plegen 
te koesteren. Deze magnetische velden 
vinden hun oorsprong beneden de foto- 
sfeer waar ze door differentiële stromin- 
gen worden versterkt. De gasdruk is hier 
namelijk zo hoog dat de velden praktisch 
geheel door de thermo-en hydrodynamica 
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van het gas worden bepaald. Omdat een 
veldbuis steeds in hydrostatisch evenwicht 
blijft met zijn omgeving wordt de gas- 
dichtheid in de buis na versterking van 
het veld geringer dan in zijn omgeving. De 
buis met het lichtere gas rijst nu omhoog, 
het coronale plasma in. Is in de fotosfeer 
en daar beneden de gasdruk van overheer- 
send belang, in de corona bepalen daaren- 
tegen de magnetische velden de dynamica 
van het gas. De oprijzende magnetische 





Fig. 1. De inhomogene structuur van de 
corona is duidelijk te zien op deze rönt- 
genopname in de golflengtebanden 3-32 
A en 44-54 À, gemaakt op 18 juni 1973 
aan boord van Skylab (American Science 
and Engineering, Cambridge, Mass. ). 





structuren bereiken hier een nieuw even- 
wicht. Een deel van de magnetische ener- 
gie komt daarbij vrij in de vorm van ver- 
hitting en versnelling van deeltjes. Op 
spectaculaire wijze geschiedt dit in een 
zonnevlam (figuur 2). 


De zonnevlam 

Het meest centrale fenomeen van de ge- 
beurtenissen die zich in een activiteits- 
centrum afspelen is de zonnevlam met 
een verscheidenheid van daarmee verbon- 
den effecten (figuur 2), in het bijzonder 
de radio-uitbarstingen (figuur 3). De zon- 
nevlam (“flare’) is een explosief proces 
waarbij binnen enkele minuten energie 
beschikbaar komt voor verhitting en deel- 
tjesversnelling (102 — 10° J). Zonder 
twijfel wordt de energie ontleend aan 
magnetische energie die tijdens de op- 
bouwfase is opgeslagen. De sterkste vlam- 
men treden bij voorkeur op in activiteits- 


centra waarin de magnetische configuratie 


(verdeling van N- en Z-magnetisme in de 
vlekkengroep en loop van de zgn. ‘neu- 
trale’ lijn) gecompliceerd is. 


Uit het gebied van de primaire energie- 
omzetting ontvangt men op aarde voor- 
namelijk röntgenstraling (remstraling) en 
radiostraling in het microgolfgebied (gyro- 
synchrotron-en plasmastraling). De ver- 
hitting en versnelling van deeltjes in een 
zonnevlam leidt verder tot verhitting en 
‘verdamping’ van de dichte fotosfeer 

aan de voetpunten van de magnetische 
veldlijnen. De waarnemingen tonen daar- 
bij lussen die afhankelijk van dichtheid 


Ke 
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Dr. A.D. Fokker jr. (53) studeerde sterren- 
kunde te Leiden. Van 1952 tot 1962 was 
hij werkzaam bij het Staatsbedrijf der PTT 
in een onderzoeksgroep die de invloed 
van de zonsactiviteit op de propagatie- 
condities in de ionosfeer bestudeerde. In 
1960 promoveerde hij op een studie naar 
de eigenschappen van verhoogde radio- 
straling van de zon. Sinds 1962 is hij ver- 
bonden aan de Sterrewacht Sonnenborgh 
te Utrecht, vanaf 1966 als lector. Zijn on- 
derzoek is gericht op de zonsactiviteit, in 
het bijzonder de zonne-radioverschijn- 
selen. 


en temperatuur voornamelijk oplichten 
in Ha, het ultraviolet of zachte röntgen- 
straling (figuur 2). Bij grotere zonnevlam- 
men is steedseen dicht filament aanwezig, 
een reusachtige elektrische stroom- 

draad in de corona **. Bij de uitbarsting 
van de zonnevlam beweegt het filament 
zich met grote snelheid (500 km*s"*) om- 
hoog. Hierbij kunnen zich schokken vor- 
men, wordt de magnetische-veldstructuur 
herordend, worden versnelde deeltjes in 
magnetische lussen opgeslagen en ont- 
snappen energetische deeltjes langs ‘open’ 
veldlijnen naar de interplanetaire ruimte. 


Als gevolg van deze verstoringen spelen 
zich in het magnetoplasma boven een ac- 
tiviteitscentrum allerlei plasmafysische 
processen af die wij gewaarworden aan 
een aantal radioverschijnselen®. Het is de 
kunst om plausibel te maken welk proces 
voor welk verschijnsel verantwoordelijk is 
en om de verschijnselen als diagnostisch 
middel te gebruiken voor de fysica van de 
actieve zon. 


Plasma- en stralingsprocessen 

Alvorens de belangrijkste solaire radiofe- 
nomenen te bespreken noemen wij enige 
processen in het coronale plasma die voor 
de interpretatie daarvan belangrijk zijn. 


1. Gyro-synchrotronstraling. Een elek- 
tron dat zich in een magnetisch veld be- 
weegt, ondervindt een lorentzkracht en 
zendt tengevolge hiervan elektromagneti- 
sche straling uit. Heeft het elektron een ex- 
treem relativistische energie, dan noemt 
men de uitgezonden continuümstraling 
synchrotronstraling. In het subrelativisti- 
sche geval wordt de straling voornamelijk 
uitgezonden bij de (lagere) harmonischen 
van de cyclotronfrequentie en noemt men 
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de straling gyrostraling. In het algemeen is 
in de corona de plasmafrequentie aanzien- 
lijk groter dan de gyrofrequentie. Omdat 
nu de bij een zonnevlam versnelde deeltjes 
subrelativistisch zijn is gyrostraling veelal 
het belangrijkste stralingsmechanisme. 


2. Plasmastraling. Indien de uitgezonden 
elektromagretische straling niet meer kan 
worden opgevat als de somstraling van in- 
dividuele deeltjes, spreekt men van plas- 
mastraling. Het plasma speelt in dat geval 
een essentiële rol in de opwekking van de 
elektromagnetische straling. Veelal ge- 
beurt dit als door het optreden van een 
instabiliteit het niveau van de altijd al wel 
aanwezige (thermische) plasmagolven aan- 
merkelijk versterkt wordt. Een deel van 
die plasma-energie kan door koppeling 
van de plasmagolven onderling of met 
andere golfsoorten in radiostraling wor- 
den omgezet die de corona verlaat. 


3. Anstabiliteit. Indien een (homogeen) 


Fig. 2. Schematische voorstelling van de 

gebeurtenissen tijdens een zonnevlam in 

dwarsdoorsnede. 

UF: uitbarstend filament 

MR: reconnectie van magnetische velden 
en deeltjesversnelling 

u: microgolfstraling 

CC: coronale condensatie (TX 3 xl 0°K ) 

UV: ultraviolette emissie 


Ha: lus met invallend koel materiaal 
(T=2x10° K) 

X: röntgenstraling 

V: ‘verdamping’ van dicht gas. 

e: elektronen 

B: magnetisch veld 





fotosfeer 





Dr.J. Kuijpers (33) studeerde theoretische 
sterrenkunde aan de RUU. In 1971 trad 
hij in dienst bij het Sterrekundig Instituut 


te Utrecht voor theoretisch onderzoek op 


het gebied van de radiostraling van de 
zon. In 1975 promoveerde hij op de ver- 
klaring van structuurlijke continuümuit- 
barstingen in het radiogebied, welke in de 
corona van de zon worden opgewekt na 
een grote zonnevlam. Daarna verbleef hij 
een jaar aan het Institute of Astronomy 
te Cambridge, U.K., om zich te bekwa- 
men in de plasma-fysica van extragalac- 
tische radiobronnen. In 1976 trad hij we- 
derom in dienst bij het Sterrekundig In- 
stituut in Utrecht om onderzoek te ver- 
richten in de extragalactische plasma- 
astrofysica. 


plasma niet in thermisch evenwicht ver- 
keert en voldoende ijl en heet is, kan het 
systeem relaxeren door de versterking van 
plasmagolven (die steeds in de ruis voor- 
komen). De kinetische energie van een be- 
paalde groep hoogenergetische deeltjes 
neemt daarbij af en komt ten goede aan 
elektromagnetische en bewegingsenergie 
van het hele plasma in de vorm van gol- 


schok (II) 


ee 
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ven. De aard van de opgewekte golven 
hangt af van de toestand van het verstoor- 
de plasma. Twee belangrijke instabilitei- 
ten zijn de zogenaamde verlieskegelinsta- 
biliteit en de ‘knobbel-in-de-staart’ insta- 
biliteit. Zo’n knobbel’ geeft aanleiding 
tot plasma-Cerenkov-straling van Lang- 
muir-golven. Dit zijn elektrostatische gol- 
ven bij de plasmafrequentie van de elek- 
tronen. Op de verlieskegel komen wij 

nog terug. 


4. Omzetting van golven. De bewegings- 
vergelijking van de plasmadeeltjes bevat 
de niet-lineaire convectieve term en de 
niet-lineaire lorentzkracht. Dat wil zeggen 
dat de respons van de deeltjesbeweging 
wordt bepaald door zowel de verstoorde 
beweging zelf als door het verstoorde veld. 
Bijgevolg kunnen de velden van meerdere 
plasmagolven samenwerken en nieuwe 
elektrische stromen (elektronenbewegin- 
gen) produceren die zelf weer als bron 
fungeren voor gelijksoortige of andere 
plasmagolven en voor (transversale) elek- 
tromagnetische straling die de (aardse) 
waarnemer kan bereiken. Is het golfniveau 
voldoende hoog, dan combineren de inter- 
agerende golven tot nieuwe golven op de 
som- en de verschilfrequenties. 


Indien een hoogfrequente golf koppelt 
met een andere laagfrequente golf of een 
fluctuatie in de deeltjesverdeling, wordt 
straling gegenereerd bij nagenoeg de fre- 
quentie zelf van die (hoogfrequente) golf. 
Is dat een plasmagolf, dan wordt radio- 
straling als fundamental’ uitgezonden. 
Treedt echter koppeling op tussen twee 
hoogfrequente golven onderling, dan kan 
combinatiestraling optreden bij de eerste 
harmonische frequentie (‘harmonic’). 


De ruisstorm 

De ruisstorm is het meest veelvuldig voor- 
komende fenomeen op metergolflengten 
(A=l à 5 m) en gaat altijd gepaard met het 
optreden van een zonnevlekkengroep. De 
straling is afkomstig van een brongebied 
op een hoogte van 0,2 à 0,6 Re (hoe kor- 
ter de golflengte hoe lager; 1 Ro is één 
zonstraal: 6,96 x 10° km) en met een uit- 
gebreidheid van 2 à 4 boogminuten van 
de aarde af gezien. De ruisstorm bestaat 
uit een relatief langzaam variërend achter- 
grondscontinuüm met daarop gesuperpo- 
neerd kortstondige (0,3 â 0,6 sec.) smal- 
bandige (3 à 6 MHz) stoten (figuur 4). 

De meeste stormen zijn sterk circulair ge- 
polariseerd. 


De stormstoten vertonen helderheidstem- 
peraturen tot 10! K,‚ overeenkomend 
met een deeltjesenergie van 100 MeV. 
Zulke energierijke deeltjes worden slechts 
bij hoge uitzondering, door de krachtigste 
zonnevlammen, geproduceerd. Van de 
radio-uitbarstingen is de ruisstorm wel het 
minst begrepen. Wegens de hoge helder- 
heidstemperatuur lijkt de verklaring te 
moeten worden gezocht in een coherent 
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Fig. 3. Schematische weergave van de 
soorten radio-uitbarstingen tijdens en na 
een zonnevlam. Naast snel driftende uit- 
barstingen (type III) komen ook langza- 
mere driften voor (type II) en continua 
zonder noemenswaardige frequentiedrift 


(type IV en u). 


stralingsmechanisme“ . De stoten zijn 
veelal gerangschikt in fijne structuren die 
een zwakke frequentiedrift vertonen over- 
eenkomend met een karakteristieke mag- 
neto-akoestische snelheid ter plaatse. Bij- 
voorbeeld kan men denken aan magneti- 
sche pulsen die plaatselijk de snelheidsver- 
deling van de elektronen instabiel maken 
en leiden tot de opwekking van plasma- 
golven op frequenties even boven de 
plasmafrequentie fo =(Ne*/e‚m)/ /2m, 
fp (Hz) =8 INN in m°,e is de la- 
ding,7 de massa van een elektron en €o 
de permittiviteit van vacuüm). 

Het continuüm zou dan worden verklaard 


door een algemeen veld van plasmagolven. 


Deze zouden in stand worden gehouden 
door energierijke elektronen die in een 
magnetisch veld worden vastgehouden. 
De ruisstorm lijkt er op te wijzen dat er 
min of meer permanent elektronen wor- 
den versneld, misschien wel op aanmerke- 
lijke hoogten (0,5 à 1 Ro). 


Bundels deeltjes 

Snel driftende (“Type HI) uitbarstingen 
zijn kortstondige (1 à 2 s) stoten, meestal 
in groepen optredend, die in de loop van 
l à 2 seconden een groot frequentiege- 
bied (van de orde van een octaaf of meer) 
in het meter- of decametergebied van 
hoog naar laag doorlopen (figuur 5). Veel 
type III-uitbarstingen gaan gepaard met 
een zonnevlam, maar vooral de decameter- 
type III’s (die eigenlijk een categorie op 
zichzelf vormen) treden onafhankelijk 
van vlammen op. 





‚ Men vat deze uitbarstingen op als het ef- 


fect van een snel langs een magnetische 
veldlijn naar buiten bewegende bundel 
elektronen. Deze bundel genereert Ceren- 
kov-plasmastraling en de zo opgewekte 
Langmuir-golven geven op hun beurt aan- 
leiding tot radio-emissie nabij de plasma- 
frequentie (‘fundamental’) en/of op de 
eerste harmonische daarvan. De meeste 
type [II-uitbarstingen corresponderen ver- 
moedelijk met de harmonische frequentie. 


In de veronderstelling dat de emissie 
plaatsvindt op hetzij de plasmafrequentie 
of het dubbele daarvan, kan men uit een 
dichtheidsmodel van de corona de snelheid 
afleiden waarmee het verstorende agens 
naar buiten beweegt. Men vindt snelhe- 
den van de orde van c/3. (c is de lichtsnel- 
heid.) 


Satellietwaarnemingen hebben deze inter- 
pretatie bevestigd. De evolutie evenwel 
van de plasmagolven die bij de knobbel- 
in-de-staart-instabiliteit worden opgewekt 
is nog steeds niet begrepen. De laatste tijd 
is er van astronomische zijde dan ook een 
sterke belangstelling voor de zogenoemde 
modulatie van Langmuirgolven en solito- 
nen. 


Schokgolven 

Bij grote zonnevlammen ziet men op 
radiogolflengten een zogenoemde type II- 
uitbarsting”. De radiostraling bestaat uit 


Fig. 4. Radiospectrogram van een ruis- 
storm op 8 januari 1971 met type l-uit- 
barstingen die gerangschikt zijn in twee 
ketens met zwakke negatieve frequentie- 
drift. De opname is gemaakt met de 60- 
kanaals radiospectrograaf (thans te 
Dwingeloo). De zwarte kanalen zijn uit- 
gevallen en het regelmatige onderbroken 
signaal is door mensen gemaakt. 


11.30.45 UT 
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twee parallelle banden die in frequentie 
een octaaf uiteenliggen (figuur 6). De 
emissie verschuift met de tijd van hoge 
naar lage frequenties. Uit interferometri- 
sche waarnemingen heeft men afgeleid 
dat het hier gaat om een verstoring die 
zich met een snelheid van meer dan 1000 
km*s”* door de corona naar buiten be- 
weegt. Directe wit-licht-opnamen vanuit 
een satelliet laten zien dat binnen de ra- 
diostructuur de dichtheid verhoogd is. 
Deze soort radio-uitbarstingen is vermoe- 
delijk van schokken afkomstig die zich 
met een snelheid van twee- tot driemaal 
de lokale Alfvén-snelheid voortplanten. 
(De Alfvén-snelheid is de snelheid waar- 
mee geringe magnetische verstoringen 
zich langs het veld voortplanten.) In het 
(botsingsloze) schokfront verandert het 
magnetische veld aanzienlijk over een 
korte afstand. 


Er zijn dus sterke elektrische stromen die 
vermoedelijk aanleiding geven tot de op- 
wekking van elektronenplasmagolven 
(Langmuir-golven). Verstrooiing en onder- 
linge koppeling van deze plasmagolven 
leiden dan tot elektromagnetische straling 
bij de fundamentele en bij de eerste har- 
monische van de lokale plasmafrequentie. 
Dat deze twee mechanismen inderdaad in 
het spel zijn vindt een fraaie bevestiging 
in de harmonische structuur van de type 
[T-uitbarsting. 


Magnetische flessen 


Continuüm (‘Type IV’) uitbarstingen zijn 
algemene verhogingen van het ruisniveau 
die zich over het grootste deel van het 
radiospectrum manifesteren. Vrijwel al- 
tijd is een type IV-uitbarsting het gevolg 
van het optreden van een grote zonne- 
vlam, waarbij plotseling elektronen tot 
(sub )relativistische snelheden worden ver- 


Fig. 5. Radiospectrogram van een groep 
type HI-uitbarstingen op 25 februari 1971 
opgenomen met de radiospectrograaf te 
Culgoora. (LD. Palmer, R.P. Lin, Proc. 
Astron. Soc. Australia 2 (1972) 101) 
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Fig. 6. Radiospectrogram van een type 
H-uitbarsting op 9 october 1969 opge- 
nomen met de radiospectrograaf te Cul- 
goora. (G.A. Dulk, Proc. Astron. Soc. 
Australia 1 (1970) 308) 


sneld. Elektronen die worden vastgehou- 
den in een magnetische lus spiraliseren 
heen en weer tussen magnetische ‘spie- 
gels’ en stralen daarbij (gyro-)synchrotron- 
straling uit. Maar tevens werken zij in 

op het plasma doordat zij dit instabiel 
maken voor plasmagolven. In het bijzon- 
der speelt hierbij een rol de verlieskegel- 
instabiliteit. 


Met verlieskegel wordt bedoeld dat de 
richtingsverdeling van de snelheden der 
snelle elektronen anisotroop is doordat 
snelle elektronen met bewegingsrichting 
binnen een zekere kegel met as evenwijdig 
aan het magnetisch veld ondervertegen- 
woordigd zijn (figuur 7). Zulke elektro- 
nen, met geringe spoedhoek, dringen na- 
melijk door tot dichtere lagen waarin ze 
in korte tijd (ongeveer 10 * s) hun energie 
door botsingen verliezen. Deze anisotrope 
verdeling van snelle elektronen is gesuper- 
poneerd op een dichte isotrope thermi- 
sche achtergrondverdeling die veel kouder 
is (T'= 2x10° K). Deze anisotropie leidt 
tot het ontstaan van plasmagolven op een 
ver-boven-thermisch energieniveau. Een 
kleine fractie van de plasma-energie lekt 
weg in de vorm van radiofrequente stra- 
lingsenergie. 


In de verschillende spectrale gebieden 

(microgolf, decimeter, meter) heeft het 
type IV-continuüm van elkaar verschil- 
lende kenmerken (figuur 3) en de stra- 
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lingsmechanismen zijn, mede vanwege de 
verschillende hoogten van oorsprong van 
de straling, niet noodzakelijk identiek. Op 
bijzonder karakteristieke wijze mani- 
festeert de type IV-uitbarsting zich in het 
decimetergebied door een hoge helder- 
heidstemperatuur (10!* K), een smalle 
bandbreedte (< 100 MHz) en door de 
merkwaardige modulaties die dikwijls in 
het continuüm optreden (figuur 8). Weder- 
om wijzen deze waarnemingen op een col- 
lectieve aard van het emissiemechanisme. 
Op radiospectrografische opnamen met 
goed tijd- en frequentiescheidend vermo- 
gen en goede contrastwerking leveren die 
een intrigerend schouwspel op. In hoofd- 
zaak zijn er drie van zulke modulaties die 
wij nu zullen bespreken. 


1. Zebrapatronen. In het frequentietijd- 
vlak zijn dit in een quasi-periodiek patroon 
gerangschikte vermeerderingen en vermin- 
deringen van de intensiteit. De typische 
afstand tussen aangrenzende ‘ribbels’ is 

3 à 8 MHz. Een zebrasysteem verschuift 
dikwijls in frequentie, waarbij telkens ver- 
anderingen in dat verloop optreden. 


Dit verschijnsel kan worden geduid als 
een verlieskegelinstabiliteit die leidt tot 
de opwekking van plasmagolven op de bo- 
venste hybride frequentie fun =({p° + 
jen op plaatsen waar de verhouding van 
de frequentie fn tot de gyrofrequentie 
fg heeltallig is (fg = (eB)/(2n m), fg (Hz) 
=28x 10!° B(T), B is de magnetische 
fluxdichtheid). Een fractie van de energie 
van deze golven wordt omgezet tn elek- 
tromagnetische straling op nagenoeg de- 
zelfde frequentie. 


Indien de variatie van fn langs de ver- 
schillende ‘resonante’ oppervlakken voor 
de verschillende heeltallige verhoudingen 
funlfg voldoende klein is (Afon S/e), 
dan zullen op verschillende frequenties 
emissiebanden optreden (figuur 8). Aan- 
gezien de lokatie van de resonante opper- 
vlakken erg gevoelig is voor de precieze 
veld- en dichtheidsstructuur in het bron- 
gebied, zal een geringe magnetohydrody- 
namische relaxatie leiden tot een merk- 
bare verschuiving van de ‘resonante’ op- 
pervlakken naar andere frequenties; van- 
daar de ‘drift’ in de zebrapatronen”. 
Indien daarentegen het brongebied niet 
homogeen genoeg is, zullen de bovenste 
hybride frequenties langs de verschillende 
oppervlakken elkaar overlappen en ziet 
men niets anders dan een continuüm. 
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Fig. 7. Schets van het brongebied in 
doorsnee van een type IV-continuüm met 
fijne structuren. De snelheidsverdeling 
van de deeltjes is aangegeven op twee 
plaatsen in de magnetische lus. De ver- 
schillende golfsoorten zijn aangeduid met: 
uh: bovenste hybride. 
w.; whistler 
t: ontsnappende elektromagnetische 
straling. 
Verder zijn twee opeenvolgende ‘resonan- 
te’ oppervlakken aangegeven, waarop fun/ 
fg heeltallig is (gelijk aan N en respectie- 
velijk N+1 ). 


2. Fiber uitbarstingen. Deze modulatie 
ziet eruit als een geïsoleerde ribbel in het 
tijd-frequentievlak met een ‘emissie’kant 
en een aangrenzende “absorptie’kant aan 
laagfrequente zijde. Daarbij ziet het er- 
naar uit dat een hoeveelheid energie is 
overgeheveld van de ene frequentie naar 
een naburige (Af = 3 MHz) hogere 
frequentie. De fiber vertoont een drift- 
snelheid die inligt tussen die van type III- 
en van type IT-uitbarstingen. 


Dit verschijnsel kan worden geduid als 
het tweede effect van de verlieskegelin- 
stabiliteit waarbij zgn. “‘whistlers’ worden 
opgewekt. 


Een whistler is een elektromagnetische 
golf op een frequentie f‚ lager dan de gy- 
rofrequentie, die een voorkeursrichting 
van voortplanting langs het magnetische 
veld heeft. Whistlers zijn welbekend van 
de karakteristieke audiofrequente fluit- 
toontjes die worden geproduceerd door 
bliksems en die langs aardmagnetische 
veldlijnen heen en weer reizen. Whistlers 
kunnen niet uit de corona ontsnappen, 
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en hun aanwezigheid in de corona kan 
dan ook niet direct worden aangetoond. 


In de corona zouden whistlers bij voorkeur 
worden opgewekt nabij de ‘voetpunten’ 
van magnetische lussen, waar de anisotro- 
pie van de elektronenrichtingsverdeling 
maximaal is. Als golfpakketjes reizen de 
whistlers langs het veld omhoog. Waar ze 
de veel hoger frequente (f = fyn) plasma- 
golven ontmoeten, zullen ze daarmee 
koppelen en elektromagnetische straling 


ema Pr dee 
bss eat els 


ets drs: 








produceren op de frequentie fun + fw. 


(N.B. Bij f‚‚ —f kan elektromagnetische 
straling zich niet voortplanten in het plas- 
ma.) Plasmagolven die combineren met 
een whistler zijn daarmee verloren voor 
het opwekken van elektromagnetische 
straling (door strooiing aan fluctuaties) 
op hun eigen fyn — frequentie; vandaar 
wellicht de met de emissie gepaard gaan- 
de intensiteitsreductie (figuur 8). 


Het frequentieverloop van de fiber zou 
een afspiegeling zijn van de naar buiten 
gerichte verplaatsing van het whistler-pak- 
ket (met een snelheid van 21 — 28 maal 
de Alfvén-snelheid), terwijl de waargeno- 
men richelafstand met de whistlerfrequen- 
tie overeenkomt (0,1 — 0,4 f,). Uit beide 
gegevens kan men onafhankelijk de plaat- 
selijke veldsterkte in de bron bepalen. In- 
derdaad blijken deze waarden met elkaar 
in overeenstemming”. 


3. Plotselinge reducties, versnelling in 
pulsen. De intensiteit van de type IV-con- 
tinuümstraling vertoont af en toe een serie 
reducties die vrijwel simultaan plaatsvin- 
den over een frequentiebereik van de orde 
van een half octaaf à een octaaf. Even plot- 
seling als de vermindering van intensiteit 
geschiedt, komt na ca. 0,3 s de intensiteit 


Fig. 8. Voorbeelden van fijne structuren 
in type IV-continua opgenomen met de 
radiospectrograaf te Dwingeloo: 
bovenaan: zebrapatronen; 

midden: fibers; 

en onder: plotselinge reducties. 
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weer terug op het uitgangsniveau. Binnen 
een serie zijn de afzonderlijke reducties 
meestal onregelmatig in de tijd verdeeld, 
maar soms herhalen ze zich quasiperio- 
diek. 


Ter verklaring denkt men aan de mogelijk- 
heid dat de verlieskegelinstabiliteit wordt 
opgeheven door de plotselinge injectie van 
snelle deeltjes (elektronen) die, met een 
bewegingsrichting ruwweg langs het mag- 
netisch veld, de verlieskegel tijdelijk op- 
heffen’. Aangezien een isotrope maar 
verder willekeurige deeltjesverdeling sta- 
biel is, zodat geen plasmagolven worden 
opgewekt, zullen de bestaande instabili- 
teiten onmiddellijk worden onderdrukt 
en de reeds aanwezige golven worden ge- 
dempt. 


Met deze interpretatie verplaatst zich het 
probleem naar de mogelijkheid van een 
mechanisme waarbij versnelling pulsge- 
wijs optreedt. Er zijn twee mogelijkheden 
geopperd. 

Volgens het ene mechanisme zouden ion- 
akoestische golven, opgewekt door elek- 
trische stromen in het vlamgebied, aanlei- 
ding geven tot Langmuir-golven?® . Ener- 
gie van die golven komt ten goede aan 
elektronen die door Landau-demping 
worden versneld. Deze energie-uitwisse- 
ling treedt op tussen golven en deeltjes 
waarvan de snelheid in de richting van de 


golfvector gelijk is aan de fasesnelheid 
van de betreffende golf. Een en ander 
wordt beschreven door een stelsel van ge- 
koppelde niet-lineaire vergelijkingen (ge- 
noemd naar Lotka en Volterra) die slechts 
periodieke oplossingen toelaten. Oor- 
spronkelijk waren die vergelijkingen op- 
gesteld ter verklaring van periodieke fluc- 
tuaties in populaties van roof- en prooi- 
dieren. 


Volgens het andere mechanisme zou ver- 
snelling in zonnevlammen tot stand ko- 
men door middel van geïnduceerde elek- 
trische velden, analoog aan de in het labo- 
ratorium geproduceerde pulsen van ‘run- 
away’ deeltjes”. 


Een aanwijzing voor het optreden van 
pulsgewijze versnelling wordt ook geleverd 
door de zeer kortstondige (< 0,02 s) stra- 
lingspieken die zijn gevonden in registra- 
ties van microgolfuitbarstingen met hoog 
tijdscheidend vermogen !®. 


Besluit 

De diagnostiek van de radiovlam zelf, in 
het microgolfgebied, staat in feite nog in 
de kinderschoenen. Voor het komende 
Solar Maximum Mission-jaar staat een 
unieke combinatie van satelliet- en grond- 
waarnemingen op stapel om de zonne- 
vlam te bestuderen. De Nederlandse ster- 
renkundige gemeenschap zal daaraan een 


tweeledige bijdrage leveren: enerzijds zal 
de radiostraling van vlammen worden be- 
studeerd met de Westerbork Synthese Ra- 
dio Telescoop (voor positie-informatie) 
en de radiospectrograaf te Dwingeloo 
(voor fijne structuren in frequentie en 
tijd). Anderzijds zal de röntgenstraling ge- 
lijktijdig worden waargenomen met de 
Hard X-ray Imaging Spectrometer in de 
SMM -satelliet. 
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Boekrecensies 


Thomas S. Kuhn, 

Black-Body Theory and the Quantum 
Discontinuity 1894-1912, 

Clarendon Press: Exford University Press, 
1978. 

XIII, 356 blz., £ 12.50, gebonden. 


In het algemeen laten fysische leerboeken 
de quantummechanica beginnen met 
Planck’s lezing op 14 december 1900 over 
de theorie van de energieverdelingsfunctie 
bij ‘zwarte stralers’. Ook vooraanstaande 
wetenschapshistorici als Kangro, Klein, 
Hermann en Jammer nemen dit stand- 
punt in. Zij interpreteren de introductie 
van de hypothese e = hv als een quantisa- 
tie van de energie E van een oscillator: 
E=0,hv, 2hv,... 

Kuhn — onder wetenschapsfilosofen en 
wetenschapshistorici geen onbekende 

— verdedigt in zijn boek de uitdagende 
stelling, dat deze visie berust op het on- 
juist lezen van Planck’s artikel. Planck wil 
namelijk Boltzmann’s waarschijnlijkheids- 
definitie van de entropie toepassen op het 
probleem van de zwarte straling. Hij telt 
daartoe het aantal oscillatoren met een 
energie ergens tussen nhv en (n + 1)hv, 
waarbij Planck volgens Kuhn uitdrukke- 
lijk niet bedoelt — maar wel zegt — dat de 
energie uitsluitend de randwaarden van 
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het interval zou kunnen aannemen. In 
essentie komt Kuhn’s stelling erop neer, 
dat Planck een zuiver klassieke, continue 
energieverandering voor ogen stond, ter- 
wijl iedereen tot riu toe gedacht heeft, 

dat Planck de energie van een oscillator 
slechts discrete waarden toekende en dat 
hij daarmee de eerste quan tum-voorwaarde 
formuleerde. | 


Zo’n standpunt schept verplichtingen. 
Het strekt tot lof van Kuhn, dat hij op 
nauwgezette wijze zijn stelling tracht te 
bewijzen. In het eerste deel van zijn boek 
gaat hij in op de weg waarlangs Planck 
kwam tot zijn energieverdelingsfunctie. 
Daarbij bespreekt hij onder meer de idee- 
en van Boltzmann en Planck over irrever- 
sibiliteit, moleculaire wanorde en het ge- 
bruik van combinatoriek voor de definitie 
van de entropie. Op zeer gedetailleerde 
wijze behandelt Kuhn de verschillende 
fases in Planck’s theorie. Hij tracht aan te 
tonen dat zowel Planck als zijn lezers in 
de periode 1900-1906 zijn theorie niet als 
een breuk met de klassieke fysica ervoe- 
ren. In het tweede deel van het boek be- 
schrijft Kuhn de groei van het besef dat 
Planck’s theorie een essentiële afwijking 
van gevestigde denkbeelden impliceert. 
Twee artikelen van Einstein en Ehrenfest 
uit 1906 vormden daartoe de eerste aan- 
zetten. Pas nadat Lorentz in 1908 het 
probleem zeer scherp had gesteld, begon 


een wat grotere groep fysici zich van de 
moeilijkheden bewust te worden. Tegen 
1910 was een aantal prominente fysici, … 
waaronder Planck zelf, ervan overtuigd 
dat op een of andere wijze een disconti- 
nuiteit ingevoerd diende te worden. In 
een afsluitend derde gedeelte beschrijft 
Kuhn Planck’s zogenaamde “tweede theo- 
rie’ van 1911-12. Deze theorie zag men 
tot nu toe als het krampachtig vasthou- 
den van een koppige oude man aan de 
klassieke fysica. Kuhn draait ook deze 
waardering 180° om: naar zijn mening is 
dit juist Planck’s eerste (gepubliceerde) 
erkenning van het niet-klassieke karakter 
van zijn theorie. 


Bij de beoordeling van het boek lijkt het 
volgende me van belang. Kuhn doet alsof 
er slechts twee, zeer scherp te bepalen, 
posities mogelijk zijn: Planck voerde op 
14 december 1900 discontinue energie op 
atomaire schaal in, of hij deed dat niet. 
Het eerste standpunt werd tot nu toe in- 
genomen, Kuhn verdedigt het tweede. 
Deze polariserende manier van voorstel- 
len bevordert de duidelijkheid van zijn 
eigen stelling en maakt een vruchtbaar 
debat mogelijk. Het moeilijk te interpre- 
teren historisch materiaal geeft evenwel 
misschien een mogelijkheid de twee posi- 
ties een zodanige ‘lijnbreedte’ te geven 
dat ze elkaar overlappen. | 

C. Hakfoort 
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Het Einstein-observatorium voor 
hoge-energieastrofysica 


J. Heise 


Op 13 november 1978 werd in de Verenigde Staten de tweede satelliet in de serie High 
Energy Astrophysical Observatories, HEAO-B, gelanceerd, die ter gelegenheid van het 
honderdste geboortejaar van Albert Einstein, het Einstein-observatorium werd gedoopt. 
De röntgentelescoop aan boord heeft een ruimtelijk scheidend vermogen en een ge- 
voeligheid, die een verbetering betekent van een factor 100, respectievelijk 1000 ten 
opzichte van het tot dan toe gebruikelijke. Hiermee is HEAO-B het eerste volwassen 
röntgen-observatorium, waarvan de resultaten een kwaliteit bezitten vergelijkbaar met 
die van de huidige grootste optische en radiotelescopen. 

De Nederlandse betrokkenheid is veelzijdig. Het Laboratorium voor Ruimteonderzoek 
te Utrecht heeft een Objective Grating Spectrometer vervaardigd, bestaande uit twee 
soorten transmissietralies van resp. 500 en 1000 lijnen/mm. De Universiteit van Leiden 
heeft een samenwerkingsprogramma voor het waarnemen van dezelfde velden aan de 
hemel met de synthese-radiotelescoop in Westerbork, en verscheidene individuele astro- 
fysici hebben gastwaarnemingsprogramma’s voor waarnemingen met het Einstein-obser- 
vatorium ingediend. De eerste resultaten zijn beschreven in een speciaal aan het Einstein- 
observatorium gewijd nummer van Astrophysical Journal Letters (november 1979). 

In dit artikel beoog ik een korte samenvatting te geven van de status van de röntgen- 
astronomie tot nu toe om daarna aan de hand van de illustraties een indruk te geven 
van de potentiële mogelijkheden van dit röntgenobservatorium. 


Röntgenastronomie 

Onder het röntgenvenster in het elektro- 
magnetische spectrum wordt meestal het 
fotonenergiegebied van 0,1 … 100 keV ver- 
staan. Hierin onderscheidt men de m.b.v. 
scherende inval te focusseren zachte rönt- 
genstraling van 0,1.…3 keV en de alleen 
met mechanische middelen te collimeren 
harde röntgenstraling van 3…..100 keV. 
De kosmos blijkt een grote verscheiden- 
heid van röntgenbronnen te kennen. Vóór 
de lancering van het Einstein-observato- 
rium waren er zo’n 300 bronnen bekend 
plus een ruimtelijk onopgeloste ‘diffuse’ 
component in de röntgenemissie. De rönt- 
genstraling van bronnen van galactische 
oorsprong blijken, voorzover geïdentifi- 
ceerd, veelal samen te hangen met de laat- 
ste stadia van sterevolutie. Jonge super- 
novarestanten emitteren over een uitge- 
strekt gebied van de orde van een lichtjaar 
synchrotronstraling afkomstig van elek- 
tronen van ongeveer 10!* eV in een mag- 
neetveld van ongeveer 10“ Gauss, die bei- 
de gegenereerd worden door een snel ro- 
terende gemagnetiseerde neutronenster, 
het geïmplodeerde restant van de super- 
nova-explosie. De roterende neutronen- 
ster zelf is als radio-röntgen- en gamma- 
stralingspulsar waarneembaar. De uitein- 
delijke energiebron is de rotatie-energie 
van de neutronenster (voorbeelden: de 
Krab-pulsar en de Vela-pulsar). 


Andere röntgenbronnen blijken samen te 
hangen met neutronensterren in nauwe 
dubbelstersystemen. Hier is de energie- 
6, 


bron van gravitationele oorsprong: mate- 
rie valt van de begeleidende ster (bijvoor- 
beeld door een sterwind) in het gravitatie- 
veld van de neutronenster, vormt ten ge- 
volge van het specifieke impulsmoment 
een dunne schijf in Kepler-banen rond de 
neutronenster. Door de differentiële rota- 
tie en viscositeit in deze schijf diffundeert 
de materie naar binnen tot een afstand 
waarop magnetische krachten van het 
magneetveld van de neutronenster domi- 
neren. Op kleinere afstand wordt de ma- 
terie langs de veldlijnen geleid tot de mag- 
netische polen. Dit proces wordt aange- 
duid als accretie, de schijf heet accretie- 
schijf. Bij accretie op een neutronenster 
komt er in de gravitatiepotentiaal een hoe- 
veelheid energie vrij overeenkomend met 
ongeveer 10% van de rustmassa-energie. 
Het proces is zelfregulerend in die zin, dat 
de stralingsdruk bij te grote accretie ver- 
dere accretie verhindert. De maximaal 
waargenomen absolute flux van 10° 
erg/s reflecteert deze grens. Er is nu 

een twintigtal röntgendubbelsterren be- 
kend, waaronder Hercules X-1. In deze 
laatste bron is de grootte van het magne- 
tisch veld aan het oppervlak van de neu- 
tronenster aan de hand van het optreden 
van cyclotronabsorptielijnen bij 60 keV 
en zijn eerste harmonische bij 120 keV 
bepaald op ongeveer 10'* Gauss. 


Een derde klasse bronnen zijn de z.g. ga- 
lactische builbronnen in en in de nabijheid 
van het galactisch centrum. Zij behoren 
zowel relatief als absoluut tot de helder- 
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ste röntgenbronnen, maar vertonen geen 
enkele aanduiding voor het lidmaatschap 
van een dubbelstersysteem. De door de 
Astronomische Nederlandse Satelliet ANS 
ontdekte z.g. röntgenbursters, uitbarstin- 
gen die binnen één seconde plaatsvinden 
en tot minuten kunnen duren, behoren tot 
deze klasse. De precieze aard van deze 
bronnen is nog niet bekend, hoewel er 
sterke aanwijzingen zijn dat ook de rönt- 
genbursters neutronensterren zijn. 
Extragalactische röntgenemissie is gede- 
tecteerd in elk type sterrenstelsel, maar is 
zeer groot voor actieve stelsels zoals Sey- 
fert-stelsels, radiostelsels en quasars. De 
precieze aard van de emissie is niet geheel 
bekend, maar er is een toenemend aantal 
aanwijzingen dat de emissie voor een be- 
langrijk deel is geconcentreerd in een 
uiterst klein gebied nabij het centrum van 
deze stelsels. Daarnaast zijn clusters van 
melkwegstelsels ook, maar nu ruimtelijk 
uitgebreide, röntgenbronnen. Mede door 
het detecteren van Fe-ijnen nabij 6 keV, 
leidt men af dat de emissie hier remstra- 
ling is door een ijl, heet (10° K) inter- 
clustergas. 

De oorsprong van de omnidirectionele- 
diffusieröntgenstraling was tot voor kort 
nog duister. Geïnterpreteerd als waarlijk 
diffuse straling afkomstig van een heet 
gas (5x105 K) tussen de clusters van de 
melkwegstelsels heeft het belangrijke 
kosmologische consequenties. Een be- 
langrijke parameter in kosmologische mo- 
dellen is namelijk de materiedichtheid in 
het heelal. Bij een gemiddelde dichtheid 
kleiner dan een kritieke waarde heeft het 
heelal een oneindig volume (‘het heelal 

is open’) en expandeert eeuwig. Bij een 
dichtheid groter dan de kritieke waarde 
is het volume eindig (‘het heelal is geslo- 
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ten’) en de huidige expansie verandert 
t.z.t. in een contractie. De dichtheid van 
de zichtbare materie is slechts ongeveer 
20% van de kritieke dichtheid. 

Eventueel overal aanwezige ijle koele ma- 
terie zou door absorptie gemakkelijk ge- 
detecteerd kunnen worden en blijkt niet 
aanwezig. De overblijvende mogelijkheid 
is dus hete materie. Een dichtheid van on- 
geveer l atoom per kubieke meter is vol- 
doende om het heelal te ‘sluiten’. 

Echter het verrassend grote aantal bron- 
nen waargenomen met het Einstein-obser- 
vatorium doet vermoeden dat een belang- 
rijke bijdrage wordt geleverd door ruim- 
telijk onopgeloste puntbronnen. 

Het vraagstuk is dus weer open, het heel- 
al waarschijnlijk ook. 


Het Einstein-röntgenobservatorium 

Het 6 meter lange in doorsnede metende 
Einstein-observatorium maakt gebruik 
van de volledige reflectie van zachte rönt- 
genstraling (4 … 120 Á) bij scherende in- 
val onder een hoek van ongeveer 1’. Een 
combinatie van achtereenvolgens hyper- 
bolisch en parabolisch gevormde spiegels 
geeft een afbeelding van de röntgenhemel 
met een gezichtsveld van 75 x 75 boog- 
minuten in het brandvlak. In dit brand- 
vlak kunnen verschillende detectoren 
worden geplaatst. De belangrijkste twee 
zijn een plaatsgevoelige proportionele tel- 
ler (de IPC, Imaging Proportional Counter) 
met een resulterend ruimtelijk scheidend 
vermogen van 1 boogminuut en een chan- 
nel array plate (HRI, High Resolution 
Images), waarbij het ruimtelijk scheidend 
vermogen door de optische eigenschappen 
van de spiegel is beperkt tot 2 boogsecon- 
den. In de stralenbundel kan een door het 
Laboratorium voor Ruimteonderzoek te 
Utrecht vervaardigde transmissietralie 
worden geplaatst, waardoor een spectrum 
verkregen kan worden met een spectraal 
scheidend vermogen van 0,5 à 1 Ä tussen 
Sen 120 Á. Zie voor een volledige be- 
schrijving van het experiment Giacconni, 
referentie 2. 


Supernova-restanten 

Een dankbaar object van studie voor het 
Einstein-observatorium is de röntgenemis- 
sie van supernova-restanten, vanwege hun 
ruimtelijke uitgebreidheid (van de orde 
van een boogminuut). De 1000 jaar oude 
Krab-nevel (figuur 1) ziet er in röntgen- 
straling ongeveer zo uit als in gepolariseerd 
licht. Duidelijk blijkt de centrale pulsar 
die pulseert met een periode van 33 milli- 
seconde en de verhoogde emissie ten 
noorden van de pulsar, die overeenkomt 
met een bekend nevelig actief gebied uit 
optische opnamen. Ook in de 10 000 jaar 
oude Vela-supernova-restant is de centrale 
pulsar, die met een periode van 89 ms al 
bekend was uit radio-, optische en gamma- 
stralingsmetingen, nu ook in röntgenemis- 
sie waargenomen. Wonderlijk genoeg 
blijkt in enkele jongere supernova-restan- 
ten (enkele honderden jaren), zoals Cas- 
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Fig. 1. Opname van de Krab-nevel met 
het Einstein-observatorium. 


siopeia A (figuur 2 uit Murray“) en in Ty- 
cho’s supernova-restant niets van een cen- 
traal object. Deze afwezigheid impliceert 
een bovengrens van 1,5x10° K van een 
mogelijke stellair restant, een factor 5 la- 
ger dan de verwachte temperatuur van 
zo’n recent gevormde neutronenster. Als 
niet in (bijna) alle supernova-explosies 


Fig. 2. Het Einstein-observatorium: de 
supernovarestant Cas A in röntgenemissie. 
Uit: Murray etal. (1979). Courtesy Astro- 
physical Journal Letters. 


neutronensterren worden gevormd dan is 
dit een extra moeilijkheid in het rijmen 
van de grote neutronensterdichtheid (af- 
geleid uit het grote aantal radiopulsars) en 
de toch al relatief lage frequentie van 
waargenomen supernova-explosies (2 tot 
4 per eeuw) in ons melkwegstelsel. 


Een andere mogelijkheid is dat er wel een 
centrale neutronenster aanwezig is, maar 
dat de oppervlaktetemperatuur veel lager 
is dan verwacht. Daar de koeling van een 
neutronenster de eerste 1000 jaar ná de 
(verondersteld hete) vorming volledig 
wordt bepaald door neutrinoverliezen uit 
de kern, kan dit betekenen dat andere tot 
nu toe niet beschouwde processen een rol 
spelen. Dit zou bijvoorbeeld het geval zijn 
als de gevormde neutronenster een kern 
bezit, die bestaat uit pionen. 





Stellaire bronnen 

Een van de verrassendste eerste resultaten 
van het Einstein-observatorium is het grote 
aantal röntgenbronnen dat gevonden is. 
Al uit waarnemingen met de Astronomi- 
sche Nederlandse Satelliet ANS bleek dat 
er behalve ‘klassieke’ galactische bronnen, 
die samenhangen met neutronensterren, 
ook andere typen röntgenbronnen be- 
staan. Deze stellaire bronnen, zoals Sirius, 
Capella, vlamsterren e.d. hebben een ab- 
solute helderheid, die 10° tot 10° keer 
kleiner zijn dan de röntgendubbelster- 
systemen. De door ANS en door sondeer- 
raketwaarnemingen gedetecteerde stellaire 
bronnen blijken nu de uiterste top van de 
ijsberg. Praktisch iedere ster emitteert een 
zekere fractie van zijn optische lichtkracht 
in röntgenstraling. Deze fractie varieert 
van ongeveer 10° voor hete sterren (op- 
pervlaktetemperatuur > 10 000 K) via 
een minimum van 10° voor sterren iets 
heter dan de zon, naar een fractie van on- 
geveer 0,1 voor de koelere sterren (N.B. 
vergelijk dit met een factor 10* voor de 
verhouding röntgenemissie tot optische 
emissie in een bron als Scorpius X-1). De 
aard van de stellaire röntgenemissie is 
waarschijnlijk remstraling uit een optisch 
dun plasma met een temperatuur van 10° 
tot 107 K, dat zich bevindt in een ijle 
corona rond de ster. Van koele sterren, 
zoals de zon, is het bekend dat een coro- 
na kan ontstaan ten gevolge van het feit 
dat de buitenste lagen in de ster in con- 
vectief evenwicht zijn. 

In hete sterren is de theorie nog niet dui- 
delijk, maar hangt het bestaan van een 
corona waarschijnlijk mede samen met de 
grote stralingsdruk en de daardoor lage 
netto versnelling aan het oppervlak van de 
ster. 


Andromeda-nevel 

De bestudering van röntgenbronnen in 
onze naaste buur, het spiraalstelsel M31 
(de Andromeda-nevel) biedt zekere voor- 
delen ten opzichte van de bronnen in ons 
eigen melkwegstelsel. De afstanden zijn 
bekend en daarom is bijvoorbeeld de ab- 
solute helderheidsverdeling nauwkeuriger 
te bepalen. De spiraalstructuur en verdere 
morfologie is duidelijker zichtbaar dan in 
ons stelsel zodat bijvoorbeeld bronposities 
eenduidiger zijn te interpreteren. 

Vóór de waarnemingen van het Einstein- 
observatorium was bekend dat de Andro- 
meda-nevel een zwakke, nog juist te detec- 
teren röntgenbron was. Figuur 3 toont de 
resultaten (Van Speybroeck” ) met de 
positiegevoelige proportionele teller (ruim- 
telijke resolutie 1’). De nevel bestaat uit 
een groot aantal puntbronnen. Het diffuse 
centrale, in dit plaatje ruimtelijke onopge- 
loste deel blijkt ook uit een groot aantal 
(45) puntbronnen te bestaan, zoals de 
opname in figuur 4 met de High Resolu- 
tion Imager (resolutie 2”) aantoont. In 
totaal zijn 69 röntgenbronnen te onder- 
scheiden. Na vergelijking met een optische 
opname is de spiraalstructuur in de posi- 
34 





Fig. 3. Het Einstein-observatorium: de 
spiraalnevel M31 in röntgenemissie (reso- 
lutie 1). Uit: Van Speybroeck et al. 
(1979). Courtesy Astrophysical Journal 
Letters. 


‘ma 


ties van röntgenbronnen duidelijk waar- 
neembaar. In totaal is zo’n 60% van de ge- 
identificeerde bronnen geassocieerd met 
z.g. populatie I-objecten (jonge sterasso- 
ciaties, gebieden met interstellaire stof, 
gebieden van geïoniseerd waterstof enz). 
De resterende 40% behoort tot de galac- 
tische builbronnen (galactic ‘bulge’ sour- 
ces) nabij het centrum van M31. Hiervan 


Fig. 4. Het Einstein-observatorium: het 
centrum van de spiraalnevel M31 in rönt- 
genemissie met een resolutie van 2”. Uit: 
Van Speybroeck et al. (1979). Courtesy 
Astrophysical Journal Letters. 


is zo’n 10% geïdentificeerd met bolvor- 
mige sterhopen. Naast de overeenkomsten, 
zijn natuurlijk ook de verschillen met ons 
eigen stelsel interessant. Opvallend is de 
bron, die binnen de nauwkeurigheid van 
2" in de kern van M31 wordt waargeno- 
men. Met een helderheid van 10° erg/s 
is deze als röntgenbron niet ongewoon 
sterk, maar is wel een factor 1000 sterker 
dan de bronnen binnen 0°.5 van het cen- 
trum van ons eigen stelsel. De radio-emis- 
sie van de kern van M31 is slechts 6% (op 
1415 MHz) van die van Sgr A, de radio- 
kern van de melkweg. Dit suggereert dat 
de röntgenemissie in M31 van thermische 
aard is, mogelijk accretie van gaswolken 
op een zwaar zwart gat (10° zonmassa’s). 
De constantheid van de emissie is een ver- 
dere ondersteuning van dit model. 


Clusters van melkwegstelsels 
Röntgenwaarnemingen zijn van vitaal be- 
lang voor het begrijpen van de structuur 
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en evolutie van clusters van melkwegstel- 
sels en is nu eerst goed mogelijk met het 
Einstein-observatorium. Een zeer groot 
deel van de waarnemingsprogramma’s is 
dan ook gericht op extragalactische rönt- 
genastronomie. Het lijkt erop dat een 
cluster begint als een grote wolk van stel- 
sels, die geleidelijk instort onder invloed 
van zijn eigen gravitatieveld om tenslotte 
een dynamisch evenwicht te bereiken 
wanneer een dichte kern en een uitgestrek- 
te halo van melkwegstelsels is gevormd. 
Men kan de snelheidsdispersie van de 
melkwegstelsels in een cluster gebruiken 
als maat voor de gravitatiepotentiaal. De 
Virgo-cluster heeft een lage dispersie en is 
een voorbeeld van een cluster in een vroeg 
evolutiestadium. Figuur 5 toont enkele 
röntgen-isofoten (Forman! ) gesuperpo- 
neerd op een optische opname van een 
deel van de groep. Bijna ieder stelsel blijkt 
een ruimtelijk uitgebreide diffuse bron 
van röntgenstraling te zijn veroorzaakt 
door remstraling van een ijl (4x10 ° deel- 
tjes per cm” ) heet gas met een tempera- 
tuur van 10’ K. Blijkbaar is dit gas gebon- 
den door de gravitatiepotentiaal van het 
stelsel. Dit geheel is gesuperponeerd op 
een diffuse harde röntgenbron (met een 
uitgebreidheid overeenkomend met de ge- 
hele Virgo-groep van stelsels, veel groter 
dan de opname in figuur 3). De harde 
röntgenmissie is afkomstig van remstraling 
van een nog ijler (5x10”* deeltjes cm ”°) 
gas met een temperatuur van 10° K dat 
wordt gebonden door de gravitatiepoten- 
tiaal van de Virgo-cluster als geheel. For- 
man et al!) interpreteren de lokale 10” 
K component als gas dat wordt ‘opge- 
veegd’ door de eigen beweging van de stel- 
sels ten opzichte van het hete 10° K-intra- 
clustergas. Een vergelijkende studie 
(Jones®) van de röntgenemissie van ver- 
schillende rijke, nabije clusters van melk- 
wegstelsels bevestigt dit beeld. De clusters 
waarin de röntgenemissie is opgebroken in 
klonten rond individuele stelsels vertonen 
een lagere temperatuur en lage snelheids- 
dispersie, terwijl de clusters met hoge snel- 
heidsdispersie een gladde, centraal gecon- 
denseerde röntgenintensiteitsverdeling 
vertonen waarin de temperatuur veel ho- 
ger is. In dit laatste geval wordt de verde- 
ling van gas meer bepaald door de gravita- 
tiepotentiaal van de cluster als geheel dan 
door die van de individuele stelsels. Deze 
clusters hebben ook een kleinere fractie 
aan spiraalstelsels, daar het interstellaire 
gas en stof in deze stelsels is opgeveegd 
door het dichtere en hetere intracluster- 
gas. 


De transmissietralies 

Het gebruik van transmissietralies is een 
nieuwe en veelbelovende techniek in de 
röntgenspectroscopie ontwikkeld in het 
Laboratorium voor Ruimteonderzoek te 
Utrecht. De spectrale resolutie van pro- 
portionele tellers voor zachtere röntgen- 
straling met een golflengte groter dan bij- 
voorbeeld 20 Á is buitengewoon slecht 
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Fig. 5. Het Einstein-observatorium: het 
centrale deel van de Virgo-cluster met 
daarop gesuperponeerd isofoten van het 
röntgenbeeld. Uit: Forman et al. (1979). 
Courtesy Astrophysical Journal Letters. 


(AN/A — 1), terwijl voor de tralies 

AN/X gelijk is aan ongeveer 0,5 Á/\.Ten 
tijde van het prepareren van de tralies 
voor het Einsteinobservatorium was de 
technologie echter nog niet zodanig dat 
tralies van voldoende dikte gemaakt kon- 
den worden, zodat de efficiencies nog erg 
laag zijn (ongeveer 20% van het theore- 





Fig. 6. Het spectrum van de röntgenbron 
Cygnus X-1 opgenomen met de Objective 
Grating Spectrometer van het Einstein- 
observatorium. 
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tisch mogelijke). Deze onvolkomenheid is 
nu overwonnen: de tralies die worden 
voorbereid voor een volgend röntgenob- 
servatorium, de Europese EXOSAT, vol- 
doen volledig aan de theoretische ver- 
wachting. Figuur 6 toont een opname van 
de heldere röntgenbron Cygnus X-1 met 
de Objective Grating Spectrometer (1000 
lijnen/mm). 

De centrale piek is de nulde orde van deze 
intense, door de spiegel- en steunrooster- 
respons schijnbaar uitgebreid lijkende 
puntvormige bron, met aan weerszijden 
de eerste orde van 5 tot 25 Á. Zwak 
zichtbaar is ook het tweede-orde-beeld, 
dat slechts 19% van de eerste orde bedraagt. 
Emissie van röntgenstraling met een golf- 
lengte langer dan 25 À is onderdrukt door 
fotoionisatie van waterstof en helium in 
het interstellaire medium tussen bron en 
waarnemer. De vorm van het continuüm- 
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spectrum samen met de op deze foto niet 
duidelijk zichtbare absorptie-K-kant van 
zuurstof bij 23,3 Á levert informatie op 
over de chemische samenstelling van het 
interstellaire medium. 


Conclusie 

Deze zes plaatjes geven slechts een globale 
indruk van de potentiële mogelijkheden 
van het Einstein-observatorium. 

Vele reeds gepubliceerde resultaten (Qua- 
sars, Deep Sky Survey, radiostelsel Cen- 


taurus A, de discrete-bronbijdrage tot de 
extragalactische achtergrond enz.) zijn on- 
besproken gebleven. Nog meer resultaten 
liggen klaar voor nadere analyse. De satel- 
liet zal werkzaam zijn tot het voorjaar van 
1981. De röntgenastronomie, ontstaan 
door de toevallige ontdekking van de 
bron Scorpius X-1 in 1962, is daarmee de 
kinderjaren ontgroeid en behoort met de 
optische, de UV-, en de radioastronomie 
tot een volwaardige tak van de observatio- 
nele sterrenkunde. 
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De dubbelsterpulsar en de algemene 


relativiteitstheorie 


H.F. Henrichs en E.P.J. van den Heuvel 


Naast de algemene relativiteitstheorie van Einstein is er een aantal ‘concurrerende’ re- 
lativistische theorieën die de zwaartekrachtwisselwerking beschrijven. Elk van deze 


theorieën doet bepaalde voorspellingen die, in principe, aan de waarnemingen getoetst 
kunnen worden. Tot voor kort waren de experimentele onzekerheden in de voorspelde 
effecten nog zodanig groot dat op grond hiervan nog vele van deze zwaartekrachttheo- 
rieën naast elkaar mogelijk bleven. 

Een recente ontdekking in de astronomie heeft er echter toe geleid dat een aantal ‘rela- 
tivistische’ effecten met verrassende nauwkeurigheid gemeten kan worden. Dit artikel 
poogt een nadere achtergrond van deze ontdekking, een radiopulsar:in een dubbelster- 
systeem, te geven. 

Alle aan deze dubbelsterpulsar gemeten verschijnselen blijken met Einsteins theorie in 
overeenstemming te zijn, doch zullen in de nabije toekomst waarschijnlijk de meeste 
andere theorieën doen afvallen. Speciaal het optreden van gravitatiestraling trekt hier- 
bij de aandacht. Men mag verwachten dat in de komende jaren de meetnauwkeurigheid 





nog aanzienlijk zal worden verhoogd. 


Newtons gravitatiewet, geïnspireerd op 
de drie wetten van Kepler, opende de 
mogelijkheid om de banen van hemellicha- 
men nauwkeurig te beschrijven en te voor- 
spellen. Sinds Newtons tijd is het plane- 
tenstelsel steeds bij uitstek het terrein ge- 
weest voor het toetsen van gravitatietheo- 
rieën. De reden hiervoor is vooral dat de 
zon het enige object in onze buurt is dat 
een zeer sterk zwaartekrachtveld bezit, 
waarin de bewegingen van testdeeltjes (de 
planeten) nauwkeurig kunnen worden be- 
studeerd. 


Ook toen Einstein in 1915 de algemene 
relativiteitstheorie had opgesteld, bleek 

al spoedig dat een aantal van de door deze 
theorie voorspelde correcties op Newtons 
theorie (zoals de periheliumdraaiing van 
de mercuriusbaan en de afbuiging van het 
licht langs de zonrand) in het zonnestelsel 
kan worden gemeten. Toch is er nog een 
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aantal door de algemene relativiteitstheo- 
rie voorspelde verschijnselen die niet in 
het zonnestelsel zijn te verifiëren omdat 
het zwaartekrachtveld van de zon hiertoe 
niet voldoende sterk is. Dit geldt in het 
bijzonder voor tweede- en hogere-orde- 
effecten, zoals het optreden van gravita- 
tiestraling. Het blijkt dat een aantal van 
deze hogere-orde-effecten in de zoge- 
naamde dubbelsterpulsar wèl gemeten 
kan worden. Deze pulsar werd in 1974 
ontdekt door Hulse en Taylor! van de 
University of Massachusetts te Amherst. 
Ze ontdekten dat de pulsar, die de naam 
PSR 1913+16 draagt (de getallen geven 
ruwweg de coördinaten aan de hemel aan), 
zich in een nauwe, zeer excentrische (e = 
0,61) baan beweegt om een (onzichtbare) 
begeleider. De omloopperiode is minder 
dan 8 uur, en de afstand tussen de twee 
sterren is van de orde van de straal van de 
zon. ’ 


Pulsars zijn zeer nauwkeurige klokken en 
het feit dat zo’n precisieklok zich in een 
sterk gravitatieveld beweegt, opende de 
mogelijkheid om in een tamelijk kort tijd- 
bestek relativistische effecten te meten. 
Om een indruk te geven van de sterkte 
van de algemeen-relativistische effecten in 
deze dubbelster: een bekend resultaat van 
de algemene relativiteitstheorie is de extra 
periheliumdraaiing van de elliptische baan 
van mercurius rond de zon. De lange as 
van de mercuriusbaan verdraait ongeveer 
43" per eeuw en een nauwkeurige meting 
van dit effect kost vele decennia. In de 
dubbelsterpulsar bedraagt daarentegen de 
waargenomen relativistische draaiing van 
de lange as meer dan 4° per jaar, en deze 
kon reeds binnen twee jaar na de ontdek- 
king van de pulsar nauwkeurig worden ge- 
meten. Dit is uiteraard het gevolg van het 
veel sterkere gravitatieveld in deze dubbel- 
ster vergeleken met dat in ons zonnestel- 
sel. Een maat voor de sterkte van het veld 
is GM/c?r, waarbij G de gravitatieconstan- 
te is, c de lichtsnelheid en r de afstand tot 
het middelpunt van het object met massa 
M. Het blijkt dat deze grootheid voor PSR 
1913+16 ruim 100 maal groter is dan 
voor mercurius. 


De draaiing van de lange as is een eerste- 
orde-correctie op de theorie van Newton. 
De relativistische correcties op de baanbe- 
weging treden omloop na omloop op en 
leiden na vele omlopen op cumulatieve 
wijze tot grote afwijkingen ten opzichte 
van de door Newtonse gravitatie voorspel- 
de posities. Vier jaar na de ontdekking 
van de pulsar waren er zoveel meer en 
nauwkeuriger waarnemingen bekend, dat 
voor het eerst inderdaad een tweede-orde- 
afwijking kon worden gemeten: de om- 
loopperiode blijkt korter te worden. 
Taylor, Fowler en McCulloch® interpre- 
teerden deze metingen aan de hand van 
Einsteins theorie. In dit artikel wordt na- 
der ingegaan op deze interpretatie. 


Radio pulsars 

Radiopulsars® zijn radiobronnen waar- 
van de emissie een duidelijk gepulst ka- 
rakter heeft (zie figuur 1). De eerste radio- 
pulsar werd in 1967 ontdekt door een 
groep in Cambridge, Engeland, geleid 
door A. Hewish, die hiervoor in 1974 de 
Nobelprijs voor Natuurkunde ontving. 
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Drs. H.F. Henrichs (30) studeerde aan de 
Universiteit van Amsterdam waar hij de 
doctoraalexamens experimentele natuur- 
kunde (1975, Zeeman Laboratorium) en 
sterrenkunde (1976) aflegde. Sindsdien 
bereidt hij zich voor op zijn proefschrift 
(aanvankelijk in ZWO-verband) onder 
leiding van prof. dr. E.P.J. van den Heu- 
vel op het Sterrenkundig Instituut van de 
Universiteit van Amsterdam. Onderzoek 
wordt verricht aan fysische processen die 
samenhangen met het massatransport in 
röntgendubbelsterren, waarbij de aandacht 
geconcentreerd ligt rond röntgenpulsars. 





Sindsdien zijn er meer dan 300 pulsars ge- 
vonden. Pulsars blijken een zeer stabiele 
pulsperiode te bezitten met een waarde 
tussen de 0,03 en 5 seconden. De stabili- 
teit van deze ‘klokken’ is van de orde 
ÖP/P = 10°. De algemeen aanvaarde 
verklaring voor de stabiliteit van deze 
klokken is dat de rotatie van een zwaar 
object de pulsperiode bepaalt. 


Bij pulsars blijkt het object een neutro- 
nenster te zijn. De massa is vergelijkbaar 
met die van de zon, terwijl de diameter 
ongeveer 20 km bedraagt. De materie van 
een neutronenster bestaat vrijwel geheel 
uit gedegenereerde neutronen. De gemid- 
delde dichtheid is dus ongeveer die van 
kernmaterie: ruwweg 1017 kg-m”®. Boven- 
dien blijken neutronensterren een extreem 
sterk magnetisch veld te bezitten met gro- 
tendeels een dipoolkarakter. Aan het op- 
pervlak bedraagt dit veld 10° T. De radio- 
emissie blijkt opgewekt te worden in ge- 
bieden in de buurt van de magnetische 
polen. Het pulserende karakter ontstaat 





Fig. 1. Intensiteit vande radiopulsar PSR 
1913+16 als functie van de tijd. Elke pul- 
sar heeft zijn eigen specifieke pulsvorm. 
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Fig. 2. Model voor een radiopulsar. De 
magnetische dipool-as valt niet samen met 
de rotatie-as. De kegelvormige radiobun- 
dels (niet getekend) zijn verbonden met 
de magnetische polen. Alleen waarnemers 
binnen het bereik van de kegelbundels zul- 
len de gepulste radiostraling zien. 





vanwege het niet samenvallen van de rota- 
tie-as en de magnetische dipool-as (zie 
figuur 2). De fysica van de magnetosfeer 
van een pulsar is zeer gecompliceerd en de 
precieze locatie van de opwekking der 
radio-emissie is nog steeds niet duidelijk. 
In elk geval liggen de bronnen ervan niet 
verder dan enkele duizenden kilometers 
van het steroppervlak. De uitgezonden 





Fig. 3. Het Arecibo Observatorium in 
Puerto Rico. Met zijn 305 m diameter is 
deze (vast opgestelde) schotel de grootste 
radiotelescoop ter wereld. De natuurlijke 
vorm van het dal wordt als ondersteuning 
gebruikt. Bijna 40000 aluminium panelen 
vormen een sferisch oppervlak dat minder 
dan 3 mm afwijkt van de ideale vorm. De 
panelen zijn geperforeerd opdat de vege- 
tatie eronder in leven blijft en daardoor 
de bodem stabiel houdt. De ontvanger in 
het brandpunt kan (zeer nauwkeurig) ver- 
plaatst worden zodat een gebied van 40° 
diameter aan de hemel bestudeerd kan 
worden. 











Prof. dr. E.P.J. van den Heuvel (39) is 
hoogleraar sterrenkunde aan de Univer- 
siteit van Amsterdam. Hij studeerde van 
1957 tot 1962 natuur-en sterrenkunde te 
Utrecht, waar hij in 1968 promoveerde 
op een studie over de rotatie van sterren. 
Hij werkte aan de RU-Utrecht, de 

Vrije Universiteit van Brussel (waar hij 
sinds 1970 Buitengewoon Hoogleraar is), 
de University of California te Santa Cruz 
en het Institute for Advanced Study te 
Princeton. Sinds 1974 is hij werkzaam te 
Amsterdam. Zijn onderzoeksterrein ligt 
op het gebied van de sterevolutie en de 
hoge-energie-astrofysica. 





radiostraling is slechts een kleine fractie 
van de totaal door de neutronenster uit- 
gezonden energieflux. Het grootste deel 
van deze flux wordt geëmitteerd in de 
vorm van laagfreqùente magnetische di- 
poolstraling opgewekt door de voortplan- 
ting van het ronddraaiende magnetische 
veld in de ruimte. 

De energiebron van een pulsar is de rota- 
tie van de neutronenster: de emissie van 
de magnetische dipoolstraling gaat ten 
koste van de rotatie-energie. We zien dan 
ook dat praktisch alle radiopulsars in de 
loop der tijd steeds langzamer gaan 
draaien en dus een steeds langere pulspe- 
riode krijgen. 


De dubbelsterpulsar 
PSR 1913+16, door Hulse en Taylor in 


juli 1974 met de reuzentelescoop van Are- 
cibo in Puerto Rico“ (zie figuur 3) ont- 
dekt, heeft een periode van slechts 0,059 
s, de op één na kortste periode van alle 
bekende radiopulsars. De metingen te 
Arecibo werden uitgevoerd bij 430 MHz 
en 1410 MHz (respectievelijk 70 en 21 
cm). De gebruikte ontvanger detecteert 
de radiostraling in 32 aan elkaar grenzende 
frequentiebanden, elk met een band- 
breedte van 250 kHz. 


De keuze van deze ontvangsttechniek is 
bepaald door (afgezien van de pulsduur 
zelf) het feit dat het pulsarsignaal disper- 
sie ondergaat in het interstellaire medium, 
als gevolg van de daar aanwezige vrije elek- 
tronen. Deze dispersie maakt dat de puls 
bij een hogere frequentie eerder arriveert 
dan bij een lagere. Door nu de ontvangen 
pulssignalen in de respectieve kanalen 
naar toenemende frequentie op de juiste 
wijze te vertragen, kan men voor de inter- 
stellaire dispersie corrigeren, en de aan- 
komsttijden bij verschillende frequenties 
laten samenvallen, waarna de signalen op- 
geteld kunnen worden. Zodoende wordt 
gecompenseerd voor het verlies aan ge- 
voeligheid dat anders door de interstellaire 
dispersie zou worden veroorzaakt. Deze 
dispersie hangt af van de afstand van de 
pulsar. Dat voor iedere pulsar de ontvan- 
ger aldus op specifieke wijze moet wor- 
den geoptimaliseerd is een bijkomstige 
moeilijkheid bij het opsporen van vooral 
zwakke pulsars. 


Reeds een paar uur na de ontdekking van 
PSR 1913+16 bleek dat de periode niet 
constant is, maar op systematische wijze 
varieert tussen 0,058967 en 0,059045 s 
(zie figuur 4). De interpretatie was onmid- 
dellijk dat deze pulsar zich in een baan 
rond een ander object beweegt. Dat de 
aankomsttijden erg slecht op een sinus- 
kromme passen maakte duidelijk dat het 
om een zeer excentrische baan gaat. De 
complete afwezigheid van eclipsen wijst 
erop dat de onzichtbare begeleider kleiner 
moet zijn dan een ster die in het stadium 
van waterstoffusie verkeert (een zoge- 
naamde hoofdreeksster, zoals de zon). 





Fig. 4. De aankomsttijden (vertaald in 
radiale snelheden) van de dubbelsterpulsar 
als functie van de baan fase. 
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Meer conclusies waren na een maand niet 
gerechtvaardigd. 


Na anderhalf jaar, in het begin van 1976, 
was het voor het eerst mogelijk om de 
verandering in de lengte van het periastron 
met redelijke nauwkeurigheid te meten 
(zie tabel I). Figuur 5 illustreert de defi- 
nities van de diverse baanelementen. 


De analyse van de aankomsttijden der 
pulsen 

In de loop der jaren heeft men de aan- 
komsttijden van de pulsen kunnen meten 
met een relatieve nauwkeurigheid van 50 
us. Van de functie die de aankomsttijden 
van de pulsen bij de aarde in de loop van 
de tijd weergeeft, verwacht men — op 
theoretische gronden — dat zij vastgelegd 
js door een 13-tal parameters. Deze zijn: 
e De coördinaten a en Ô van de pulsar 
aan de hemel. Het dopplereffect, als ge- 
volg van de baanbeweging van de aarde 
rond de zon en de rotatie van de aarde 
om haar as, veroorzaakt veranderingen in 
de aankomsttijden, welke in de loop van 
het jaar variëren. De wijze van deze varia- 
tie hangt van a en Ô af. 

e De pulsperiode P en diens afgeleide 
naar de tijd, P. 

e De vijf keplerbaan-parameters (in fi- 
guur 5 is hun betekenis weergegeven). 
Deze veronderstellen newtonse mecha- 
nica en bepalen de hoofdmoot van de 
periodieke dopplerverschuivingen die het 
pulssignaal ondergaat. Tenslotte: 


e Vier parameters van relativistische oor- 


sprong, die de correcties op de newtonse 
baanbeweging weergeven. (Het meenemen 
van nog meer relativistische parameters is 
niet nodig daar deze in alle theorieën or- 
des van grootte kleiner zijn dan deze vier) 


Verder moet men nog een (bekende) alge- 
meen-relativistische klok-correctie toe pas- 
sen om in rekening te brengen dat het op 
aarde gevoelde zwaartekrachtveld van de 
zon in de loop van het jaar varieert. 


Elk der 13 parameters beïnvloedt de aan- 
komsttijden van de pulsen op zijn eigen 
wijze en, door de waarden van deze para- 





Fig. 5. De baan van de dubbelsterpulsar 
in de ruimte. De symbolen stemmen over- 
een met die in tabel I. De lengte van het 
periastron is wg. De halve lange as van 

de elliptische baan isa, 





meters te variëren, bepaalt men (met een 
kleinste-kwadratenmethode) die combi- 
natie van parameterwaarden waarvoor een 
optimale aanpassing bij de waargenomen 
aankomsttijden wordt verkregen. In tabel 
l zijn de aldus verkregen waarden van de 
13 parameters weergegeven met hun bij- 
behorende meetonzekerheden zoals be- 
paald door Taylor en medewerkers. De 
met deze 13 waarden verkregen functie, 
die de aankomsttijden der pulsen voor- 
spelt, is zo goed dat men hiermee de puls- 
aankomsttijden 8 maanden van tevoren 
voorspellen kan met een nauwkeurigheid 
beter dan 100 microseconden! 


De relativistische parameters 

De vier relativistische parameters zijn: 

e War, de algemeen -relativistische hoek- 
snelheid van periastrondraaiing; 

e ‘y,welke een maat is voor de (variabele) 
vertraging die de puls ondergaat in het 
dubbelstersysteem als gevolg van de gravi- 
tatieroodverschuiving gecombineerd met 
het transversale dopplereffect. 

e sini, waarbij i de inclinatiehoek van 
het baanvlak is ten opzichte van de hemel- 
bol (zie figuur 5). Deze beinvloedt de ver- 
traging van het pulssignaal in het zwaarte- 
krachtveld van de begeleider. 

e P,de verandering van de dubbelster- 
periode (als gevolg van emissie van gravi- 
tatiestraling). 

Deze vier parameters zijn thans alle met 
een 95% of hogere betrouwbaarheid be- 





Tabel l. De 13 parameters voor de dubbelsterpulsar PSR 191 3+16. 


gf 


Rechte klimming (1950,0) 

Declinatie 

Pulsperiode 

Tijdsafgeleide van periode 

Geprojecteerde halve lange as _ 

Excentriciteit 

Baanperiode 

Lengte van het periastron 

Tijd van periastronpassage 

Transversaal dopplereffect en 
gravitationele roodverschuiving 

Inclinatiehoek 

Tijdsafgeleide van baanperiode 


Snelheid van periastrondraaiing 


a =19"13"125474(4) 
8 _=16°01'08/02(6) 

P =0,059029995269%2) s 
P =(8,64*0,02) x 107° 


a, sini =2,3424(7) lichtseconde 


e _=0,617155(7) 

Po =27906,98172(5) s 

wo = 178°864(2) 

To =JD 2442321 433206(1) 


\ =0,0047(7)s 

sini =0,81(16) 
P, =(-32*+0,6)x 107% 
O =4°226(2) jaar" 
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y = 0,002951 (Mp) (Mp +2 Mo) Mp +M‚)” 213 


sini __ =0,5083 (M‚) * (Mp +M)2/3 
P, = 1,70 Xx 102 (M‚M) (My + My) 1/3 
WAR == 211 (Mp + Mo )2/3 jaar 





Tabel II. De massa-afhankelijkheid van de relativistische parameters. 


(alle massa’s in zonne-eenheden) 





paald. Elk van deze vier parameters blijkt 
een functie te zijn van de onbekende 
massa’s M‚ en My van de pulsar en zijn 
begeleider (en van de waargenomen kep- 
lerbaanparameters, die zeer nauwkeurig 
bekend zijn). De precieze aard van deze 
vier functionele verbanden hangt nog af 
van de gebruikte gravitatietheorie. Voor 
het geval van de algemene relativiteits- 
theorie zijn deze vier functies in tabel [I 
weergegeven. Meting van de vier parame- 
ters levert vier vergelijkingen op waaruit 
men de onbekende massa’s M‚, en My, zal 
moeten bepalen. 


Daar er meer vergelijkingen dan onbeken- 
den zijn is het de vraag of de M‚- en M- 
waarden die door oplossing van de ver- 
schillende vergelijkingen verkregen wor- 
den binnen de grenzen van de meetnauw- 
keurigheid met elkaar in overeenstemming 
zijn. Indien dit het geval is, is de toege- 
paste theorie bruikbaar. Het blijkt dat dit 
voor de algemene relativiteitstheorie inder- 
daad het geval is. Figuur 6 toont de ge- 
bieden in het (M‚, My) vlak welke door 


Fig. 6. De vier gemeten relativistische 
parameters leveren vier vergelijkingen op 
voor de massa’s M), en My, van de pulsar 
en zijn begeleider (zie tabel II). Van deze 
vier is AR met zeer grote nauwkeurig- 
heid gemeten; daar deze alleen van de 
som M= M), + My, afhangt, levert hij de 
teruglopende lijn op. Voor de andere 3 
(hogere-orde relativistische) parameters 
is het 95%-betrouwbaarheidsinterval aan- 
gegeven, hetgeen gezamenlijk het gear- 
ceerde gebiedje levert. Men mag stellen 
dat deze twee gebieden niet significant 
verschillen, gezien de onzekere grenzen 
van het gearceerde gebied. Samen levert 
dit M,=1,39*0,15 Mo en My =144+ 
0,15 Mo. 


leider (M ) 
S 


ND 
o 





Massa van de bege 


Ó 1.0 30 


2.0 
Massa van de Pulsar (Mo ) 


Nederlands Tijdschrift voor Natuurkunde A46 (1) 1980 


de gemeten betrouwbaarheidsgebieden 
van de vier parameters worden vastgelegd. 
Het eerste-orde-effect vo is met zeer 

grote nauwkeurigheid gemeten en levert 
de teruglopende rechte op; in dit gebied 
is de onzekerheid kleiner dan de dikte van 
de lijn. Het gecombineerde gebied, opge- 
leverd door de drie hogere-orde-parameters 
is gearceerd. De grenzen van dit gebiedje 
zijn echter nogal onzeker. Men neemt dan 
ook aan dat deze twee gebieden niet sig- 
nificant verschillen. Alles samen levert: 
My =1,39+ 0,15 Mo en Mo =1,44# 
0,15 Mo.(Mo is de massa van de zon.) 
Merk op dat de pulsar en zijn begeleider 
vrijwel dezelfde massa hebben. 


De aard van de begeleider 

In bovenstaande analyse is veronderstelt 
dat de begeleider zelf geen bijdrage levert 
tot de correcties op de keplerbaanpara- 
meters. Dit is alleen het geval indien deze 
begeleider zelf een compacte ster is, dat 
wil zeggen een neutronenster of een witte 
dwerg. Wanneer de begeleider een minder 
compacte ster is, zal door getijdewerking 
een zodanige vervorming van deze ster op- 
treden dat het bijbehorende gravitationele 
quadrupoolmoment een aanzienlijke bij- 
drage tot w gaat geven. Op deze gronden 
kan men een gewone waterstofrijke hoofd- 
reeksster (zoals de zon) direct uitsluiten 
omdat deze een periastrondraaiing van 
meer dan 3° per omloop zou leveren. Het 
quadrupoolmoment van een heliumster 
zou de gemeten & nog wel kunnen op- 
brengen. Echter, omdat in dat geval de re- 
lativistische bijdrage tot w slechts een 
kleine fractie daarvan kan zijn, zou een 
dergelijke heliumster slechts een geringe 
massa mogen hebben (slechts enkele tien- 
den van een zonmassa). Dit levert dan 
weer problemen op met de gemeten waar- 
de van P‚, die meer dan zesmaal groter is 
dan men voor de zwaarst mogelijke he- 
liumster kan verwachten. De enige moge- 
lijkheden die voor de begeleider overb lij- 
ven zijn dan: een neutronenster en (welis- 
waar minder waarschijnlijk) een zware, 
langzaam roterende witte dwerg. (Een 


snel roterende afgeplatte witte dwerg zou | 


door zijn quadrupoolmoment een positie- 
ve P‚ geven.) 


Een recente optische identificatie” sluit 
een neutronenster uit. Hieraan wordt 
door de meeste astronomen geen geloof 
gehecht daar deze identificatie blijkt te 
berusten op een aanwijsbare, echter nog 
niet door iedereen begrepen fout in de ge- 
bruikte coördinatentransformatie. 





Fig. 7. Mogelijke evolutieweg naar een 
dubbelsterpulsar®. In een dubbelster met 
massa’s 20en 8Mo zal de zwaarste ster 
het snelst evolueren (a) en massa overdra- 
gen aan de begeleider (b). De achtergeble- 
ven heliumster (c) explodeert als een super- 
nova (d). Het restant is een neutronen- 
ster die eerst een röntgenpulsar zal wor- 
den (e) door de vangst van sterrewindma- 
terie en uiteindelijk als radiopulsar zal 
achterblijven (h). Intussen is de begeleider 
ook als heliumster aan zijn eind gekomen 
(g)en indien het systeem deze tweede 
supernova-explosie overleeft zal een dub- 


belsterpulsar kunnen overblijven (h), waar- 


van de jonge zowel als de oude neutronen- 
ster een radiopulsar kan zijn. 
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dubbelster pulsar 


Er zijn ook evolutionaire overwegingen 
die een neutronenster als begeleider het 
meest aannemelijk maken (zie figuur 7). 
Berekeningen uitgevoerd in het Sterren- 
kundig Instituut van de Universiteit van 
Amsterdam tonen dat de voorloper van 
een dubbelsterpulsar een röntgendubbel- 
ster geweest zou kunnen zijn, bestaande 
uit een neutronenster die röntgenstraling 
uitzendt (de tegenwoordige radiopulsar) 
en een zware, zichtbare ster. Wanneer 
deze laatste ster aan het eind van zijn 
leven explodeert als een supernova zal er 
een tweede neutronenster achterblijven. 
Men verwacht dat in de meeste gevallen. 
de dubbelster als gevolg van de supernova- 
explosie ontbonden zal worden. Maar àls 
de twee neutronensterren bijeen blijven 
wordt inderdaad een zeer excentrische 
baan verwacht. Tot nu toe is dit de enig 
mogelijke evolutieweg gebleken volgens 
welke een dubbelsterpulsar gevormd kan 
worden. 


Gravitatiestraling 

Taylor en zijn medewerkers voerden nog 
een test uit om de significantie van de af- 
name van de baanperiode als gevolg van 
gravitatiestraling na te gaan. Hiertoe 
maakten ze de veronderstelling dat de 
baanperiode P‚, in de loop van de tijd niet 
veranderd zou zijn. Onder deze veronder- 
stelling losten ze wederom alle baanpara- 
meters op, waaronder ook het tijdstip van 
periastrondoorgang. Figuur 8 toont de al- 
dus uit de waarnemingen afgeleide mo- 
menten van periastrondoorgang. De rech- 
te lijn stelt het ware tijdstip voor van 
doorgang in 1974. De parabool stelt het 
teruglopen van de berekende momenten 
van doorgang voor indien de baanperiode 
afneemt als gevolg van gravitatiestraling 
volgens Einsteins theorie (indien beide 
sterren neutronensterren zijn met massa 
1,41 Mo ). De figuur toont duidelijk dat 
de afname van Pp in fraaie overeenstem- 


Fig. 8. De periastrondoorgang als functie 
van de tijd. De punten geven het moment 
van periastrondoorgang zoals die geweest 
zou zijn als men bij de analyse de door- 
gang van 1974 (de rechte lijn) had ge- 
bruikt en de baanperiode constant had 
gehouden. De parabool komt overeen met 
de algemeen-relativistische voorspelde af- 
name in de baanperiode als gevolg van 
gravitatiestraling. 
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ming is met de algemene relativiteits- 
theorie. 


Men bedenke dat deze gravitatiestraling 
quadropoolstraling is. Diverse gravitatie- 
theorieën voorspellen ook het optreden 
van gravitationele dipoolstraling die (in 
het algemeen) sterker is dan de quadru- 
poolstraling. Helaas is in deze theorieën 
de intensiteit van de dipoolstraling pre- 
cies nul indien de massa’s van de twee 
sterren aan elkaar gelijk zijn, hetgeen hier 
niet uitgesloten is. Evenwel, de massabe- 
palingen zijn verricht onder de vooronder- 
stelling dat de algemene relativiteitstheo- 
rie (die geen dipoolstraling voorspelt) gel- 
dig is. Een analyse op grond van een an- 
dere theorie zal in het algemeen andere 
parameterwaarden opleveren, maar is nog 
niet gedaan. Het lijkt echter weinig waar- 
schijnlijk dat zo’n analyse precies dezelf- 
de waarden voor de massa’s zal opleveren. 


Kanttekeningen 

Er is kritiek geuit op de juistheid van de 
formule voor gravitatiestraling die ten 
grondslag ligt aan de door Tayloret al. uit- 
gevoerde analyse. In het bijzonder voert 
men aan dat de formule (Landau en Lif- 
shitz geven de afleiding ervan in hun boek 
als opgave) niet exact is, doch slechts een 
benadering. De benadering ligt in het feit 
dat men in de berekening aanneemt dat 
de baan niet krimpt tijdens het uitzenden 
van de gravitatiestraling. Men brengt ver- 
volgens dit energieverlies in rekening door 
de baan te laten krimpen. Dit is natuurlijk 
geen consistente aanpak van het probleem, 
maar tot nu toe ondervindt een exacte be- 
rekening zowel conceptuele als mathema- 
tische moeilijkheden. Toch geloven de 
meeste auteurs dat het slechts om een ho- 
gere-orde correctie gaat, veel kleiner dan 
het voorspelde en waargenomen effect 
zelf. Hier ligt nog een uitdaging voor de 
mathematische fysici. 


Verder kan worden geconstateerd dat in 
figuur 8 het eerste punt een nogal domi- 
nerend effect heeft op de kromming. Een 
paar jaar extra meten zal hier echter de 
twijfel kunnen wegnemen. 


Op grond van bovenstaande kan men ver- 
wachten dat de dubbelsterpulsar hèt stan- 
daard testobject zal worden voor gravita- 
tietheorieën. Tot nu toe is Einsteins alge- 
mene relativiteitstheorie in ieder geval in 
staat gebleken om alle gemeten niet-new- 
tonse effecten kwalitatief èn kwantitatief 
te beschrijven. De komende jaren zullen 
leren welke andere gravitatietheorieën de 
toets aan de dubbelsterpulsar overleven. 


De auteurs zijn drs. P. Hut zeer erkente- 
lijk voor een verhelderende discussie over 
gravitatiestraling. 
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Clusters van sterrenstelsels 


N. Roos 


Clusters van sterrenstelsels zijn de grootste objecten in ons expanderend heelal. Op een 
schaal groter dan 10 megaparsec (1 Mpc = 3,1 x 10° m) kan de verdeling van materie 
als homogeen worden beschouwd. De snelheid waarmee sterrenstelsels (St’s ) met onder- 
linge afstand d (> 10 Mpc) van elkaar vandaan vliegen wordt dan gegeven door Hubble’s 
wet: v= Hod waarin H, de zgn. Hubble-constante is*. Op kleinere schaal blijken ster- 
renstelsels bij voorkeur in groepen of clusters van twee tot enkele honderden voor te 
komen. Het gemiddeld aantal heldere stelsels** per Mpc° bedraagt ongeveer 2x10”*. 
Wanneer deze willekeurig in de ruimte waren verdeeld zou de gemiddelde afstand van 
een stelsel tot zijn naaste buur ongeveer 4,5 Mpc bedragen. In werkelijkheid is die af- 
stand ongeveer 1 Mpc. De afstand van ons melkwegstelsel tot zijn naaste buur de An- 
dromeda-nevel (M31) bijvoorbeeld is 0,7 Mpc; bovendien vallen ons melkwegstelsel en 
M31 naar elkaar toe met een snelheid van 50-100 km/s en niet van elkaar af zoals zou 
moeten in een uniform expanderend heelal. Samen met nog wat kleinere stelsels vor- 


men ze een groep, de zgn. Lokale Groep van sterrenstelsels. 
In dit artikel zal het gaan om de vraag hoe groepen of clusters van sterrenstelsels ont- 
staan en hoe hun vormingsgeschiedenis samenhangt met die van sterrenstelsels. 


Vorming van sterrenstelsels en clusters 
van St’s 

Sinds de ontdekking van de expansie van 
het heelal door Hubble, de ontdekking en 
interpretatie van de 3 K-achtergrondstra- 
ling, en de juiste voorspelling van de kos- 
mische heliumabundantie staat het ‘Big 
Bang’ model stevig in zijn schoenen! 

In dit kosmologisch model is het heelal 
geëvolueerd uit een zeer dichte, hete toe- 
stand, waarin straling en materie met 
elkaar in thermisch evenwicht waren. 
Door de expansie koelden straling en 
materie snel af. Het thermische evenwicht 
tussen hen werd verbroken bij 

ze 1000*** toen elektronen en pro- 
tonen (re)combineerden. Men veronder- 
stelt dat op dat moment fluctuaties met 
kleine amplitude aanwezig waren in de 
verdeling van materie. De expansie in zo’n 
dichtheidsverhoging wordt sterker afge- 
remd dan gemiddeld. De fluctuaties heb- 
ben daardoor de neiging te groeien (gravi- 
tatie-instabiliteit) en in sommige gebieden 
zal de expansie zelfs worden gestopt, 
waarna het gebied samentrekt totdat een 
evenwichtssituatie is bereikt. De mate van 


*_We zullen in dit artikel aannemen dat H, = 
50 km-s"*. Mpcr*. 

** Onder ‘heldere’ sterrenstelsels zullen we ver- 
staan: St's met optische helderheid 

Le 310 Ee 

Lois de helderheid van de zon (= 3,86x10° 
watt). De massa van de zon Mo = 2 x 10° kg. 
*** In Big Bang-kosmologieën hangt de rood- 
verschuivingsparameter z = AA/X (AA = XA, — À, 
A is de golflengte van een foton op het moment 
van uitzendingen A, de golflengte bij ontvangst) 
direct samen met de schaal R van het heelal vol- 
gens AIR, e(iteh*, 

Het verband tussen R en de leeftijd van het 
heelal f hangt af van de dichtheid in het heelal, 
maar wordt (in ieder geval in het vroege heelal) 
goed beschreven door R/R, = (t/t,)? ”. De 
huidige leeftijd van het heelal, f,, wordt geschat 
op 10 à 20 miljard jaar. 
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contractie wordt bepaald door de hoeveel- 
heid energie die het gebied kwijt heeft 
kunnen raken. Zo kunnen massaconcen- 
traties als sterren, sterrenstelsels en clus- 
ters van sterrenstelsels zijn ontstaan. 

Het is mogelijk dat gaswolken met de af- 
metingen van grote clusters (M ongeveer 
10 # # M,) zich als eerste losmaakten 
uit de kosmische expansie, en vervolgens 
instortten. De contractie hoeft zeker niet 
in 3 onafhankelijke richtingen tegelijker- 
tijd te gebeuren. In de afgeplatte gebieden 
van hoge dichtheid (“pancakes”) die zo 
ontstonden zouden de sterrenstelsels (St’s) 
gevormd kunnen zijn? 


In het alternatieve scenario gaat men er- 
van uit dat ook dichtheidsfluctuaties die 
minder massa omvatten kunnen groeien. 
Men neemt aan dat het spectrum van 
dichtheidsfluctuaties ôp/p na recombina- 
tie voldeed aan 
ope Ml, well3=23. 

In dit geval zullen de fluctuaties die de 
kleinste massa omvatten het eerst omke- 
ren en instorten. Eerst worden de St’s ge- 
vormd. Vervolgens ontstaan de groepen 
en clusters uit de fluctuaties in de verde- 
ling van St’s. 

Een belangrijke vraag is die naar het on- 
derscheid tussen St’s en clusters van 

St’s. Waarom komen er alleen St’s voor 
met MS 10'? Mo en RS 50 kpc*. Het 
blijkt nu dat deze karakteristieke massa 
en straal direct volgen uit de eis dat een 
gaswolk door middel van vrij-vrij-inter- 


* Uitzonderingen op deze regel vormen zgn. 
cD (central Dominant)stelsels, welke waar- 
schijnlijk een bijzondere vormingsgeschiedenis 
achter de rug hebben (zie volgende paragraaf). 
De massa van deze stelsels is ongeveer 10'° Mo. 





Drs. Nico Roos (29) studeerde theoreti- 
sche natuurkunde aan de VU te Amster- 
dam. 

Sinds 1977 is hij verbonden aan de Rijks- 
universiteit te Leiden (vakgroep sterren- 
kunde) waar hij bezig is met een promotie- 
onderzoek aan clusters van sterrenstelsels. 


acties tussen protonen en elektronen zijn 
thermische energie zo snel kan uitzenden 
dat ongehinderde contractie onder in- 
vloed van zijn zwaartekracht kan plaats- 
vinden. De zo gevonden massa en straal 
van 10!* Me resp. 75 kpc zijn volledig uit 
te drukken in fundamentele fysische con- 
stanten?. Dissipatie van energie heeft 
zeer waarschijnlijk een belangrijke rol ge- 
speeld bij de vorming van St’s, terwijl de 
vorming van clusters alleen is bepaald 
door de zwaartekracht die de St’s op 
elkaar uitoefenen. 


Een belangrijk succes voor dit model is 
dat in numerieke simulaties van een ex- 
panderend heelal bestaande uit ongeveer 
1000 massapunten de waargenomen clus- 
tering goed wordt gereproduceerd. Een - 
belangrijke statistische maat voor de clus- 
tering van St’s is de zgn. twee-puntscorre- 
latiefunctie £(r). Deze functie kan als 
volgt worden gedefinieerd: 


n(r)=n 41 +E(r)}, 


waarbij n de gemiddelde dichtheid aan 
St’s is en n(r) de gemiddelde dichtheid op 
een afstand r van een St. Uit de waarne- 
mingen volgt dat £(r) over het hele gebied 
waarover zij is bepaald, 10 kpc <r <10 
Mpe, goed wordt beschreven door een 
eenvoudige machtswet® : 


E(r) = (ro/r)"°®, 


Het is dus niet mogelijk over een karakte- 
ristieke grootte van clusters te spreken. 
Bij de kwantitatieve beschrijving van de 
clustering van St’s vormt de zgn. ‘multi- 
plicity function’ een goede aanvulling op 
de twee-puntscorrelatiefunctie. Deze ver- 
delingsfunctie geeft het aantal clusters met 
totale optische helderheid L (hij neemt 
af met L volgens L2*0°3). 

De correlatiefunctie die ontstaat in de nu- 
merieke simulatie van een expanderend 
heelal heeft dezelfde vorm als de waarge- 
nomen functie, en ook de waargenomen 


en berekende multiplicity-functies komen 
41 


Yo = 10 Mpc. 


goed overeen“ . De mate van clustering 
aan het eind van de simulatie, bepaald 
door r‚,‚ hangt af van de waarde van z 
waarop de berekening wordt gestart, of- 
wel de roodverschuiving waarop St’s wor- 
den gevormd. Dat blijkt ongeveer bij z = 
20 te moeten gebeuren. 

Begincondities zoals de vorm van het fluc- 
tuatiespectrum en de dichtheid van het 
heelal blijken weliswaar geen sterke, maar 
toch meetbare invloed te hebben op de 
eindsituatie, en het ziet ernaar uit dat 
verdere studie op dit gebied veel kan ver- 
duidelijken over de omstandigheden in 
het vroege heelal en de ontstaansgeschie- 
denis van St’s en clusters van St’s. 


De rijke clusters 


Classificatie 

Abell heeft in 1958 een catalogus samen- 
gesteld van 2712 rijke clusters van sterren- 
stelsels, geïdentificeerd op platen van de 
Palomar Sky Survey. Deze clusters heb- 
ben roodverschuiving z = 0,02 - 0,2 en 
bevatten 20 tot 100 heldere leden binnen 
een afstand van 3 Mpc tot het clustercen- 
trum. Heldere sterrenstelsels kunnen wor- 
den onderscheiden in 3 typen: Spiraalstel- 
sels (Sp), SO-stelsels, welke beide bestaan 
uit een platte roterende schijf met een 
rond of ellipsvormig lichaam (de ‘bulge’) in 
het centrum, en elliptische stelsels (E). 
SO-stelsels bevatten evenals elliptische 
stelsels zo goed als geen (neutraal) gas 

en hebben over het algemeen wat grotere 
bulges dan spiraalstelsels. Er bestaat een 
verband tussen de structuur van rijke 
clusters en de fractie spiraalstelsels die ze 
bevatten”. Sommige clusters bevatten 
overwegend SP-stelsels (ongeveer 60%), 
25% SO-stelsels en 15% E-stelsels (deze 


Fig. 1. Een deel van de Hercules-cluster. 
Dit is een typisch voorbeeld van een 
(open) onregelmatige, spiraalrijke cluster. 
Opvallend veel leden zijn vervormd door 
wisselwerking met een ander stelsel. De 
foto is 33’ breed. Op de afstand van 
Hercules (-- 200 Mpc) komt dat overeen 
met — 2 Mpc. 


percentages vindt men ook in de gebieden 
buiten de rijke clusters) De structuur van 
deze clusters is meestal onregelmatig en 
er zijn vaak verschillende massaconcentra- 
ties in te onderscheiden. Hun tegen- 
hangers hebben meestal een regelmatige, 
bijna sferisch-symmetrische vorm met een 
goed gedefinieerd centrum waarin de 
dichtheid aan sterrenstelsels heel hoog kan 
zijn (1000 Mpc ® ). Deze clusters bestaan 
overwegend uit SO- en Estelsels en bevat- 
ten weinig Sp-stelsels. In alle clusters 
neemt de fractie elliptische stelsels toe 
met de dichtheid aan sterrenstelsels. 

In de buurt van het centrum van veel re- 
gelmatige clusters met een hoog percentage 
E-stelsels, vindt men vaak een enorm 
groot en helder elliptisch stelsel: een zo- 
genaamd cD (central Dominant)-stelsel. 
De massa van deze stelsels kan gemakke- 
lijk 10/° M bedragen. 


Instorting 

In numerieke simulaties van een expan- 
derend heelal, zoals in de vorige paragraaf 
genoemd, manifesteert clustering zich het 
eerst op kleine en vervolgens op steeds 
grotere schaal. Simulaties van de evolutie 
van rijke clusters laten ook sterke subclus- 
tering zien tijdens de expansie en de daar- 
op volgende instorting van de cluster. 
Gedurende deze periode, welke ongeveer 
te = (Gp) “duurt (G = gravitatieconstan- 
te, p = dichtheid van de cluster) krijgen 
de sterrenstelsels steeds grotere ‘random’ 
snelheden door wisselwerking met het 
voortdurend veranderend gravitatieveld. 
Dit proces staat bekend als violent rela- 
xation. Al na ongeveer 3/2 t, is er een 
evenwicht tussen kinetische (7) en poten- 
tiële energie (V/) van de cluster ontstaan, 
beschreven door het viriaaltheorema: 27 
+ V = 0. De subclustering is dan uitge- 
wist en de cluster heeft een regelmatig 
dichtheidsverloop. Dit suggereert dat 
twee extreme clustertypen uit de vorige 
paragraaf verschillende evolutiestadia re- 
presenteren. 


Missing mass 
Uit metingen van de snelheden van een 
aantal sterrenstelsels in een cluster kan 





met behulp van het viriaaltheorema de 
clustermassa worden geschat. Deze blijkt 
dan ongeveer een factor 10 groter te zijn 
dan de som van de massa’s van de ster- 
renstelsels. Alleen de massa in het optisch 
zichtbare deel van de stelsels wordt dan 
meegerekend. Metingen van de rotatie- 
snelheid van neutraal waterstof buiten het 
zichtbare deel van de schijf van spiraalstel- 
sels wijzen erop dat deze stelsels zijn om- 
ringd door donkere zware halo’s. Deze 
halo’s zouden kunnen bestaan uit objec- 
ten die al in een vroeg stadium van de vor- 
ming van de stelsels ontstonden. De massa 
M‚ van de stelsels blijkt evenredig met de 
straal R, toe te nemen. Wanneer nu de 
‘missing mass’ van de clusters aan de ster- 
renstelsels vastzit, zijn deze stelsels onge- 
veer 10 maal zo groot en zwaar als men 
aanvankelijk dacht. 


Botsingen, en de vorming van elliptische 
stelsels 

Zeker wanneer St’s zijn omringd door 
massieve halo’s zullen botsingen tussen 
St’s veel voorkomen in rijke clusters. De 
gemiddelde tijd tj die verloopt tussen 
2 opeenvolgende botsingen van een stel- 
sel kan worden geschat met 

fest Zr TR; Nei Oc waarin nej de dicht- 
heid aan St’s en o,; de snelheidsdispersie 
van de stelsels in de cluster is. In het cen- 
trale deel van een dichte, regelmatige clus- 
ter (naj — 10° Mpc *, 0.1 = 10° km/s) 
is dat typisch 10° tot 10’ jaar voor een 
stelsel met R‚ — 10 kpc. Voor grotere 
stelsels (halo’s) is die tijd aanzienlijk kor- 
ter. Deze botsingstijd moet natuurlijk 
worden vergeleken met de leeftijd van 
zo’n cluster welke — 5 x 10° jaar be- 


draagt. 


Tijdens een ontmoeting tussen 2 sterren- 
stelsels met relatieve snelheid op grote af- 
stand v en afstand van dichtste nadering 
p, wordt op de sterren een getijde- 

kracht uitgeoefend, welke de inter- 

ne energie U van het stelsel ver- 

hoogt. Wanneer v >> 09, de interne snel- 
heidsdispersie van het stelsel, en p >R 
volgt met behulp van de aanname dat de 
sterren zich tijdens de ontmoeting niet 
ten opzichte van het centrum van het stel- 
sel verplaatsen (de zgn. impulsbenadering) 
dat AU/U « (Re/p)' (og/v)* . Botsingen 
met lage v en p hebben dus het grootste 
effect. Bij lage snelheden is de impuls- 
benadering niet betrouwbaar meer. Er tre- 
den dan gecompliceerde vervormingen 
van het stelsel op tijdens de botsing welke 
alleen numeriek goed gesimuleerd kunnen 
worden. De toename van de interne ener- 
gie gaat ten koste van de baanenergie en 
de botsingen worden sterk inelastisch bij 
lage snelheden. Uit numerieke berekenin- 
gen blijkt dat twee identieke stelsels, die 
botsen met p <R, en v < og, elkaar bin- 
nen korte tijd na de eerste ontmoeting 
invangen en een nieuw stelsel vormen met 
OZ Og. 


en 
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De inelasticiteit van sterrenstelsels kan ge- 
makkelijk op een schematische manier 
worden ingebouwd in numerieke simula- 
ties. Dergelijke berekeningen tonen aan 
dat tijdens de vorming van clusters, wan- 
neer Ocj nog laag is, een aanzienlijke frac- 
tie van de sterrenstelsels naar zijn meest 
nabije buur toevalt en wordt ingevangen. 
In Arp’s ‘Atlas of peculiar galaxies’, ko- 
men veel dubbelstelsels voor die blijkens 
hun sterk afwijkende uiterlijk grote ge- 
tijdekrachten op elkaar uitoefenen. In 
1972 werd voor het eerst door de gebroe- 
ders Toomre gesuggereerd dat deze stel- 
sels elkaar binnen zo’n 5x10® jaar zouden 
invangen en dan gezamenlijk een elliptisch 
stelsel zouden vormen. De numerieke ex- 
perimenten hebben het inmiddels op ver- 
scheidene gronden aannemelijk gemaakt 
dat alle elliptische stelsels (behalve mis- 
schien de zgn. dwerg-elliptische stelsels) 
het resultaat zijn van invangst® ). De sfe- 
rische component (SC, de ‘bulge’) van 
spiraalstelsels en SO’s vertonen veel over- 
eenkomsten met elliptische stelsels. Re- 
cente waarnemingen hebben aangetoond 
dat de grootte van de SC gecorreleerd is 
met de morfologische klasse waartoe een 
St behoort (SO-stelsels hebben grotere 
bulges dan spiraalstelsels). Het is moge- 
lijk dat alle heldere St’s zijn gevormd als 
schijfvormige (spiraal-)stelsels met nauwe- 
lijks of geen SC en dat vervolgens, door 
middel van invangst van kleine stelsels, de 
SC’s en daarmee de huidige typen St’s 
ontstonden. Het feit dat invangst (en 
waarschijnlijk het ontstaan van SO- en 
E-stelsels) vooral plaatsvindt tijdens het 
ontstaan (t <t) van clusters vormt dan 
weer een aanwijzing dat de hierboven ge- 
noemde clustertypen verschillende evolu- 
tiestadia representeren. 

De cD-stelsels in de centra van regelma- 
tige, spiraal-arme clusters lijken hun enor- 
me afmetingen te danken te hebben aan 
het feit dat ze een stuk of tien St’s heb- 
ben opgeslokt. Door het grote aantal bot- 
singen zijn de meeste stelsels in clusters 
het grootste deel van hun halo’s kwijt. Ze 
bewegen nu door een zee van donker halo- 
materiaal. Door hun wisselwerking met 
deze materie verliezen vooral de grote 
stelsels voortdurend baanenergie. Dit pro- 
ces staat bekend onder de naam ‘“dynam- 
ical friction’ en is het gevolg van de ten- 
dens equipartie van energie te bewerkstel- 
ligen tussen de stelsels en het lichtere 
achtergrondmateriaal. Zo zullen de zware 
stelsels in een rijke cluster de neiging 
hebben naar de bodem van de potentiaal- 
put te spiraleren om daar opgegeten te 
worden door het grote centrale cD-stelsel. 


Intergalactisch gas 

De meerderheid van geïdentificeerde ex- 
tragalactische röntgenbronnen zijn geas- 
socieerd met rijke clusters van sterrenstel- 
sels. Het belangrijkste emissiemechanisme 
is remstraling afkomstig van geïoniseerd 
gas met een temperatuur van 

T 107 à® K en dichtheid — 10° 
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Fig. 2. Twee rijke clusters (in negatief 
afgedrukt) met hun röntgenisofoten. De 
röntgenwaarnemingen zijn onlangs gedaan 
met de HEAO-2 (High Energy Astrophys- 
ical Observatory) satelliet. A1 367 is een 
onregelmatige, spiraalrijke cluster, terwijl 
A85 een centraal geconcentreerde cluster 
is met een hoog percentage aan elliptische 
stelsels. Ook de röntgenwaarnemingen la- 
ten een duidelijk verschil in structuur tus- 
sen de twee clusters zien. 


protonen/cm°. 

Deze temperatuur is ongeveer gelijk aan 
de viriaaltemperatuur: kTvi; = GM il Rei 
(Mei en Roi zijn de massa en straal van de 
cluster). De totale massa van het gas is 
0,1 Mj. Het bestaan van dit gas blijkt 
ook uit radiowaarnemingen van sterren- 
stelsels in clusters. Radiostelsels zenden in 
twee tegenoverliggende richtingen wolken 
uit welke een mengsel van relativistische 
deeltjes en magneetvelden bevatten en 
daardoor zichtbaar worden in het radio- 
gebied (synchrotronstraling!). Wanneer 
deze stelsels zich met hoge snelheid door 
het intergalactisch gas bewegen worden 
deze wolken naar achteren gedrukt (‘ram 
pressure’) en we zien een zgn. ‘kop-staart- 
stelsel’ of ‘radio-spoor’. De druk die het 
gas uitoefent langs de bewegingsrichting 
van de sterrenstelsels kan het gas uit spi- 
raal stelsels ‘wegvegen’. Zo’n schoonge- 
veegd stelsel zal als een SO-stelsel worden 
geclassificeerd. 

Er zijn aanwijzingen dat dit inderdaad ge- 
beurt in de centra van rijke clusters. On- 
langs is de aanwezigheid van een emissie- 
lijn afkomstig van 1 5x geïoniseerd Fe ont- 


dekt in het röntgenspectrum van het 
clustergas. De Fe-abundantie die hieruit 
is af te leiden geeft aan dat het gas groten- 
deels van stellaire herkomst is (elementen 
zwaarder dan He zijn in sterren geprodu- 
ceerd, niet in de Big Bang). Het is goed 
mogelijk dat de halo’s van St’s hierbij een 
belangrijke rol hebben gespeeld. Het ga- 
lactische (halo)gas kan door de processen 
die samenhangen met de evolutie van de 
cluster, zoals instorting van de cluster, 
botsingen en ram pressure, uit de stelsels 
zijn vrijgemaakt en verhit tot zijn huidige 
temperatuur. 


Waarnemingen nu en in de nabije 
toekomst 

Het in november 1978 gelanceerde Ein- 
stein-observatorium of HEAO-2 (High 
Energy Astrophysical Observatory) is in 
staat röntgenwaarnemingen te doen met 
een hoog ruimtelijk en spectraal schei- 
dend vermogen (zie figuur 2). De verwach- 
ting is dat deze röntgentelescoop indivi- 
duele clusters zal kunnen onderscheiden 
tot aan een roodverschuiving van 

ze (1 — 3). De evolutie van het clus- 
tergas kan hiermee direct worden waarge- 
nomen. Dit is van belang voor het ant- 
woord op de vraag welke het eerst ont- 
stonden: sterrenstelsels of clusters van 
sterrenstelsels, en hóe ze ontstonden. 


Recente optische waarnemingen van verre 
clusters (z = 1) laten zien dat een groot 
aantal sterrenstelsels in die clusters rela- 
tief erg blauw zijn. In de centrale delen 
van nabije, rijke clusters komen overwe- 
gend E-stelsels voor, welke een rodere 
kleur hebben. De kleur van de stelsels in 
de clusters bij z > 1 is veel meer 
karakteristiek voor spiraalstelsels. Ook ra- 
diostelsels lijken een sterke kleurevolutie 
te ondergaan. Het ligt voor de hand dat er 
een verband is tussen de evolutie van ster- 
renstelsels en de sterke z-afhankelijkheid 
die wordt waargenomen in het aantal 
radiobronnen en quasars in het heelal. 


In 1983 wordt de Space Telescope gelan- 
ceerd. Met dit instrument zullen sterren- 
stelsels tot aan z -” 1 goed bestudeerd kun- 
nen worden over een breed gedeelte van 
hun spectrum (UV tot aan het nabije IR). 
Het ziet ernaar uit dat we binnen afzien- 
bare tijd de vorming van sterrenstelsels en 
clusters van sterrenstelsels direct kunnen 
waarnemen. | 
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